We need to back up a little, because there's an error in the critical point calculation. We have
$$ \frac{\partial f}{\partial x} \ = \ 4x^3 \ - \ 16x \ = \ 4x \ (x^2 \ - \ 4 ) \ = \ 0 \ \ , $$
$$ \frac{\partial f}{\partial y} \ = \ 4y^3 \ - \ 36y \ = \ 4y \ (y^2 \ - \ 9 ) \ = \ 0 \ \ . $$
The possible "pairings" of coordinate solutions for these equations give us nine critical points: $ \ (0, 0) \ , \ (0, \ \pm 3) \ , \ (\pm 2, \ 0) \ , \ (\pm 2, \ \pm 3) \ $ .
There is also an error in the Hessian matrix, since one of the first partial derivatives was incorrect; the second partial derivatives are
$$ \frac{\partial^2 f}{\partial \ x^2} \ = \ 12x^2 \ - \ 16 \ \ \ , \ \ \frac{\partial^2 f}{\partial \ y^2} \ = \ 12y^2 \ - \ 36 \ \ \ , \ \ \frac{\partial^2 f}{\partial x \ \partial y} \ = \ \frac{\partial^2 f}{\partial y \ \partial x} \ = \ 0 \ \ . $$
So the Hessian is
$$ \mathbf{H} \ = \ \begin{pmatrix} 12x^2-16 & 0 \\0&12y^2-36 \\ \end{pmatrix} \ \ . $$
Because the entries in the matrix are even functions, it saves us some work to evaluate it for the critical points:
(0, 0) -- $ \begin{pmatrix} -16 & 0 \\0 &-36 \\ \end{pmatrix} \ \ , $ which is negative definite, making this location a local maximum of the function;
(0, ±3) -- $ \begin{pmatrix} -16 & 0 \\0 &72 \\ \end{pmatrix} \ \ , $ which is indefinite;
(±2, 0) -- $ \begin{pmatrix} 32 & 0 \\0 &-36 \\ \end{pmatrix} \ \ , $ which is also indefinite;
(±2, ±3) -- $ \begin{pmatrix} 32 & 0 \\0 &72 \\ \end{pmatrix} \ \ , $ which is positive definite, hence these points are local minima.
We can verify this using the $ \ D-$ index often learned on our first encounter with multivariate optimization, $ \ D \ = \ f_{xx} \ f_{yy} \ - \ f^2_{xy} \ $ :
$$ D \ \vert_{(0,0)} \ = \ (-16) \ (-36) \ - \ 0^2 \ > \ 0 \ \ , \ \ f_{xx} \ < \ 0 \ \ , $$
so this is a local maximum;
$$ D \ \vert_{(0,\pm 3)} \ = \ (-16) \ (72) \ - \ 0^2 \ < \ 0 \ \ , \ \ D \ \vert_{(\pm 2,0)} \ = \ (32) \ (-36) \ - \ 0^2 \ < \ 0 \ \ , $$
which indicates that these are "saddle points"; and
$$ D \ \vert_{(\pm 2, \pm 3)} \ = \ (32) \ (72) \ - \ 0^2 \ > \ 0 \ \ , \ \ f_{xx} \ > \ 0 \ \ , $$
hence these are local minima.
A three-dimensional graph of $ \ f(x,y) \ $ is a bit difficult to read because of the large range of relevant values of the function (as we shall see), and because of the complicated character of part of the surface. Instead, we will view "slices" through the surface to investigate the critical points.
A slice through $ \ y \ = \ 0 \ $ shows the local maximum at $ \ ( 0, \ 0, \ 0) \ $ and what appear to be minima at $ \ x \ = \ \pm 2 \ $ . However, if we take slices nearby at $ \ y \ = \ \pm \frac{1}{2} \ , $ we see that the cross-sectional curve "shifts downward" (computing the function at these values of $ \ y \ $ confirms this). So $ \ (0, \ 0) \ $ indeed marks a local maximum, but $ \ (0, \ \pm 2) \ $ are saddle points (maxima in the $ \ y-$ direction, but minima in the $ \ x-$ direction).
Similarly, we can examine the rest of the critical points by using slices in the vicinity of $ \ y \ = \ \pm 3 \ $ . The cross-sectional curves at $ \ y \ = \ \pm \frac{5}{2} \ $ and $ \ y \ = \ \pm \frac{7}{2} \ $ "shift upward" , so we find that the points $ \ (\pm 2, \ \pm 3) \ $ are true local minima, but the points $ \ (0, \ \pm 3) \ $ are saddle points (maxima in the $ \ x-$ direction, minima in the $ \ y-$ direction).
For questions (c) and (d), then, we have identified one local maximum at $ \ (0, \ 0) $ and four local minima at $ \ (\pm 2, \ \pm 3) \ $ . The degree of the polynomial in each "dimension" of our function is even, and the leading coefficient is positive in both cases ( $ \ x^4 \ - \ 8x^2 \ $ and $ \ y^4 \ - \ 18y^2 \ $ ) , so the surface "opens upward" in both dimensions. We can conclude that the local minimal are in fact global minima; hence, the global minimum value of our function is $ \ f(\pm 2, \ \pm 3) \ $ $= \ 16 \ - \ 32 \ + \ 81 \ - \ 162 \ = \ -97 \ $ . On the other hand, these polynomials have no upper bound, so our function grows without limit when we "get far from the origin". Thus, the local maximum we found at the origin is only that; the function has no global maximum.
One way to understand how the test works is by looking at the Taylor Series of the function $f(x)$ centered around the critical point, $x = c$:
$$
f(x) = f(c) + f'(c)(x-c) + \frac{f''(c)}{2}(x-c)^2 + \cdots
$$
Note: In your question you said that the n-th derivative is non-zero. Here I'm assuming the n+1-st derivative is the first to be non-zero at $x=c$. It doesn't make a difference, it's just the way I learned it.
If $f'(c) = \cdots = f^{(n)}(c) = 0$ and $f^{(n+1)} \ne 0$, then the Taylor Series ends up looking like this:
$$
f(x) = f(c) + \frac{f^{(n+1)}(c)}{(n+1)!}(x-c)^{n+1} + \frac{f^{(n+2)}(c)}{(n+2)!}(x-c)^{n+2} + \cdots
$$
Consider what happens when you move $f(c)$ to the other side of the equation:
$$
f(x) - f(c) = \frac{f^{(n+1)}(c)}{(n+1)!}(x-c)^{n+1} + \frac{f^{(n+2)}(c)}{(n+2)!}(x-c)^{n+2} + \cdots
$$
What does $f(x) - f(c)$ mean?
- If $f(x) - f(c) = 0$, then $f(x)$ has the same value as it does at $x = c$.
- If $f(x) - f(c) < 0$, then $f(x)$ has a value less than it has at $x = c$.
- If $f(x) - f(c) > 0$, then $f(x)$ has a value greater than it has at $x = c$.
We expect $f(x) - f(c) = 0$ at $x = c$ (the equation reflects this), but we're more interested in what it does on either side of $x = c$. When $x$ is really close to $c$, i.e. $(x-c)$ is a really small number, we can say:
$$
f(x) - f(c) \approx \frac{f^{(n+1)}(c)}{(n+1)!}(x-c)^{n+1}
$$
because the higher powers of a small number "don't matter" as much.
Concerning local extrema
If $n$ is odd, then our approximation of $f(x) - f(c)$ is an even-power polynomial. That means $f(x)$ has the same behavior - is either less than or greater than $f(c)$ - on both sides of $x = c$. Therefore it's a local extreme. If $f^{(n+1)}(c) > 0$, then $f(x)$ is greater than $f(c)$ on both sides of $x = c$. Otherwise, if $f^{(n+1)}(c) < 0$, then $f(x)$ is less than $f(c)$ on both sides of $x = c$
If, on the other hand, $n$ is even, then our approximation of $f(x) - f(c)$ is an odd-power polynomial centered around $x = c$. Therefore $f(x)$ will be greater than $f(c)$ on one side of $x = c$, and less on the other. That means $x = c$ isn't a local extreme.
Concerning saddle points
Note that if you differentiate both sides of our approximation twice, you get:
$$
f''(x) \approx \frac{f^{(n+1)}(c)}{(n-1)!}(x-c)^{n-1}
$$
If $n$ is even, this is another odd-power polynomial centered around $x = c$. It therefore has opposite behavior on each side of $x = c$, giving you a saddle point.
Best Answer
This should bring you forward:
However, you should - at least once - really understand those criteria (e.g. graphically) rather than just apply them.