Complex and Repeated Roots

What function, added to its own second derivative, gives zero? Why does that produce a wave? This question forces us past the comfortable case of distinct real roots. It pushes us into the two situations that make second-order linear equations genuinely interesting. These are complex conjugate roots, which produce oscillation, and repeated roots, which produce resonance-like growth against a decaying background. Both cases arise directly from applying the same characteristic-equation ansatz used for Homogeneous Equations with Constant Coefficients. But each breaks the naive guess y=erty = e^{rt} in a different way. One breaks it because rr leaves the real numbers. The other breaks it because rr collides with itself. These two cases are not exotic edge cases to be tolerated. They are the mathematical description of every damped and undamped oscillator. They describe every critically damped shock absorber, and every RLC circuit at resonance. A mass on a spring with no friction oscillates forever. This happens because its characteristic roots are pure imaginary. A door closer returns the door to rest as fast as possible without overshooting. Such a closer is tuned to a repeated real root. Understanding these cases completely means more than just pattern-matching a formula. It lets you read the physical behavior of a system directly off its characteristic polynomial. The big picture is this: a second-order linear homogeneous equation with constant coefficients always has a two-dimensional solution space. The characteristic polynomial's roots always parametrize a basis for it. This holds no matter what kind of roots they are. Complex roots come in conjugate pairs by the Fundamental Theorem of Algebra applied to a real polynomial. Euler's formula converts the resulting complex exponentials into real sines and cosines. Repeated roots require an extra idea, reduction of order. This idea manufactures the missing second solution that a naive exponential ansatz cannot supply. This article proves both extensions from first principles. It shows why eαtcos⁡βte^{\alpha t}\cos\beta t and eαtsin⁡βte^{\alpha t}\sin\beta t solve the equation when the roots are α±iβ\alpha \pm i\beta. It shows why tertt e^{rt} solves the equation when rr is a repeated root. It builds a complete general theory. It checks that theory with two asymmetric worked examples. It shows exactly what fails if the discriminant hypothesis is dropped.

Key takeaways
  • A second-order linear ODE with constant coefficients, ay′′+by′+cy=0ay'' + by' + cy = 0, has characteristic equation ar2+br+c=0ar^2+br+c=0.
  • If the discriminant b2−4ac<0b^2-4ac<0, roots are complex conjugates α±iβ\alpha\pm i\beta; solutions are eαtcos⁡βte^{\alpha t}\cos\beta t and eαtsin⁡βte^{\alpha t}\sin\beta t.
  • If the discriminant is zero, there is one repeated real root r=−b/2ar=-b/2a; solutions are erte^{rt} and tertte^{rt}.
  • Reduction of order is the general machine that manufactures a second solution from one known solution.
  • Complex roots encode oscillation frequency (β\beta) and exponential envelope rate (α\alpha) simultaneously.
  • The Wronskian test confirms linear independence in both cases, guaranteeing the general solution is complete.
  • Key idea: the characteristic polynomial's roots, however they are arranged, always parametrize a full basis of solutions — complex roots via Euler's formula, repeated roots via a tt-multiplied twin.

Where this lives today

Complex characteristic roots describe every underdamped mechanical oscillator: a mass–spring–damper system mx′′+cx′+kx=0mx''+cx'+kx=0 has complex roots exactly when damping is light, giving the familiar decaying-sine motion seen in shock absorbers and seismographs. In electrical engineering, the same mathematics governs RLC circuits, where α\alpha is the decay rate of transient current and β\beta is the ringing frequency; resonance and filter design in signal processing depend on locating these roots in the complex plane. Repeated real roots correspond to critical damping, the precise boundary between oscillatory and non-oscillatory return to equilibrium, used deliberately in the design of door closers, analog meter needles, and control systems tuned to settle without overshoot as fast as possible. In control theory, the location of characteristic (here, closed-loop pole) roots in the complex plane — real part negative for stability, imaginary part for oscillation frequency — is the entire content of root-locus and pole-placement design, a direct descendant of this classification. In quantum mechanics, the time-independent Schrödinger equation for bound states produces eigenvalue problems whose characteristic-type analysis (though for boundary value problems, not IVPs) still hinges on the same real/complex root dichotomy. The reduction-of-order technique used to build tertte^{rt} generalizes directly to constructing missing solutions for Higher-Order Linear Equations with roots of higher multiplicity, and both root types reappear as the "homogeneous part" one must add to any particular solution found by the Method of Undetermined Coefficients.

Ask the tutor: How does this classification of characteristic roots extend to third- and fourth-order constant-coefficient equations with mixed real, repeated, and complex roots?

