Clearly (0,0) is a minima. From the graph of the function, we
suspect that there are maxima along the 45-degree diagonals
in the (x,y)-plane. If we parametrize the diagonal by (t,t), then
solve for the critical points along that line, we find
that there are critical points at
(sqrt(5)/2, sqrt(5)/2)) and (-sqrt(5)/2, -sqrt(5)/2)).
The other diagonal is parametrized by (t, -t) and along this line the
critical points are (sqrt(3)/2, sqrt(3)/2) and
(-sqrt(3)/2, -sqrt(3)/2).
Non-symmetric solutions (saddle points) occur somewhere near
(0.259, -0.966), (-0.259, 0.966),
(-0.966, 0.259), (0.966, -0.259).