Solution (source code)

= Solution

The <gradient> equations are $y=x^3$ and $x=y^3$. Their real solutions are $(0,0)$, $(1,1)$ and $(-1,-1)$. The <Hessian matrix> is $\begin{pmatrix}-12x^2&4\\4&-12y^2\end{pmatrix}$, so the origin has <eigenvalues> $\pm4$ and is a <saddle point>, with value zero. At either other point the <eigenvalues> are $-8,-16$, giving strict maxima with value two. They are global maxima because
$$
4xy-x^4-y^4\le4xy-2x^2y^2=2-2(xy-1)^2\le2,
$$
and equality occurs precisely at those two points.

In <polar coordinates>, the zero contour is
$$
0=2r^2\sin2\theta-r^4(\cos^4\theta+\sin^4\theta).
$$
Besides the saddle $r=0$, its nonzero branches satisfy
$$
\boxed{r^2=\frac{2\sin2\theta}{\cos^4\theta+\sin^4\theta}
=\frac{8\sin2\theta}{3+\cos4\theta}.}
$$
The denominator is positive; real radii require $\sin2\theta\ge0$. The branches form two closed lobes in the first and third quadrants, meeting at the origin. Positive contours lie inside these lobes around the maxima, and negative contours extend outside. Along the coordinate axes the function is negative away from the saddle, whereas along $y=x$ it is initially positive; these signs show the crossing geometry in the sketch.