DefinitionSecond-order linear homogeneous constant-coefficient equation, its characteristic equation, and its solution basis

Let a,b,ca,b,c be real constants with a≠0a \neq 0. Consider the differential equation

ay′′+by′+cy=0,t∈R.a y'' + b y' + c y = 0, \qquad t \in \mathbb{R}.

The characteristic equation is the algebraic equation obtained by substituting the trial solution y=erty=e^{rt}:

ar2+br+c=0.ar^2 + br + c = 0.

Its roots are

r1,2=−b±b2−4ac2a.r_{1,2} = \frac{-b \pm \sqrt{b^2-4ac}}{2a}.

Three cases are classified by the discriminant Δ=b2−4ac\Delta = b^2-4ac:

  • Δ>0\Delta > 0: two distinct real roots r1≠r2r_1 \neq r_2 (covered by the constant-coefficient basics).
  • Δ<0\Delta < 0: a complex conjugate pair r1,2=α±iβr_{1,2} = \alpha \pm i\beta with α=−b2a\alpha = -\frac{b}{2a}, β=4ac−b22a>0\beta = \frac{\sqrt{4ac-b^2}}{2a} > 0.
  • Δ=0\Delta = 0: a single repeated real root r=−b2ar = -\frac{b}{2a} of multiplicity two.

A fundamental set of solutions is a pair y1,y2y_1,y_2 solving the equation whose Wronskian W(y1,y2)(t)=y1y2′−y1′y2W(y_1,y_2)(t) = y_1 y_2' - y_1' y_2 is nonzero for all (equivalently, for one) tt. Given a fundamental set, the general solution is

y(t)=C1y1(t)+C2y2(t),C1,C2∈R arbitrary,y(t) = C_1 y_1(t) + C_2 y_2(t), \qquad C_1, C_2 \in \mathbb{R} \text{ arbitrary},

and every solution of the equation arises this way for a unique choice of C1,C2C_1, C_2.

Intuition

Start from what already works: when the characteristic equation has two distinct real roots r1,r2r_1, r_2, the functions er1te^{r_1 t} and er2te^{r_2 t} are visibly independent — one grows or decays at one rate, the other at a different rate, and no constant multiple of one can equal the other. The whole difficulty of complex and repeated roots is that this visible independence disappears. When roots are complex, er1te^{r_1t} and er2te^{r_2t} are complex-valued functions of a real variable, and we need genuinely real solutions to describe a real physical system. When roots collide, er1te^{r_1t} and er2te^{r_2t} become the same function, and we have lost half our solution space. The fix for complex roots is Euler's formula, already familiar: eiθ=cos⁡θ+isin⁡θe^{i\theta} = \cos\theta + i\sin\theta. If r=α+iβr=\alpha+i\beta solves the characteristic equation, then ert=eαt(cos⁡βt+isin⁡βt)e^{rt} = e^{\alpha t}(\cos\beta t + i \sin \beta t) solves the ODE as a complex-valued function. Because the ODE has real coefficients, the real and imaginary parts of any complex solution are themselves real solutions — a fact that converts one complex exponential into two real, independent, oscillating solutions. The picture is a spiral: α\alpha controls how fast the spiral expands or contracts (the envelope eαte^{\alpha t}), and β\beta controls how fast it winds around (the oscillation frequency). The fix for repeated roots is different in spirit: it is not an algebraic trick but a systematic method, reduction of order, that manufactures a second solution out of thin air once one solution is known. The intuition is that near a repeated root, the two nearby-but-distinct roots rr and r+ϵr+\epsilon of a perturbed problem give solutions erte^{rt} and e(r+ϵ)te^{(r+\epsilon)t}; the combination e(r+ϵ)t−ertϵ\frac{e^{(r+\epsilon)t}-e^{rt}}{\epsilon} stays finite as ϵ→0\epsilon \to 0 and limits to tertte^{rt}, which is exactly the "extra" solution needed to fill the two-dimensional gap.

Figure
Complex roots as a spiral: envelope eαte^{\alpha t} modulating oscillation of frequency β\beta
t y y = e^{αt} cos βt envelope ±e^{αt} α < 0: decaying oscillation (spiral sink)
Figure
Repeated root: erte^{rt} and t ertt\,e^{rt} as two genuinely different curves for r=−1r=-1
t y y₁ = e^{-t} y₂ = t e^{-t} (peaks at t=1)
TheoremGeneral solution for complex and for repeated characteristic roots

