I found the quadratic approximation as $9 + \frac{1}{2}(-9x^2 - 9y^2)$
The problem is that the triple derivatives all end up 0 at (0,0), so I get that the error approximation is 0. Wolfram alpha calculates the triple derivatives having sin(y) or sin(x) in them making them 0. I know you are supposed to plug in the values into $1/6(f_(xxx)(x^3) + ...)$ but each triple derivative ends up as 0 which turns out to be the wrong answer.