 Solution

Set both first partial derivatives equal to zero:Thus $x=0,\pm 1/\sqrt{2}$ and $y=0,\pm \sqrt{2}$, independently, giving exactly nine stationary points. The Hessian matrix is diagonal everywhere:Its diagonal entries are its eigenvalues, so the second-derivative test immediately classifies the points:

In particular,All nine stationary points are nondegenerate, so no higher-order test is needed. Neither type of local extremum is global: $f(x,0)=x^4-x^2\to +\infty$ as $|x|\to \infty$, whereas $f(0,y)=-y^4+4y^2\to -\infty$ as $|y|\to \infty$.

To organize the contour plot, complete the squares:The level set through all four outer saddle points of a scalar function factors exactly:It is therefore the union of the circle $x^2+y^2=5/2$ and the hyperbola $y^2-x^2=3/2$, which cross at those four saddle points of a scalar function. Near the origin, $f=-x^2+4y^2+O((x^2+y^2{)}^2)$, so the zero contours have tangents $y=\pm x/2$. These two exact structural features, together with the nested contours around the four local extrema, determine the requested sketch. They are an example of a separable quartic landscape with nine stationary points.