Let a,b,c∈Ra,b,c \in \mathbb{R}, a≠0a \neq 0, and consider

ay′′+by′+cy=0.ay'' + by' + cy = 0.

(Complex case.) Suppose Δ=b2−4ac<0\Delta = b^2-4ac < 0, so the characteristic equation has roots r1,2=α±iβr_{1,2} = \alpha \pm i\beta with α=−b/2a∈R\alpha = -b/2a \in \mathbb{R} and β=4ac−b2/(2a)>0\beta = \sqrt{4ac-b^2}/(2a) > 0. Then

y1(t)=eαtcos⁡(βt),y2(t)=eαtsin⁡(βt)y_1(t) = e^{\alpha t}\cos(\beta t), \qquad y_2(t) = e^{\alpha t}\sin(\beta t)

are solutions of the ODE, they are linearly independent on R\mathbb{R}, and every solution has the unique form

y(t)=eαt(C1cos⁡βt+C2sin⁡βt),C1,C2∈R.y(t) = e^{\alpha t}\big(C_1 \cos\beta t + C_2 \sin \beta t\big), \qquad C_1,C_2 \in \mathbb{R}.

(Repeated case.) Suppose Δ=b2−4ac=0\Delta = b^2-4ac = 0, so the characteristic equation has the single root r=−b/2ar = -b/2a of multiplicity two. Then

y1(t)=ert,y2(t)=t erty_1(t) = e^{rt}, \qquad y_2(t) = t\,e^{rt}

are solutions of the ODE, they are linearly independent on R\mathbb{R}, and every solution has the unique form

y(t)=(C1+C2t) ert,C1,C2∈R.y(t) = (C_1 + C_2 t)\,e^{rt}, \qquad C_1, C_2 \in \mathbb{R}.
Proof

We prove both cases. Throughout, we use the standard existence–uniqueness theorem for linear ODEs (a linear equation with continuous, here constant, coefficients has a unique solution for each initial condition y(t0)=y0, y′(t0)=y0′y(t_0)=y_0,\ y'(t_0)=y_0', defined for all tt), and the fact that if y1,y2y_1,y_2 solve a linear homogeneous equation with nonvanishing Wronskian at one point, then their Wronskian never vanishes and {y1,y2}\{y_1,y_2\} spans the full two-dimensional solution space. This last fact is proved from Abel's identity, W′=−baWW' = -\frac{b}{a}W, whose solution W(t)=W(t0)e−ba(t−t0)W(t) = W(t_0)e^{-\frac{b}{a}(t-t_0)} is never zero if it is nonzero once; we take this as an already-established structural fact about second-order linear equations, not specific to root type.

Step 1 (Complex case): verifying y1,y2y_1, y_2 solve the equation. Work first with the complex exponential z(t)=ertz(t) = e^{rt} where r=α+iβr=\alpha+i\beta, treated as a C\mathbb{C}-valued function of the real variable tt, differentiated term by term: z′(t)=rertz'(t) = re^{rt}, z′′(t)=r2ertz''(t) = r^2 e^{rt}. Substituting into the ODE (whose coefficients a,b,ca,b,c are real, so this substitution is valid formally regardless of rr's being complex):

az′′+bz′+cz=(ar2+br+c)ert=0⋅ert=0,az'' + bz' + cz = (ar^2+br+c)e^{rt} = 0 \cdot e^{rt} = 0,

since rr satisfies the characteristic equation by hypothesis. So z(t)=ertz(t)=e^{rt} is a (complex-valued) solution. By Euler's formula, ert=e(α+iβ)t=eαteiβt=eαt(cos⁡βt+isin⁡βt)e^{rt} = e^{(\alpha+i\beta)t} = e^{\alpha t}e^{i\beta t} = e^{\alpha t}(\cos\beta t + i \sin \beta t). Write z(t)=u(t)+iv(t)z(t) = u(t) + iv(t) where u(t)=eαtcos⁡βtu(t) = e^{\alpha t}\cos\beta t and v(t)=eαtsin⁡βtv(t) = e^{\alpha t}\sin\beta t are real-valued. Since a,b,ca,b,c are real and zz solves the equation:

