Linear Second-Order ODEs

One-line statement: A linear second-order ODE has the form

a (x) y^{''} + b (x) y^{'} + c (x) y = f (x)

, where the unknown function and its derivatives appear only to the first power and are not multiplied together. Key result in words: The general solution is the sum of the general solution of the homogeneous equation (a two-dimensional vector space spanned by two linearly independent solutions) plus any single particular solution of the full equation. Where it applies: Mechanical vibrations, electrical circuits, quantum mechanics, heat conduction, wave propagation, beam deflection, and virtually every linear physical system with inertia. Who/when: Foundations laid by Euler, d'Alembert, and Lagrange (18th century); existence theory by Cauchy, Lipschitz, and Picard (19th century).

A linear second-order ordinary differential equation is an equation of the form $a (x) y^{''} + b (x) y^{'} + c (x) y = f (x)$ in which the highest derivative is the second, and the dependent variable $y$ together with its derivatives enters linearly. These equations are the workhorses of applied mathematics because Newton's second law, $F = ma$ , produces a second derivative whenever forces depend on position and velocity. The structure of their solution set is remarkably rigid: solutions of the homogeneous version form a vector space of exactly dimension two, so just two well-chosen solutions encode all others. A surprising fact is that even though most such equations cannot be solved in elementary closed form, the qualitative behavior — oscillation, decay, resonance — is completely dictated by the coefficients, allowing precise prediction without explicit formulas.

The Core Insight

Linearity converts a differential equation into linear algebra. Because the operator

L [y] = a y^{''} + b y^{'} + cy

satisfies

L [α y_{1} + β y_{2}] = α L [y_{1}] + β L [y_{2}]

, the homogeneous solution set is the kernel of a linear map and is therefore a vector space. For a second-order equation that kernel is exactly two-dimensional. Everything else — superposition, the role of the Wronskian, the general solution = homogeneous + particular structure — is a direct consequence of this single algebraic fact.

Formal Statement

Linear Second-Order ODE

A linear second-order ODE on an interval $I$ is any equation expressible as

a (x) y^{''} + b (x) y^{'} + c (x) y = f (x),

⟳This is the most general form in which

y, y^{'}, y^{''}

appear to the first power with no products among them; coefficients may depend on

x

but not on

y

where:

$y = y (x)$ is the unknown function,
$a (x), b (x), c (x)$ are the coefficient functions (with $a (x) \neq = 0$ on $I$ ),
$f (x)$ is the forcing term (or source/inhomogeneity).

The equation is homogeneous if $f (x) \equiv 0$ , and inhomogeneous (nonhomogeneous) otherwise. Dividing by $a (x)$ gives the standard form

y^{''} + p (x) y^{'} + q (x) y = g (x), p = \frac{b}{a}, q = \frac{c}{a}, g = \frac{f}{a} .

⟳Valid wherever

a (x) \neq = 0

; points where

a

vanishes are singular points requiring separate analysis.

If the coefficients $p, q$ are constants, the equation has constant coefficients; otherwise it has variable coefficients.

Differential Operator and Linearity

Define the linear differential operator

L \equiv \frac{d ^{2}}{d x ^{2}} + p (x) \frac{d}{d x} + q (x), L [y] = y^{''} + p y^{'} + q y .

⟳

L

acts on the vector space

C^{2} (I)

; linearity follows because differentiation is linear and multiplication by a fixed function is linear.

Then $L [α y_{1} + β y_{2}] = α L [y_{1}] + β L [y_{2}]$ for all scalars $α, β$ and twice-differentiable $y_{1}, y_{2}$ .

Existence and Uniqueness (Cauchy–Picard)

If $p (x)$ , $q (x)$ , and $g (x)$ are continuous on an open interval $I$ containing $x_{0}$ , then for any prescribed values $y_{0}, y_{0}^{'}$ the initial value problem

y^{''} + p y^{'} + q y = g, y (x_{0}) = y_{0}, y^{'} (x_{0}) = y_{0}^{'}

⟳Reduce to a first-order system

u^{'} = A (x) u + b (x)

with

u = (y, y^{'})^{⊤}

; apply the Picard–Lindelöf theorem to the Lipschitz right-hand side.

has a unique solution existing on the entire interval $I$ .

Structure of the Solution Set

The homogeneous equation $L [y] = 0$ has a solution set that is a two-dimensional vector space. If $y_{1}, y_{2}$ are two linearly independent solutions (a fundamental set), then every homogeneous solution is

y_{h} = c_{1} y_{1} (x) + c_{2} y_{2} (x), c_{1}, c_{2} \in R .

⟳The map

y \mapsto (y (x_{0}), y^{'} (x_{0}))

from the kernel of

L

R^{2}

is a linear isomorphism by existence-uniqueness; hence the kernel has dimension 2.

The general solution of the inhomogeneous equation is

y = y_{h} + y_{p} = c_{1} y_{1} + c_{2} y_{2} + y_{p},

⟳If

L [y_{p}] = g

and

L [Y] = g

, then

L [Y - y_{p}] = 0

, so

Y - y_{p}

lies in the homogeneous space; rearrange.

where $y_{p}$ is any particular solution of $L [y] = g$ .

Wronskian and Linear Independence

The Wronskian of two solutions is

W (y_{1}, y_{2}) (x) = y_{1} y_{1}^{'} y_{2} y_{2}^{'} = y_{1} y_{2}^{'} - y_{2} y_{1}^{'} .

⟳The determinant tests whether the constants

c_{1}, c_{2}

matching given

(y, y^{'})

data are uniquely solvable; nonzero determinant ⇒ invertible coefficient matrix.

Two solutions of the same equation are linearly independent on $I$ if and only if $W \neq = 0$ at some (equivalently every) point of $I$ . The Wronskian obeys Abel's identity:

W (x) = W (x_{0}) exp (- \int_{x_{0}}^{x} p (t) d t) .

⟳Differentiate

W

, use the ODE to replace

y_{i}^{''}

- p y_{i}^{'} - q y_{i}

; cancellation yields

W^{'} = - p W

, a separable first-order ODE.

A profound consequence: the Wronskian is either identically zero or never zero — it cannot vanish at isolated points for solutions of the same linear equation.

Why It Holds / Derivation

Constant-Coefficient Homogeneous Solution via the Characteristic Equation

Consider $y^{''} + p y^{'} + q y = 0$ with constant $p, q$ .

Step 1. Seek exponential solutions

y = e^{r x}

, motivated by the fact that

e^{r x}

reproduces itself under differentiation. Then

y^{'} = r e^{r x}

y^{''} = r^{2} e^{r x}

. Step 2. Substitute:

(r^{2} + p r + q) e^{r x} = 0.

⟳Factor out the never-zero

e^{r x}

; the bracket must vanish.

Step 3. Since

e^{r x} \neq = 0

, we require the characteristic (auxiliary) equation

r^{2} + p r + q = 0, r_{1, 2} = \frac{- p \pm p ^{2} - 4 q}{2} .

⟳Quadratic formula; the discriminant

Δ = p^{2} - 4 q

governs the case structure.

Step 4 — Case analysis.

Distinct real roots ( $Δ > 0$ ): $y = c_{1} e^{r_{1} x} + c_{2} e^{r_{2} x}$ . Independence: $W = (r_{2} - r_{1}) e^{(r_{1} + r_{2}) x} \neq = 0$ .

Repeated root ( $Δ = 0$ , $r_{1} = r_{2} = r$ ): one solution $e^{r x}$ ; a second independent solution is $x e^{r x}$ , found by reduction of order:

y_{2} = x e^{r x}, L [x e^{r x}] = (2 r + p) e^{r x} = 0 since r = - p /2.

⟳Try

y_{2} = v (x) e^{r x}

; substitution forces

v^{''} = 0

, so

v = x

(taking the nontrivial part).

Complex conjugate roots ( $Δ < 0$ , $r = α \pm i β$ ): using Euler's formula $e^{i β x} = cos β x + i sin β x$ , the real fundamental set is

y = e^{α x} (c_{1} cos β x + c_{2} sin β x) .

⟳Take real and imaginary parts of the complex solution

e^{(α + i β) x}

; both solve the real ODE because

L

has real coefficients.

This exhausts all cases and produces a two-dimensional solution space, confirming the Structure Theorem. $■$

Variation of Parameters — Finding a Particular Solution

Given a fundamental set ${y_{1}, y_{2}}$ of $L [y] = 0$ , seek $y_{p} = u_{1} (x) y_{1} + u_{2} (x) y_{2}$ .

Step 1. Impose the convenient constraint

u_{1}^{'} y_{1} + u_{2}^{'} y_{2} = 0

to keep the algebra clean. Then

y_{p}^{'} = u_{1} y_{1}^{'} + u_{2} y_{2}^{'}

. Step 2. Differentiate again and substitute into

L [y_{p}] = g

. The homogeneous terms cancel, leaving

u_{1}^{'} y_{1}^{'} + u_{2}^{'} y_{2}^{'} = g .

⟳

L [y_{1}] = L [y_{2}] = 0

eliminates all undifferentiated

u_{i}

terms; only the derivatives of

u_{i}

survive.

Step 3. Solve the linear system (constraint + [13]) by Cramer's rule using the Wronskian:

u_{1}^{'} = - \frac{y _{2} g}{W}, u_{2}^{'} = \frac{y _{1} g}{W} .

⟳The

2 \times 2

coefficient matrix is exactly the Wronskian matrix; invert it.

Step 4. Integrate to obtain the particular solution:

y_{p} (x) = - y_{1} \int \frac{y _{2} g}{W} d x + y_{2} \int \frac{y _{1} g}{W} d x .

⟳Antidifferentiate [14]; constants of integration only add homogeneous solutions, already absorbed into

y_{h}

This always works whenever ${y_{1}, y_{2}}$ is known, regardless of the form of $g$ . $■$

Leonhard Euler | 1707–1783

Euler was the architect of the exponential-substitution method for linear ODEs with constant coefficients. In his Institutiones calculi integralis (1768–1770) he systematized the use of $y = e^{r x}$ and explicitly handled repeated and complex roots, employing the formula $e^{i θ} = cos θ + i sin θ$ that bears his name to convert complex exponentials into real trigonometric solutions. He also introduced the technique now called the integrating factor and studied the equidimensional (Euler–Cauchy) equation $x^{2} y^{''} + a x y^{'} + b y = 0$ . His correspondence with d'Alembert and Daniel Bernoulli on the vibrating string drove much of the theory forward, and his notation for differentials and functions remains standard.

Building Intuition

Think of a linear second-order ODE as describing a mass on a spring with possible friction and an external push. The term $y^{''}$ is acceleration (inertia), $p y^{'}$ is velocity-proportional damping, $q y$ is a restoring force, and $g$ is the external driving force. The characteristic roots are the system's "DNA": their real parts tell you whether disturbances grow or decay; their imaginary parts tell you whether the system oscillates.

Two negative real roots → sluggish, overdamped return to rest (think a door closer with stiff hydraulics).
A repeated negative root → critically damped, the fastest non-oscillating return.
Complex roots with negative real part → underdamped ringing that fades.
Pure imaginary roots → undamped perpetual oscillation.

The two-dimensional solution space corresponds to the two pieces of information a Newtonian system needs to predict the future: initial position and initial velocity. Once you fix those, the trajectory is locked in forever — that is uniqueness made physical.

Damping regimes of $y'' + 2\zeta\omega_0 y' + \omega_0^2 y = 0$ for fixed initial displacement

Worked Examples

Underdamped Mechanical Oscillator

Given: A mass–spring–damper with

m = 2 kg

, damping

c = 4 N\cdots/m

, stiffness

k = 50 N/m

, released from rest at displacement

x (0) = 0.1 m

x^{'} (0) = 0

. Equation:

m x^{''} + c x^{'} + k x = 0 \Rightarrow 2 x^{''} + 4 x^{'} + 50 x = 0 \Rightarrow x^{''} + 2 x^{'} + 25 x = 0

. Step 1 — Characteristic equation:

r^{2} + 2 r + 25 = 0

. Discriminant

Δ = 4 - 100 = - 96 < 0

. Step 2 — Roots:

r = \frac{- 2 \pm - 96}{2} = - 1 \pm 24 i = - 1 \pm 26 i

. So

α = - 1 s^{- 1}

β = 26 \approx 4.899 rad/s

. Step 3 — General solution:

x (t) = e^{- t} (c_{1} cos β t + c_{2} sin β t)

. Step 4 — Apply initial conditions:

x (0) = c_{1} = 0.1

x^{'} (t) = e^{- t} [(- c_{1} + β c_{2}) cos β t + (- c_{2} - β c_{1}) sin β t]

; setting

x^{'} (0) = 0

- c_{1} + β c_{2} = 0 \Rightarrow c_{2} = c_{1} / β = 0.1/4.899 = 0.0204

. Answer:

x (t) = e^{- t} (0.100 cos (4.899 t) + 0.0204 sin (4.899 t)) m

The amplitude decays as $e^{- t}$ (time constant $1 s$ ) while ringing at $\approx 0.78 Hz$ .

Forced RLC Circuit — Variation of Parameters

Given: Solve

y^{''} + 4 y = 8 sin (t)

with

y (0) = 0

y^{'} (0) = 0

(a circuit driven off-resonance; natural frequency

2

, drive frequency

1

). Step 1 — Homogeneous solution:

r^{2} + 4 = 0 \Rightarrow r = \pm 2 i

, so

y_{h} = c_{1} cos 2 t + c_{2} sin 2 t

. Here

y_{1} = cos 2 t, y_{2} = sin 2 t

W = y_{1} y_{2}^{'} - y_{2} y_{1}^{'} = 2 cos^{2} 2 t + 2 sin^{2} 2 t = 2

. Step 2 — Particular solution by method of undetermined coefficients (faster here): try

y_{p} = A sin t + B cos t

. Then

y_{p}^{''} = - A sin t - B cos t

, so

y_{p}^{''} + 4 y_{p} = 3 A sin t + 3 B cos t = 8 sin t

. Match:

3 A = 8 \Rightarrow A = 8/3

B = 0

. Step 3 — General solution:

y = c_{1} cos 2 t + c_{2} sin 2 t + \frac{8}{3} sin t

. Step 4 — Initial conditions:

y (0) = c_{1} = 0

y^{'} (t) = - 2 c_{1} sin 2 t + 2 c_{2} cos 2 t + \frac{8}{3} cos t

;

y^{'} (0) = 2 c_{2} + \frac{8}{3} = 0 \Rightarrow c_{2} = - \frac{4}{3}

. Answer:

y (t) = - \frac{4}{3} sin 2 t + \frac{8}{3} sin t

The motion is a superposition of the natural mode (frequency 2) and the steady driven response (frequency 1) — a clean beating pattern with no growth because the drive is off resonance.

Common Misconceptions

The sum of two solutions of an inhomogeneous equation is generally NOT a solution. Superposition holds only for the homogeneous equation. If

L [y_{1}] = g

and

L [y_{2}] = g

, then

L [y_{1} + y_{2}] = 2 g \neq = g

. Superposition of forcing terms does work in this sense: if

L [y_{a}] = g_{a}

and

L [y_{b}] = g_{b}

, then

L [y_{a} + y_{b}] = g_{a} + g_{b}

. Keep these two statements distinct. "Linear" does not mean "solvable in closed form." Linearity guarantees structural theorems (two-dimensional solution space, superposition), but variable-coefficient equations like Bessel's or Legendre's require special functions or series — there is no elementary formula in general. Two solutions that vanish together but are still independent? No — and here is why the Wronskian matters. Consider the functions

f (x) = x^{2}

and

g (x) = x ∣ x ∣

(- 1, 1)

. Their Wronskian is identically zero, yet they are linearly independent as abstract functions. This does not contradict the theory because they are not both solutions of the same linear second-order ODE with continuous coefficients on that interval — no such equation exists for which

x ∣ x ∣

is a

C^{2}

solution at the origin. The Wronskian's "all-or-nothing" property (Abel's identity) holds only for genuine solutions of one linear equation; it is precisely what distinguishes solution functions from arbitrary functions.

Connections

Linear second-order ODEs sit at a crossroads of analysis and algebra. They are the simplest nontrivial case of Linear Differential Equations of arbitrary order, and reduce to a planar Systems of First-Order ODEs via $u = (y, y^{'})^{⊤}$ , linking them to matrix exponentials and eigenvalue theory. The constant-coefficient case is the natural home of the Laplace Transform, which algebraizes the equation into a rational function of $s$ . Variable-coefficient cases with singular points open into the Frobenius Method and the theory of Special Functions (Bessel, Legendre, Hermite). The eigenvalue problem $L [y] = λ y$ with boundary conditions is the gateway to Sturm–Liouville Theory and the spectral decomposition underlying Fourier Series.

Real-World Applications

Mechanical vibrations / suspension design: $m \overset{x}{¨} + c \overset{x}{˙} + k x = F (t)$ models car suspensions, building seismic dampers, and machine isolation mounts; engineers tune $c$ for critical damping to suppress oscillation.
RLC electrical circuits: $L \overset{q}{¨} + R \overset{q}{˙} + q / C = V (t)$ governs charge in a series circuit; resonant tuning of radios exploits the underdamped peak.
Quantum mechanics (time-independent Schrödinger equation): $- \frac{ℏ ^{2}}{2 m} ψ^{''} + V (x) ψ = E ψ$ is a linear second-order ODE whose eigenvalues are quantized energy levels.
Beam deflection (Euler–Bernoulli, reduced): simply-supported beams under load yield second-order boundary value problems for transverse displacement.
Control systems: the transfer function of a second-order plant determines step-response overshoot and settling time; PID tuning is literally choosing $p$ and $q$ .
Population & pendulum dynamics (linearized): the small-angle pendulum $\ddot{θ} + (g / L) θ = 0$ and linearized predator–prey models reduce to constant-coefficient form.
Heat & wave equations via separation of variables: PDEs reduce to ODEs like $X^{''} + λ X = 0$ , whose eigenvalues set the allowed spatial modes.

Historical Development

1690s — The Bernoulli family and Leibniz solve early special second-order equations arising from the catenary and elastica problems.
1739 — Euler develops the exponential-substitution method for linear constant-coefficient equations and treats repeated and complex roots systematically.
1743 — Euler publishes the characteristic-equation technique; the term "characteristic" enters the lexicon. d'Alembert's work on the vibrating string (1747) supplies key physical motivation.
1760s–1770s — Lagrange invents variation of parameters and adjoint methods, generalizing particular-solution construction to arbitrary forcing.
1827 — Niels Henrik Abel proves Abel's identity for the Wronskian, clarifying the all-or-nothing nature of linear independence among solutions.
1836–1837 — Charles Sturm and Joseph Liouville develop the spectral theory of $(p y^{'})^{'} + (q + λ w) y = 0$ , founding Sturm–Liouville theory and oscillation theorems.
1838 — Friedrich Bessel's functions and the Bessel equation become central to variable-coefficient problems in physics and astronomy.
1868 — The Frobenius series method (refined later by Frobenius, 1873) provides solutions near regular singular points.
1890s — Picard and Lindelöf establish the modern existence-uniqueness theory via successive approximations, placing the entire subject on rigorous footing.

Subtleties & Limits

While the linear theory is essentially complete, several deep issues persist at its boundaries. Closed-form solvability is governed by differential Galois theory (Picard–Vessiot), which determines exactly when a linear ODE has solutions expressible in elementary or Liouvillian terms — but deciding this algorithmically for high-order or parametric families remains intricate. Singular points where $a (x) = 0$ (irregular singularities especially) generate solutions with essential singularities and Stokes phenomena, whose asymptotic behavior connects to resurgence theory and is still actively studied. The transition from linear to weakly nonlinear systems (e.g., the Duffing or van der Pol equations) shows that linear superposition is a fragile idealization; perturbation series can diverge, and chaos appears. Finally, the spectral theory of singular Sturm–Liouville operators on unbounded domains — central to quantum scattering — raises subtle questions about continuous spectra, self-adjoint extensions, and limit-point/limit-circle classification that blend ODE theory with functional analysis.

Linear Second-Order ODEs

The Core Insight

Formal Statement

Linear Second-Order ODE

Differential Operator and Linearity

Existence and Uniqueness (Cauchy–Picard)

Structure of the Solution Set

Wronskian and Linear Independence

Why It Holds / Derivation

Constant-Coefficient Homogeneous Solution via the Characteristic Equation

Variation of Parameters — Finding a Particular Solution

Leonhard Euler | 1707–1783

Building Intuition

Damping regimes of $y'' + 2\zeta\omega_0 y' + \omega_0^2 y = 0$ for fixed initial displacement

Worked Examples

Underdamped Mechanical Oscillator

Forced RLC Circuit — Variation of Parameters

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also

Before this

Where this leads