Skip to content
Part IA Michaelmas, Lent, Easter Term

Stationary Points

Definition

A stationary point of f(x,y)f(x, y) occurs where:

f=0\nabla f = \mathbf{0}

That is: fx=0\frac{\partial f}{\partial x} = 0 and fy=0\frac{\partial f}{\partial y} = 0.

At a stationary point, the function is “flat” in all directions.


Types of Stationary Points

  1. Local minimum: ff increases in all directions
  2. Local maximum: ff decreases in all directions
  3. Saddle point: ff increases in some directions, decreases in others

3D surface showing three types: local minimum (bowl), local maximum (hill), and saddle point (hyperbolic paraboloid)


The Hessian Matrix

At a stationary point, the second-order behaviour is determined by the Hessian matrix:

H=(fxxfxyfyxfyy)H = \begin{pmatrix} f_{xx} & f_{xy} \\ f_{yx} & f_{yy} \end{pmatrix}

Since fxy=fyxf_{xy} = f_{yx} (for smooth functions), HH is symmetric.


Classification Criterion

Let D=det(H)=fxxfyy(fxy)2D = \det(H) = f_{xx}f_{yy} - (f_{xy})^2

ConditionClassification
D>0D > 0 and fxx>0f_{xx} > 0Local minimum
D>0D > 0 and fxx<0f_{xx} < 0Local maximum
D<0D < 0Saddle point
D=0D = 0Inconclusive

Why D>0D > 0 implies extremum: Both eigenvalues of HH have same sign.

Why D<0D < 0 implies saddle: Eigenvalues have opposite signs.

Why D=0D = 0 is inconclusive: One eigenvalue is zero; higher-order terms needed.


Taylor Series Interpretation

Near stationary point (a,b)(a, b):

f(a+h,b+k)f(a,b)+12(h2fxx+2hkfxy+k2fyy)f(a+h, b+k) \approx f(a,b) + \frac{1}{2}(h^2 f_{xx} + 2hk f_{xy} + k^2 f_{yy})

In matrix form: ff(a,b)+12hTHhf \approx f(a,b) + \frac{1}{2}\mathbf{h}^T H \mathbf{h}

The quadratic form hTHh\mathbf{h}^T H \mathbf{h} determines whether we have:

  • Minimum: quadratic positive in all directions (HH positive definite)
  • Maximum: quadratic negative in all directions (HH negative definite)
  • Saddle: quadratic can be positive or negative depending on direction

Example

Find and classify stationary points of f(x,y)=x34x2+2xyy2f(x, y) = x^3 - 4x^2 + 2xy - y^2.

Find stationary points:

  • fx=3x28x+2y=0f_x = 3x^2 - 8x + 2y = 0
  • fy=2x2y=0f_y = 2x - 2y = 0

From fyf_y: y=xy = x.

Substituting into fxf_x: 3x28x+2x=03x^2 - 8x + 2x = 0 3x26x=03x^2 - 6x = 0 x(3x6)=0x(3x - 6) = 0

So x=0x = 0 or x=2x = 2.

Stationary points: (0,0)(0, 0) and (2,2)(2, 2).

Compute Hessian:

  • fxx=6x8f_{xx} = 6x - 8
  • fxy=2f_{xy} = 2
  • fyy=2f_{yy} = -2

At (0,0)(0, 0): H=(8222)H = \begin{pmatrix} -8 & 2 \\ 2 & -2 \end{pmatrix}

D=(8)(2)4=12>0D = (-8)(-2) - 4 = 12 > 0, and fxx=8<0f_{xx} = -8 < 0.

Local maximum.

At (2,2)(2, 2): H=(4222)H = \begin{pmatrix} 4 & 2 \\ 2 & -2 \end{pmatrix}

D=(4)(2)4=12<0D = (4)(-2) - 4 = -12 < 0.

Saddle point.


Global vs Local

  • Local minimum: f(x,y)f(a,b)f(x,y) \geq f(a,b) for all nearby (x,y)(x,y)
  • Global minimum: f(x,y)f(a,b)f(x,y) \geq f(a,b) for all (x,y)(x,y) in the domain

Multiple stationary points may exist; the global minimum/maximum is found by comparing function values.


Connection to Eigenvalues

The Hessian is real symmetric, so has real eigenvalues λ1,λ2\lambda_1, \lambda_2.

  • Both positive → minimum (D>0,fxx>0D > 0, f_{xx} > 0)
  • Both negative → maximum (D>0,fxx<0D > 0, f_{xx} < 0)
  • One positive, one negative → saddle (D<0D < 0)

D=λ1λ2D = \lambda_1\lambda_2 and fxxf_{xx} determines the sign of λ1\lambda_1 (for convention).


Why Stationary Points Matter

  • Optimisation problems (minimise cost, maximise efficiency)
  • Machine learning (gradient descent finds minima of loss functions)
  • Physics (equilibrium states are stationary points of potential energy)
  • Economics (profit maximisation, cost minimisation)