0=az′′+bz′+cz=(au′′+bu′+cu)+i(av′′+bv′+cv).0 = az'' + bz' + cz = (au''+bu'+cu) + i(av''+bv'+cv).

A complex number is zero if and only if its real and imaginary parts are both zero. Hence

au′′+bu′+cu=0andav′′+bv′+cv=0,au'' + bu' + cu = 0 \quad \text{and} \quad av'' + bv' + cv = 0,

i.e. y1=uy_1 = u and y2=vy_2 = v both solve the ODE. This is the general principle: for a real linear operator, the real and imaginary parts of a complex solution are real solutions.

Step 2 (Complex case): linear independence via the Wronskian. Compute at t=0t=0:

y1(0)=1,y1′(0)=α,y2(0)=0,y2′(0)=β.y_1(0) = 1, \quad y_1'(0) = \alpha, \quad y_2(0) = 0, \quad y_2'(0) = \beta.

(These follow directly by differentiating u,vu,v using the product rule and evaluating at t=0t=0, where cos⁡0=1,sin⁡0=0\cos 0=1,\sin 0=0.) So

W(y1,y2)(0)=y1(0)y2′(0)−y1′(0)y2(0)=(1)(β)−(α)(0)=β.W(y_1,y_2)(0) = y_1(0)y_2'(0) - y_1'(0)y_2(0) = (1)(\beta) - (\alpha)(0) = \beta.

Since Δ<0\Delta<0 forces β=4ac−b2/(2a)≠0\beta = \sqrt{4ac-b^2}/(2a) \neq 0 (indeed β>0\beta>0 by our sign convention), W(y1,y2)(0)≠0W(y_1,y_2)(0)\neq 0. By the Wronskian criterion, y1,y2y_1,y_2 are linearly independent and form a fundamental set, so the general solution is y=C1y1+C2y2y = C_1 y_1 + C_2 y_2 as claimed, and this representation is unique for each solution by the existence–uniqueness theorem (matching C1,C2C_1,C_2 to initial data y(0),y′(0)y(0),y'(0) via the nonzero Wronskian, i.e. invertible linear system).

Step 3 (Repeated case): verifying y1=erty_1=e^{rt} solves the equation. When Δ=0\Delta=0, the characteristic equation ar2+br+c=0ar^2+br+c=0 factors as a(r−r0)2=0a(r-r_0)^2=0 where r0=−b/2ar_0=-b/2a is the repeated root, so c=ar02c = ar_0^2 and b=−2ar0b=-2ar_0. Direct substitution of y1=er0ty_1=e^{r_0t} gives ay1′′+by1′+cy1=(ar02+br0+c)er0t=0ay_1''+by_1'+cy_1=(ar_0^2+br_0+c)e^{r_0t}=0, exactly as for any characteristic root.

Step 4 (Repeated case): constructing y2=ter0ty_2 = t e^{r_0t} by reduction of order. Since we only have one solution y1=er0ty_1=e^{r_0t}, we seek a second solution of the form y2(t)=v(t) y1(t)=v(t)er0ty_2(t) = v(t)\,y_1(t) = v(t)e^{r_0t} for an unknown function vv, and derive what vv must be — this is the method of reduction of order. Compute:

y2′=v′er0t+r0ver0t,y2′′=v′′er0t+2r0v′er0t+r02ver0t.y_2' = v'e^{r_0t} + r_0 v e^{r_0t}, \qquad y_2'' = v''e^{r_0t} + 2r_0v'e^{r_0t} + r_0^2ve^{r_0t}.

Substitute into ay2′′+by2′+cy2=0ay_2''+by_2'+cy_2=0 and divide through by er0te^{r_0t} (never zero):

a(v′′+2r0v′+r02v)+b(v′+r0v)+cv=0.a(v'' + 2r_0v' + r_0^2 v) + b(v'+r_0v) + cv = 0.

Group by derivative order of vv:

av′′+(2ar0+b)v′+(ar02+br0+c)v=0.a v'' + (2ar_0+b)v' + (ar_0^2+br_0+c)v = 0.

The coefficient of vv vanishes because r0r_0 is a characteristic root: ar02+br0+c=0ar_0^2+br_0+c=0. The coefficient of v′v' vanishes because r0=−b/2ar_0=-b/2a is exactly the repeated root, so 2ar0+b=2a(−b2a)+b=−b+b=02ar_0+b = 2a\left(-\frac{b}{2a}\right)+b = -b+b=0. (This double vanishing is precisely what makes r0r_0 a repeated, not merely a, root — it is a double root of the characteristic polynomial exactly when it is also a root of the derivative 2ar+b2ar+b.) We are left with

av′′=0⟹v′′=0(since a≠0).av'' = 0 \quad\Longrightarrow\quad v''=0 \quad (\text{since } a\neq 0).

Integrating twice, v(t)=C1+C2tv(t) = C_1 + C_2 t for arbitrary constants. Any nonconstant choice gives a genuinely new solution; the simplest is v(t)=tv(t)=t, giving

y2(t)=ter0t.y_2(t) = t e^{r_0 t}.

Substituting back confirms directly: y2′=er0t(1+r0t)y_2' = e^{r_0t}(1+r_0t), y2′′=er0t(2r0+r02t)y_2'' = e^{r_0t}(2r_0+r_0^2t), and

ay2′′+by2′+cy2=er0t[t(ar02+br0+c)+(2ar0+b)]=er0t[t⋅0+0]=0,ay_2''+by_2'+cy_2 = e^{r_0t}\Big[t(ar_0^2+br_0+c) + (2ar_0+b)\Big] = e^{r_0t}[t\cdot 0 + 0] = 0,

using both vanishing coefficients found above.

Step 5 (Repeated case): linear independence. Compute the Wronskian at t=0t=0: y1(0)=1y_1(0)=1, y1′(0)=r0y_1'(0)=r_0, y2(0)=0y_2(0)=0, y2′(0)=e0(1+r0⋅0)=1y_2'(0) = e^{0}(1+r_0\cdot 0) = 1. So

W(y1,y2)(0)=(1)(1)−(r0)(0)=1≠0.W(y_1,y_2)(0) = (1)(1) - (r_0)(0) = 1 \neq 0.

Hence y1,y2y_1,y_2 are linearly independent, form a fundamental set, and the general solution is y=(C1+C2t)er0ty=(C_1+C_2t)e^{r_0t}, uniquely determined by initial data via the same argument as Step 2. This completes the proof of both cases. ■\blacksquare

∎

Worked examples

Example 1 (complex roots, asymmetric coefficients). Solve the initial value problem

2y′′−4y′+10y=0,y(0)=3,y′(0)=1.2y'' - 4y' + 10y = 0, \qquad y(0) = 3, \quad y'(0) = 1.

Characteristic equation: 2r2−4r+10=02r^2-4r+10=0, i.e. r2−2r+5=0r^2-2r+5=0. Discriminant: Δ=(−2)2−4(1)(5)=4−20=−16<0\Delta = (-2)^2-4(1)(5) = 4-20=-16<0, confirming the complex case. Roots:

r=2±−162=1±2i.r = \frac{2 \pm \sqrt{-16}}{2} = 1 \pm 2i.

So α=1\alpha=1, β=2\beta=2. By the theorem, a fundamental set is y1=etcos⁡2ty_1=e^{t}\cos 2t, y2=etsin⁡2ty_2=e^t\sin 2t, and the general solution is

y(t)=et(C1cos⁡2t+C2sin⁡2t).y(t) = e^{t}\big(C_1\cos 2t + C_2 \sin 2t\big).

Applying initial conditions. At t=0t=0: y(0)=C1=3y(0)=C_1 = 3, so C1=3C_1=3. Differentiate: y′(t)=et(C1cos⁡2t+C2sin⁡2t)+et(−2C1sin⁡2t+2C2cos⁡2t)y'(t) = e^t(C_1\cos2t+C_2\sin2t) + e^t(-2C_1\sin2t+2C_2\cos2t). At t=0t=0:

y′(0)=C1+2C2=1.y'(0) = C_1 + 2C_2 = 1.

With C1=3C_1=3: 3+2C2=1⇒C2=−13+2C_2=1 \Rightarrow C_2 = -1. Solution:

y(t)=et(3cos⁡2t−sin⁡2t).y(t) = e^{t}\big(3\cos 2t - \sin 2t\big).

Check: y(0)=1⋅(3−0)=3y(0)=1\cdot(3-0)=3 ✓. Compute y′(t)y'(t) directly: y′=et(3cos⁡2t−sin⁡2t)+et(−6sin⁡2t−2cos⁡2t)=et(cos⁡2t−7sin⁡2t)y' = e^t(3\cos2t-\sin2t) + e^t(-6\sin2t-2\cos2t) = e^t(\cos2t - 7\sin2t); at t=0t=0, y′(0)=1y'(0)=1 ✓. One can also verify 2y′′−4y′+10y=02y''-4y'+10y=0 symbolically by substituting y=et(3cos⁡2t−sin⁡2t)y=e^t(3\cos2t-\sin2t) and using that both etcos⁡2t,etsin⁡2te^t\cos2t,e^t\sin2t individually solve the equation, so any real linear combination does too, by linearity of the ODE operator.

Example 2 (repeated roots, asymmetric coefficients). Solve

y′′+6y′+9y=0,y(−1)=2,y′(−1)=−5.y'' + 6y' + 9y = 0, \qquad y(-1) = 2, \quad y'(-1) = -5.

Note the initial condition is posed at t=−1t=-1, not t=0t=0, to test the general solution honestly. Characteristic equation: r2+6r+9=0r^2+6r+9=0, i.e. (r+3)2=0(r+3)^2=0. Discriminant Δ=36−36=0\Delta = 36-36=0, confirming the repeated case, with repeated root r=−3r=-3. General solution (by the theorem): y(t)=(C1+C2t)e−3ty(t) = (C_1+C_2t)e^{-3t}. Apply y(−1)=2y(-1)=2:

y(−1)=(C1−C2)e3=2⟹C1−C2=2e−3.y(-1) = (C_1 - C_2)e^{3} = 2 \quad\Longrightarrow\quad C_1 - C_2 = 2e^{-3}.

Differentiate: y′(t)=C2e−3t−3(C1+C2t)e−3t=e−3t[(C2−3C1)−3C2t]y'(t) = C_2e^{-3t} - 3(C_1+C_2t)e^{-3t} = e^{-3t}\big[(C_2-3C_1) - 3C_2t\big]. At t=−1t=-1:

y′(−1)=e3[(C2−3C1)+3C2]=e3(4C2−3C1)=−5.y'(-1) = e^{3}\big[(C_2-3C_1)+3C_2\big] = e^3(4C_2-3C_1) = -5.

So 4C2−3C1=−5e−34C_2 - 3C_1 = -5e^{-3}. Solve the linear system. From the first equation, C1=C2+2e−3C_1 = C_2 + 2e^{-3}. Substitute:

4C2−3(C2+2e−3)=−5e−3  ⟹  C2−6e−3=−5e−3  ⟹  C2=e−3.4C_2 - 3(C_2+2e^{-3}) = -5e^{-3} \;\Longrightarrow\; C_2 - 6e^{-3} = -5e^{-3} \;\Longrightarrow\; C_2 = e^{-3}.

Then C1=e−3+2e−3=3e−3C_1 = e^{-3}+2e^{-3} = 3e^{-3}. Solution:

y(t)=(3e−3+e−3t)e−3t=e−3(3+t)e−3t=(3+t)e−3(t+1).y(t) = \big(3e^{-3} + e^{-3}t\big)e^{-3t} = e^{-3}(3+t)e^{-3t} = (3+t)e^{-3(t+1)}.

Check at t=−1t=-1: y(−1)=(3−1)e0=2y(-1) = (3-1)e^{0} = 2 ✓. Compute y′(t)=e−3(t+1)−3(3+t)e−3(t+1)=e−3(t+1)[1−3(3+t)]=e−3(t+1)(−8−3t)y'(t) = e^{-3(t+1)} - 3(3+t)e^{-3(t+1)} = e^{-3(t+1)}\big[1-3(3+t)\big] = e^{-3(t+1)}(-8-3t). At t=−1t=-1: y′(−1)=e0(−8+3)=−5y'(-1) = e^0(-8+3) = -5 ✓. This confirms the general repeated-root solution correctly absorbs an off-origin initial point.

CounterexampleDropping the "repeated root exactly at r0=−b/2ar_0=-b/2a" hypothesis: why tertte^{rt} fails for a merely nearby, non-repeated root

The proof's Step 4 crucially used 2ar0+b=02ar_0+b=0, which holds only when r0r_0 is exactly the double root −b/2a-b/2a, not an approximation to it. To see the hypothesis is not cosmetic, take the equation

y′′−5y′+6y=0,y'' - 5y' + 6y = 0,

whose characteristic equation r2−5r+6=(r−2)(r−3)=0r^2-5r+6=(r-2)(r-3)=0 has two distinct real roots r=2,3r=2,3 — discriminant Δ=25−24=1≠0\Delta = 25-24=1\neq0, so this is emphatically not the repeated-root case. Suppose someone, misremembering the theorem, guesses that because 22 is a root, y2=te2ty_2 = te^{2t} should also be a solution (treating 22 as if it were a repeated root). Test it directly: with y2=te2ty_2=te^{2t},

y2′=e2t(1+2t),y2′′=e2t(4+4t).y_2' = e^{2t}(1+2t), \qquad y_2'' = e^{2t}(4+4t).

Substitute into the left side:

y2′′−5y2′+6y2=e2t[(4+4t)−5(1+2t)+6t]=e2t[(4−5)+(4t−10t+6t)]=e2t(−1+0)=−e2t≠0.y_2''-5y_2'+6y_2 = e^{2t}\big[(4+4t) - 5(1+2t) + 6t\big] = e^{2t}\big[(4-5) + (4t-10t+6t)\big] = e^{2t}(-1+0) = -e^{2t} \neq 0.

So te2tte^{2t} is not a solution of this equation. The reduction-of-order computation in the proof shows exactly why: substituting y=v(t)e2ty=v(t)e^{2t} into this equation would produce, in place of Step 4's clean av′′=0av''=0, an equation with a nonzero coefficient on v′v' (since 2ar0+b=2(2)+(−5)=−1≠02a r_0 + b = 2(2)+(-5) = -1 \neq 0 here), forcing v′v' to satisfy a nontrivial first-order linear ODE whose solution is v(t)=C1+C2etv(t) = C_1 + C_2 e^{t} rather than v(t)=C1+C2tv(t)=C_1+C_2t. Indeed the true fundamental set here is {e2t,e3t}\{e^{2t}, e^{3t}\} (distinct-root case), and te2tte^{2t} is genuinely outside the solution space. The lesson: the multiplier tt is licensed only by the double vanishing of both characteristic-polynomial and its derivative at r0r_0, a feature unique to genuine repeated roots.

Scope, edge cases & common traps

SituationCorrect treatmentCommon error
Δ<0\Delta<0Use eαtcos⁡βt, eαtsin⁡βte^{\alpha t}\cos\beta t,\ e^{\alpha t}\sin\beta tLeaving the answer as complex exponentials e(α±iβ)te^{(\alpha\pm i\beta)t} without converting to real form
Δ=0\Delta=0Use ert, terte^{rt},\ te^{rt}Writing C1ert+C2ertC_1e^{rt}+C_2e^{rt} (collapses to one constant — not a valid basis)
Δ>0\Delta>0, roots close togetherStill use er1t,er2te^{r_1t}, e^{r_2t}Prematurely invoking tertte^{rt} before checking Δ=0\Delta=0 exactly
Sign of β\betaβ=4ac−b2/(2a)\beta=\sqrt{4ac-b^2}/(2a); taking −β-\beta just swaps sign of C2C_2, same solution spaceForgetting β\beta must use 4ac−b24ac-b^2 (positive) inside the root, not b2−4acb^2-4ac
Higher-order repeated rootsA root of multiplicity kk contributes ert,tert,…,tk−1erte^{rt}, te^{rt}, \dots, t^{k-1}e^{rt}Assuming only t1t^1 ever appears, even for multiplicity ≥3\geq 3
Initial conditions not at t=0t=0Substitute the actual t0t_0 into general solution and its derivativeSilently assuming initial data is always given at t=0t=0

The mathematician behind it

Origins · 1 minLeonhard Euler and how this was discovered
LE
Leonhard Euler
1707–1783 · 18th-century analysis, mechanics, and the identity connecting exponentials to trigonometry

Euler's identity eiθ=cos⁡θ+isin⁡θe^{i\theta}=\cos\theta+i\sin\theta, introduced in his 1748 Introductio in analysin infinitorum, is the single algebraic fact that converts complex characteristic roots into real oscillatory solutions. Euler developed the calculus of complex exponentials alongside his systematic study of linear differential equations with constant coefficients, effectively inventing the characteristic-equation method used throughout this article. His work on solving linear recurrences and differential equations via auxiliary polynomial equations directly anticipates the modern treatment of repeated and complex roots.

History
- 1740: Euler and Daniel Bernoulli correspond on solving linear ODEs via exponential trial solutions, an early form of the characteristic equation.
- 1743: Euler publishes a systematic method for nnth-order linear equations with constant coefficients, including the repeated-root case.
- 1748: Euler's Introductio establishes eiθ=cos⁡θ+isin⁡θe^{i\theta}=\cos\theta+i\sin\theta, later called Euler's formula.
- 1799: Gauss proves the Fundamental Theorem of Algebra, guaranteeing every real polynomial's complex roots occur in conjugate pairs.
- 1830s: Cauchy formalizes existence–uniqueness theory for linear ODEs, placing the fundamental-set argument on rigorous footing.
Why do complex characteristic roots always come in conjugate pairs?

Because a,b,ca,b,c are real numbers, the characteristic polynomial ar2+br+car^2+br+c has real coefficients. If r=α+iβr=\alpha+i\beta (with β≠0\beta\neq0) is a root, take the complex conjugate of the equation ar2+br+c=0ar^2+br+c=0: since conjugation distributes over sums and products and fixes real numbers, arˉ2+brˉ+c=ar2+br+c‾=0ˉ=0a\bar r^2+b\bar r+c = \overline{ar^2+br+c}=\bar 0=0. So rˉ=α−iβ\bar r=\alpha-i\beta is also a root. This is a general fact about polynomials with real coefficients, not specific to second order, and it guarantees the two roots of a quadratic with negative discriminant are exactly a conjugate pair, never two unrelated complex numbers.

Can the constant α\alpha in eαtcos⁡βte^{\alpha t}\cos\beta t be zero, and what does that mean physically?

Yes: α=−b/2a=0\alpha=-b/2a=0 exactly when b=0b=0, i.e. there is no first-derivative (damping) term. Then the solutions reduce to pure cos⁡βt\cos\beta t and sin⁡βt\sin\beta t, describing undamped simple harmonic motion that neither grows nor decays — the idealized frictionless spring or lossless LC circuit. Whenever α<0\alpha<0 (which happens exactly when a,ba,b have the same sign, i.e. positive damping in the usual physical sign convention), the oscillation decays; α>0\alpha>0 would describe an unstable, growing oscillation, which arises in some driven or actively-forced systems but not in a passive damped oscillator.

Is reduction of order only useful for the repeated-root case?

No — reduction of order is a general technique for turning any one known solution y1y_1 of a second-order linear homogeneous equation (even with non-constant coefficients) into a second, independent solution y2=v(t)y1(t)y_2=v(t)y_1(t), by substituting and solving the resulting first-order equation for v′v'. The repeated-root case is the cleanest illustration because the resulting equation for vv simplifies all the way to v′′=0v''=0. In more general variable-coefficient equations, the equation for v′v' is a nontrivial but always first-order (hence directly integrable or separable) linear ODE, which is why the method is called "reduction of order": it reduces a second-order problem to a first-order one.

What happens for a third-order or higher equation with a repeated root, say multiplicity 3?

The pattern generalizes: a real root rr of multiplicity kk in the characteristic equation of an nnth-order constant-coefficient linear ODE contributes exactly kk independent solutions,

ert, tert, t2ert, …, tk−1ert.e^{rt},\ te^{rt},\ t^2e^{rt},\ \dots,\ t^{k-1}e^{rt}.

This is proved by an extension of the same reduction-of-order argument (or, more efficiently, by noting the characteristic polynomial and its first k−1k-1 derivatives all vanish at rr, matching the vanishing coefficients seen in the second-order proof). This full statement is developed further under Higher-Order Linear Equations.

Why is the Wronskian test, rather than "eyeballing" independence, necessary to trust the general solution?

Because independence can be subtle to see by inspection, especially for oscillatory or nearly-identical-looking functions — for instance eαtcos⁡βte^{\alpha t}\cos\beta t and eαtsin⁡βte^{\alpha t}\sin\beta t are not obviously "different rates" the way er1t,er2te^{r_1t}, e^{r_2t} are for distinct real roots. The Wronskian gives an algorithmic, rigorous test: if it is nonzero at even one point, the pair is guaranteed linearly independent everywhere (by Abel's identity, since the Wronskian of two solutions of a linear equation either vanishes everywhere or nowhere). This removes any need to guess or visually confirm independence.

How do complex and repeated roots interact with initial value problems compared to the distinct real root case?

The mechanics are identical in spirit: write the general two-parameter solution from the fundamental set, then solve a 2×2 linear system for C1,C2C_1, C_2 using y(t0)=y0y(t_0)=y_0 and y′(t0)=y0′y'(t_0)=y_0'. The system is always solvable uniquely because the Wronskian at t0t_0 is nonzero (guaranteed by the theorem's proof), regardless of whether the roots are real, complex, or repeated. The only extra bookkeeping is more careful differentiation — using the product rule for tertte^{rt}, or the product-and-chain rule for eαtcos⁡βte^{\alpha t}\cos\beta t — which is exactly what both worked examples demonstrate explicitly.

Does the sign convention for β\beta (choosing +⋯+\sqrt{\cdots} vs. −⋯-\sqrt{\cdots}) matter for the final answer?

No. If one instead labels the roots α±iβ\alpha \pm i\beta with β<0\beta<0, the resulting basis eαtcos⁡βt,eαtsin⁡βte^{\alpha t}\cos\beta t, e^{\alpha t}\sin\beta t spans exactly the same two-dimensional space, since cos⁡\cos is even and sin⁡\sin is odd — the effect is just to swap the sign of C2C_2 in the general solution. By convention this article takes β>0\beta>0 to avoid ambiguity, but the solution space, and hence every physical prediction, is identical either way.

Before this

Where this leads