I'm working on the following problem:
Consider the $f:\mathbb{R}^2 \mapsto \mathbb{R}$ defined by $$f(x,y) = \frac{xy(x^2-y^2)}{x^2+y^2}$$ for everywhere outside the origin and $f(0,0)=0$ at the origin is. Show that the second derivative of $f$ exists but the mixed partials are not equal at the origin.
My work:
I'm getting for $\nabla f = <\frac{yx^4+4x^2y^3-y^5}{(x^2+y^2)^2},\frac{x^5-4x^3y^2-xy^4}{(x^2+y^2)^2}>$. My analysis show that the partials are continuous at zero and take the value $0$ at the origin. I'm having trouble analyzing the second derivative:
$$ A = \begin{bmatrix}f_{xx}&f_{xy}\\f_{yx}&f_{yy}\end{bmatrix} = \begin{bmatrix}\frac{-4xy^3(x^2-3y^2)}{(x^2+y^2)^3}&\frac{x^6+9x^4y^2-9x^2y^4-y^6}{(x^2+y^2)^3}\\\frac{x^6+9x^4y^2-9x^2y^4-y^6}{(x^2+y^2)^3}&\frac{4yx^3(y^2-3x^2)}{(x^2+y^2)^3}\end{bmatrix} $$
$$\lim_{\|h\|\rightarrow 0} \frac{\nabla f((0,0)+h)-\nabla f((0,0))- A \cdot h\|}{\|h\|} \\ = \lim_{\|h\|\rightarrow 0} \frac{\|\nabla f(h)-\langle f_{xx}h_x+f_{xy}h_y,f_{yx}h_x+f_{yy}h_y\rangle\|}{\|h\|} \\= \lim_{\|h\|\rightarrow 0} \frac{\|\nabla f(h)-\langle\frac{-4xy^3(x^2-3y^2)}{(x^2+y^2)^3}h_x+\frac{x^6+9x^4y^2-9x^2y^4-y^6}{(x^2+y^2)^3}h_y,\frac{x^6+9x^4y^2-9x^2y^4-y^6}{(x^2+y^2)^3}h_x+\frac{4yx^3(y^2-3x^2)}{(x^2+y^2)^3}h_y\rangle\|}{\|h\|}$$
When I do a polar conversion on this thing and send $r \rightarrow 0$ (uniform for all $\theta$), I get all my $r$'s canceling out -- leaving a limit dependent on theta. But according to the problem, this ratio should go to zero though. Does anyone know what's going on?
I've seen other problems like this that show partials are not continuous at the the origin. That's different though from just establishing existance of the total derivative, right?