Characteristic Equation Method for Linear ODEs

Core Idea

For constant-coefficient linear ODEs, assume a solution y = e^(rx) and substitute to obtain a characteristic equation. For y'' + py' + qy = 0, the characteristic equation is r² + pr + q = 0. The roots r determine the solution form: distinct real roots give y = c₁e^(r₁x) + c₂e^(r₂x); complex roots give oscillatory solutions; repeated roots require x factors. This algebraic approach elegantly solves a wide class of equations.

Explainer

You know from second-order linear ODEs that the general solution to y″ + py′ + qy = 0 is a linear combination of two independent solutions. The characteristic equation method is the systematic algorithm for finding those two solutions when p and q are constants. The key insight is an inspired guess, or ansatz: if we try y = e^(rx), then y′ = re^(rx) and y″ = r²e^(rx). Substituting into the equation gives r²e^(rx) + pre^(rx) + qe^(rx) = 0. Factor out e^(rx) — which is never zero — and you get r² + pr + q = 0. The differential equation has become a quadratic.

This quadratic, called the characteristic equation, is solved with the quadratic formula you already know. The roots r₁ and r₂ determine the solution form, and there are three cases based on the discriminant p² − 4q:

• Two distinct real roots (p² − 4q > 0): The two independent solutions are e^(r₁x) and e^(r₂x), giving y = c₁e^(r₁x) + c₂e^(r₂x). Each root contributes its own exponential.
• **Complex conjugate roots** r = α ± βi (p² − 4q < 0): The complex exponentials e^((α+βi)x) and e^((α−βi)x) are valid, but we can use Euler's formula to rewrite them in real form: y = e^(αx)(c₁cos(βx) + c₂sin(βx)). This is the oscillatory case — the solution is a sinusoid with exponentially growing or decaying amplitude.
• Repeated root r₁ = r₂ = r (p² − 4q = 0): Two copies of e^(rx) are not independent, so e^(rx) alone cannot span the solution space. The second independent solution is xe^(rx), giving y = (c₁ + c₂x)e^(rx).

The elegance of this method is that it converts a calculus problem into a purely algebraic one. The structure of the ODE's solutions — exponential growth or decay, oscillation, polynomial growth — is entirely encoded in the location of the characteristic roots in the complex plane. Roots with large negative real parts decay fast; purely imaginary roots oscillate without damping; roots with positive real parts grow without bound. This geometric picture of roots predicting behavior is the foundation for understanding stability in differential equations and control systems.