What Are We Solving?
A second-order linear ODE with constant coefficients looks like:
Let's break down every single symbol so nothing is mysterious:
| Symbol | Meaning |
|---|---|
| y | The unknown function we want to find — \(y = y(t)\), a function of time. |
| y′ | First derivative \(\frac{dy}{dt}\) — the rate of change of \(y\). Think: velocity. |
| y″ | Second derivative \(\frac{d^2y}{dt^2}\) — how fast the rate of change is changing. Think: acceleration. |
| a, b, c | Real constants. They "weight" each term. They don't change with time. |
| = 0 | Homogeneous — no external forcing. The equation is self-contained. |
A second-order ODE needs two initial conditions to pin down a unique solution — typically \(y(0)\) (initial position) and \(y'(0)\) (initial velocity). This is why the general solution will always have exactly two free constants.
Why Guess \(y = e^{rt}\)?
Before any formal machinery, ask a deeper question:
For the equation \(a y'' + b y' + c y = 0\) to be satisfied everywhere, all three terms must live in the same family — so they can cancel each other. You need a function whose derivatives are just scaled copies of itself.
The only function with that property is the exponential:
Each derivative gives back \(e^{rt}\), just scaled by a power of \(r\). This isn't a magical guess — it's the only structural choice that makes the cancellation possible. The ansatz (educated guess) \(y = e^{rt}\) is logically forced.
The Characteristic Equation
Dividing by \(e^{rt}\) converts our ODE into a plain quadratic in \(r\):
This is the characteristic equation. The ODE — a calculus problem — has become an algebra problem. The roots \(r\) tell you the growth rates or oscillation frequencies baked into the system.
The quadratic formula gives the roots:
The quantity \(\Delta = b^2 - 4ac\) is the discriminant. It's the fork in the road:
Two real distinct roots
Pure exponential behavior. No oscillation. The solution grows or decays monotonically (or a mix of both).
One repeated real root
The boundary case. The second independent solution picks up a factor of \(t\).
Complex conjugate roots
Oscillatory behavior. The imaginary part drives the wiggle; the real part controls decay or growth.
Real Distinct Roots (\(\Delta > 0\))
When \(\Delta > 0\), two distinct real roots \(r_1\) and \(r_2\) pop out of the quadratic formula. Each one gives a valid solution:
Why can you add them? — Superposition
The ODE is linear and homogeneous. That means if \(y_1\) satisfies the equation and \(y_2\) satisfies the equation, then so does any combination \(c_1 y_1 + c_2 y_2\). Let's verify this directly. Plug \(c_1 y_1 + c_2 y_2\) into the ODE:
Are they actually independent? — The Wronskian
Two solutions are linearly independent if one is not just a scalar multiple of the other. We check via the Wronskian determinant:
Plugging in:
Since \(r_1 \neq r_2\), we have \(r_2 - r_1 \neq 0\), and \(e^{(r_1+r_2)t}\) is never zero. So \(W \neq 0\) everywhere — the solutions are linearly independent. ✓
Two free constants \(c_1, c_2\) are determined by the initial conditions \(y(0)\) and \(y'(0)\).
What does the physics look like?
Each exponential either grows or decays depending on the sign of its root. If both roots are negative: both terms decay and the system settles. If one root is positive: the positive term eventually dominates and the solution blows up. This is an overdamped (or unstable) system.
Solve \(y'' - 5y' + 6y = 0\).
Euler's Formula: Circles from Exponentials
When \(\Delta < 0\), the quadratic formula gives complex roots. To handle them, we need the most important identity in mathematics:
What does this mean? Let's go step by step. The exponential \(e^{i\theta}\) takes an imaginary number as input. Normally, \(e^x\) for real \(x\) just stretches or compresses along the real line. But when the input is imaginary, the exponential rotates in the complex plane.
At angle \(\theta\):
- The real part \(\cos\theta\) oscillates between −1 and +1 horizontally.
- The imaginary part \(\sin\theta\) oscillates between −1 and +1 vertically.
- Together, they trace a unit circle in the complex plane.
Two immediately useful consequences:
These let us run the conversion in both directions. We'll need the forward direction next: turning a complex exponential into sines and cosines.
Complex Conjugate Roots (\(\Delta < 0\))
When \(\Delta = b^2 - 4ac < 0\), the square root produces an imaginary number. Define:
So the two roots are \(r = \alpha \pm \beta i\). Let's name what these mean physically before we do anything else:
| Symbol | Physical role |
|---|---|
| \(\alpha\) | Decay/growth rate. If \(\alpha < 0\): solution decays to zero (damped). If \(\alpha > 0\): solution grows unboundedly. If \(\alpha = 0\): amplitude stays constant forever. |
| \(\beta\) | Oscillation frequency. How many radians per unit time the solution wiggles. Higher \(\beta\) = faster wiggling. |
Step 1 — Write the complex solutions
Each root gives a solution \(e^{(\alpha \pm \beta i)t}\). Using Euler's formula:
Step 2 — Build real solutions by combining
Both complex exponentials are valid solutions, but they're complex-valued. We want real solutions. Use superposition: add them and subtract them:
These two new functions are:
- Real-valued — no imaginary parts remain.
- Both solutions to the ODE — linear combinations of solutions are solutions (superposition).
- Linearly independent — their Wronskian is \(\beta e^{2\alpha t} \neq 0\) (since \(\beta > 0\)).
So we have our two independent real solutions.
where \(\alpha = \operatorname{Re}(r)\) and \(\beta = \operatorname{Im}(r)\). The exponential \(e^{\alpha t}\) is the envelope; the trig functions are the oscillation.
Geometric picture: a spiral
Think of this as two parts working together. The \(e^{\alpha t}\) is an envelope — an exponential shell that either squeezes the solution toward zero (if \(\alpha < 0\)) or stretches it outward (if \(\alpha > 0\)). Inside that envelope, \(c_1\cos\beta t + c_2\sin\beta t\) is pure oscillation, spinning at frequency \(\beta\). Together: a damped spiral in phase space.
Underdamped
Wiggles while decaying to zero. Classic RLC circuit, shock absorber. Envelope shrinks.
Undamped
Pure oscillation forever. No energy loss. Ideal spring in a vacuum. Amplitude constant.
Negatively damped
Oscillation grows unboundedly. Unstable system. Envelope expands.
Solve \(y'' + 2y' + 5y = 0\).
The Full Map
Given \(ay'' + by' + cy = 0\), write \(ar^2 + br + c = 0\) and compute \(\Delta = b^2 - 4ac\):
| Discriminant \(\Delta\) | Roots | General Solution | Behavior |
|---|---|---|---|
| \(\Delta > 0\) | \(r_1, r_2 \in \mathbb{R}\), distinct | \(c_1 e^{r_1 t} + c_2 e^{r_2 t}\) | Pure exponential |
| \(\Delta = 0\) | \(r = -b/2a\) (repeated) | \((c_1 + c_2 t)\,e^{rt}\) | Exponential × polynomial |
| \(\Delta < 0\) | \(r = \alpha \pm \beta i\) | \(e^{\alpha t}(c_1\cos\beta t + c_2\sin\beta t)\) | Oscillation with envelope |
Interactive Solver
Dial in the coefficients and initial conditions. The solver computes the roots, classifies the case, shows the closed-form solution with your exact constants, and plots \(y(t)\) — all in real time.