You probably mean the function to be

We use the second (partial) derivative test for this.
Compute the first partial derivatives of f :


Find the critical points of f, where both first derivatives vanish. Since
for all x and y, we're left with

which means either x = 0 or x² + 2y² = 1, and

so y = 0 or x² + 2y² = 2.
• If x = 0 and y = 0, then one critical point is (0, 0).
• If x ≠ 0 and y = 0, then x² + 2y² = x² = 1, which means x = ±1. So we have two additional critical points, (-1, 0) and (0, 0).
• If x = 0 and y ≠ 0, then x² + 2y² = 2y² = 2, so y² = 1 and y = ±1, and we get two more critical points at (0, -1) and (0, 1).
All five of these critical points belong to D.
Compute the Hessian matrix for f and check the sign of its determinant at each critical point. Remember what the second partial derivative test concludes:
• If det H(x, y) is negative, then that point is a saddle point
• If det H(x, y) is positive, and ∂²f/∂x² is negative, then that point is a maximum
• If det H(x, y) is positive, and ∂²f/∂x² is positive, then that point is a minimum
The Hessian matrix for f is


Note that ∂²f/∂x² is the (1, 1) entry. Compute the determinant and evaluate it at each critical point.


⇒ local minimum at (0, 0) of f(0, 0) = 0

⇒ saddle point (-1, 0)

⇒ saddle point at (1, 0)


⇒ local maximum at (0, -1) of f(0, -1) = 2/e ≈ 0.736


⇒ local maximum at (0, 1) of f(0, 1) = 2/e ≈ 0.736
Check for extrema along the boundary. We can parameterize the boundary by setting
x(t) = 2 cos(t)
y(t) = 2 sin(t)
with 0 ≤ t < 2π. Then f(x, y) = f(x(t), y(t)) reduces to a function of only t, which we call F(t) :


Find the critical points of F on the interval [0, 2π) :

sin(2t) = 0
2t = arcsin(0) + 2nπ or 2t = π - arcsin(0) + 2nπ
(where n is any integer)
2t = 2nπ or 2t = π + 2nπ
t = nπ or t = π/2 + nπ
Over [0, 2π), critical points occur at t = 0, t = π/2, t = π, and t = 3π/2.
Check the sign of the second derivative of F at each critical point. The test in the one-variable case says
• if F''(t) < 0, then that point is a maximum
• if F''(t) > 0, then that point is a minimum


⇒ local minimum when t = 0 of F(0) = 4/e⁴ ≈ 0.0733. This t corresponds to the point (2, 0), since x(0) = 2 and y(0) = 0.

⇒ local maximum when t = π/2 of F(π/2) = 8/e⁴ ≈ 0.147, referring to the point (0, 2)

⇒ local minimum when t = π of F(π) = 4/e⁴ ≈ 0.0733 at the point (-2, 0)

⇒ local maximum when t = 3π/2 of F(3π/2) = 8/e⁴ ≈ 0.147 at (0, -2)
So, the absolute extrema of f(x, y) over D are
• minimum : f(0, 0) = 0
• maximum : f(0, -1) = f(0, 1) = 2/e ≈ 0.736