To find critical points we need to solve the system of equations
and
. Let
.
Let's take a look at a plot of
.
See Figure BLANK.
Now let's take a look at a contour plot of
. We do it just like the plot of the surface, except we use the argument contour = map.
See Figure BLANK.
Now let's find the first order partial derivatives of
, set them equal to zero, and solve for values of
. We do this with the solve function, which assumes that the expressions are set equal to zero by default. Keep in mind that solve finds all real and complex solutions; we only care about the real valued solutions, however, and will ignore the rest.
Critical points are at
,
, and
. (Note that both partial derivatives exist everywhere.) We need to see what the Hessian says at those locations.
We do the Second Derivative Test by plugging in the above three points into
and
. (Of course we can do it mentally but let's try it with Maxima.) A quick way to plug numbers into expressions is with the subst function.
From the above we conclude that
is a local minimum of
, and the value is
Doing the same for
and
shows that
is a saddle point and
is also a local minimum (we know this already by symmetry).
G. Jay Kerns 2009-12-01