Multivariable calculus

Critical Points of xy e^(-x^2-y^2)

f(x,y) = xy e^(-x^2-y^2) has exactly five critical points. The origin is a saddle with D = -1. The two points where x and y are both 1/sqrt2 or both -1/sqrt2 are local maxima with D = 4/e^2, about 0.5413, and value 1/(2e). The two points where x and y have opposite signs are local minima with the same D.

f(x,y)=xyex2y2f(x,y) = xy\,e^{-x^{2}-y^{2}}
  • (0,0)\left(0, 0\right) saddle pointdiscriminant D = -1
  • (0.7071067811865475,0.7071067811865475)\left(0.7071067811865475, 0.7071067811865475\right) local maximumdiscriminant D = 0.5413411329464508
  • (0.7071067811865475,0.7071067811865475)\left(-0.7071067811865475, -0.7071067811865475\right) local maximumdiscriminant D = 0.5413411329464508
  • (0.7071067811865475,0.7071067811865475)\left(0.7071067811865475, -0.7071067811865475\right) local minimumdiscriminant D = 0.5413411329464508
  • (0.7071067811865475,0.7071067811865475)\left(-0.7071067811865475, 0.7071067811865475\right) local minimumdiscriminant D = 0.5413411329464508

Factor first, then enumerate the cases

Each partial needs the product rule, and the two resulting terms share a common factor that is worth pulling out immediately.

fx=yex2y2(12x2),fy=xex2y2(12y2)\frac{\partial f}{\partial x} = y\,e^{-x^{2}-y^{2}}\left(1 - 2x^{2}\right), \qquad \frac{\partial f}{\partial y} = x\,e^{-x^{2}-y^{2}}\left(1 - 2y^{2}\right)

The exponential is never zero, so it divides out of both equations and leaves two factored conditions.

y(12x2)=0andx(12y2)=0y\left(1 - 2x^{2}\right) = 0 \qquad \text{and} \qquad x\left(1 - 2y^{2}\right) = 0

Now enumerate. The first equation gives y=0y = 0 or x=±12x = \pm\frac{1}{\sqrt{2}}. The second gives x=0x = 0 or y=±12y = \pm\frac{1}{\sqrt{2}}. Pairing them off:

  • y=0y = 0 with x=0x = 0 gives the origin.
  • y=0y = 0 with y=±12y = \pm\frac{1}{\sqrt{2}} is a contradiction, so nothing.
  • x=±12x = \pm\frac{1}{\sqrt{2}} with x=0x = 0 is a contradiction, so nothing.
  • x=±12x = \pm\frac{1}{\sqrt{2}} with y=±12y = \pm\frac{1}{\sqrt{2}} gives four points, one in each quadrant.

Five points in total, and because the enumeration was exhaustive, that is the complete list.

Classifying all five at once

The second partials are more pleasant than they look, especially the mixed one, which factors completely.

fxx=2xy(32x2)ex2y2,fyy=2xy(32y2)ex2y2,fxy=(12x2)(12y2)ex2y2f_{xx} = -2xy\left(3 - 2x^{2}\right)e^{-x^{2}-y^{2}}, \quad f_{yy} = -2xy\left(3 - 2y^{2}\right)e^{-x^{2}-y^{2}}, \quad f_{xy} = \left(1-2x^{2}\right)\left(1-2y^{2}\right)e^{-x^{2}-y^{2}}

At the origin the first two vanish and fxy=1f_{xy} = 1, so D=1D = -1 and the point is a saddle. Near (0,0)(0,0) the exponential is close to 1 and ff behaves like xyxy, positive in two opposite quadrants and negative in the other two.

At the four outer points, x2=y2=12x^{2} = y^{2} = \frac{1}{2}, so both factors in fxyf_{xy} are zero and the mixed partial disappears. The exponent is 1-1 there, and 32x2=23 - 2x^{2} = 2, so the diagonal entries reduce to fxx=fyy=4xye1f_{xx} = f_{yy} = -4xy\,e^{-1}.

D=(4xye1)20=16x2y2e2=4e20.5413D = \left(-4xy\,e^{-1}\right)^{2} - 0 = 16x^{2}y^{2}e^{-2} = \frac{4}{e^{2}} \approx 0.5413

The discriminant is the same positive number at all four, so fxxf_{xx} alone decides. Where xx and yy agree in sign, xy=12xy = \frac{1}{2} and fxx=2e<0f_{xx} = -\frac{2}{e} < 0: local maxima with f=12e0.184f = \frac{1}{2e} \approx 0.184. Where they differ in sign, xy=12xy = -\frac{1}{2}, fxx=2e>0f_{xx} = \frac{2}{e} > 0, and those are local minima with f=12ef = -\frac{1}{2e}.

Those extremes are global. Since xyx2+y22|xy| \le \frac{x^{2}+y^{2}}{2}, writing r2=x2+y2r^{2} = x^{2}+y^{2} gives fr22er2|f| \le \frac{r^{2}}{2}e^{-r^{2}}, and the right side is largest at r2=1r^{2} = 1, where it equals 12e\frac{1}{2e}.

The mistake students make

The most common loss is the origin itself. Students factor the two equations, jump to the interesting case x=±12x = \pm\frac{1}{\sqrt{2}} and y=±12y = \pm\frac{1}{\sqrt{2}}, report four critical points and stop. The case x=y=0x = y = 0 satisfies both equations just as legitimately, and it is the only saddle on the surface. Enumerating all four pairings on paper, including the ones that die, is what prevents this.

The second error is dropping a term from the product rule. Differentiating xyex2y2xy\,e^{-x^{2}-y^{2}} in xx means differentiating a product of xyxy and the exponential, so both pieces contribute.

fx=yex2y2+xy(2x)ex2y2=yex2y2(12x2)\frac{\partial f}{\partial x} = y\,e^{-x^{2}-y^{2}} + xy\left(-2x\right)e^{-x^{2}-y^{2}} = y\,e^{-x^{2}-y^{2}}\left(1 - 2x^{2}\right)

Keep only the second term and you get 2x2yex2y2-2x^{2}y\,e^{-x^{2}-y^{2}}, with 2xy2ex2y2-2xy^{2}e^{-x^{2}-y^{2}} from the same slip in yy. Both vanish on the whole xx-axis and the whole yy-axis, so the answer becomes two full lines of critical points in place of the single saddle at the origin, and the 12x21 - 2x^{2} factor that locates the four real extrema never appears.

A cheap check: the four outer critical points must sit at x2=y2=12x^{2} = y^{2} = \frac{1}{2}, distance 1 from the origin. If your answer puts them anywhere else, or puts them on an axis, the product rule is the first place to look.

Frequently asked questions

Why do all four outer points share the same discriminant?

Because the surface has a symmetry. Replacing xx with x-x turns ff into f-f, and so does replacing yy with y-y. Sign changes flip maxima into minima but leave curvature magnitudes alone, so fxxf_{xx} and fyyf_{yy} change sign together while fxyf_{xy} stays zero. The product fxxfyyf_{xx}f_{yy} is then unchanged, and D=4e2D = \frac{4}{e^{2}} at every one of the four.

Is the list of five points really complete?

Yes. Both partials factor as a product of a never zero exponential with two simple algebraic factors, and setting a product to zero means setting one factor to zero. The four pairings of those factors are exhaustive, two are contradictory, and the remaining two yield the origin and the four outer points. No approximation or search is involved, so nothing can be hiding elsewhere in the plane.