The finite form of Taylor's theorem with the first $n$ terms iswhere $\xi$ lies strictly between $a$ and $x$ for $x\ne a$. A sufficient hypothesis is that $f$ has continuous derivatives through order $n-1$ on the interval joining the endpoints and an $n$th derivative in its interior. The displayed polynomial has degree at most $n-1$; its Taylor remainder uses the $n$th derivative, not the $(n+1)$st. At $x=a$ the remainder is zero. This is a finite Taylor polynomial identity and does not require an infinite Taylor series to converge.

The gradient is orthogonal to each regular level set and points towards increasing $f$. Near either minimum, the linear approximation is $(8\delta x,\frac{7}{3}\delta y)$, so arrows point away from the point. Near the maximum at the origin it is $(-4\delta x,-\frac{5}{3}\delta y)$, so arrows point inwards. Near either saddle point it is $(-2\delta x,\frac{24}{7}\delta y)$: arrows point towards the vertical axis horizontally and away from the saddle vertically.

On the horizontal axis the gradient is $(4x(x^2-1),0)$. Its arrows point left on $x<-1$, right on $-1<x<0$, left on $0<x<1$, and right on $x>1$. On the vertical axis it isConsequently its arrows point down on $y<-1/\sqrt{2}$, up on $-1/\sqrt{2}<y<0$, down on $0<y<1/\sqrt{2}$, and up on $y>1/\sqrt{2}$. These directions, and representative arrows near each stationary point, appear in the right-hand panel. The gradient itself is zero at all five stationary points.

The contour plot is symmetric under either $x\mapsto -x$ or $y\mapsto -y$, because $f$ depends only on their squares. Denote the three critical heights by $m=\frac{7}{2}\ln 3$, $h=\frac{1}{4}+\frac{7}{2}\ln (7/2)$ and $M=1+\frac{7}{2}\ln 3$, where $m<h<M$.

For $m<c<h$, the level set $f=c$ consists of separate closed curves around the two minima. At $c=h$ the branches meet at the two saddle points. For $h<c<M$, there is an inner closed curve surrounding the local maximum and an outer closed curve surrounding all five stationary points. At $c=M$ the inner curve shrinks to the origin; for $c>M$ only the outer closed curve remains. At each saddle point the critical level set crosses in the two directions determined byThe figure plots representative level sets, including the critical saddle height; its second panel adds the gradient directions requested below.

The Hessian matrix has entriesAt each stationary point the mixed entry vanishes, so the signs of the diagonal entries directly classify the local quadratic change:In particular the two local minima are also global minima: the squared term is nonnegative and $\ln (3+y^2)\ge \ln 3$, with equality in both bounds exactly at $(\pm 1,0)$. The local maximum at the origin is not a global maximum, since $f$ grows without bound as $x^2+y^2\to \infty$. The opposite diagonal signs at the other two points establish their saddle point character.

Put $s=x^2+y^2$. The gradient components areA stationary point requires both components to vanish. If $x\ne 0$, then $s=1$; the second component subsequently forces $y=0$, giving $(\pm 1,0)$. If $x=0$, either $y=0$ or, writing $u=y^2$,Only $u=1/2$ is nonnegative. Thus there are five stationary points, with valuesThe cases $x=0$ and $s=1$ exhaust the first partial derivative equation, so no other stationary points have been lost by division.

The volume element is the absolute Jacobian determinant times the coordinate differentials. In the Cartesian coordinate system, cylindrical coordinates $(r,\theta ,z)$ and spherical coordinates $(\rho ,\vartheta ,\phi )$, respectively,Here $\rho$ is distance from the origin, $\vartheta$ is the polar angle from the positive $z$ axis, and $\phi$ is the azimuthal angle. The spherical coordinates convention fixes which angle enters the sine factor.

For polar coordinates, $x=r\cos \theta$ and $y=r\sin \theta$. The Jacobian determinant is $r$, so the area element isThis accounts for the increasing length of an angular strip as its distance from the origin increases.