Unit 6, lesson 45

Second-Order Linear Equations

Why a second-order linear homogeneous equation has a two-parameter family of solutions, how the Wronskian tests independence, and how the characteristic roots give every solution when the coefficients are constant.

  • Theory part 89
  • Practice part 90
  • Second-Order Linear Equations

A mass on a spring is described by its position yy, its velocity y′y' and its acceleration y′′y'', and many such systems obey a law tying the three together linearly. Knowing where the mass starts and how fast it is moving picks out exactly one motion. This lesson explains why two numbers are always enough, and how to write down every solution when the coefficients are constants.

Key terms

Linear homogeneous equation

Definition A second-order linear differential equation has the form

P(t) y′′+Q(t) y′+R(t) y=G(t),P(t)\,y'' + Q(t)\,y' + R(t)\,y = G(t),

with PP, QQ, RR and GG functions of tt alone. It is homogeneous when G(t)=0G(t) = 0 for all tt. On an interval where P(t)≠0P(t) \ne 0 we may divide by PP and get the normalized form

y′′+p(t) y′+q(t) y=g(t).y'' + p(t)\,y' + q(t)\,y = g(t).

Points where P(t)=0P(t) = 0 are singular points, and the theory below is applied on intervals that avoid them.

In plain words Linear means yy, y′y' and y′′y'' appear only to the first power, multiplied by functions of tt and added. Homogeneous means nothing is left over on the right side.

Example The equation t2y′′+ty′−y=0t^2 y'' + t y' - y = 0 is linear and homogeneous. Its leading coefficient vanishes at t=0t = 0, so on the interval t>0t > 0 we normalize to

y′′+1t y′−1t2 y=0.y'' + \frac{1}{t}\,y' - \frac{1}{t^2}\,y = 0.

The function y=ty = t solves it, since 0+t−t=00 + t - t = 0 after multiplying back by t2t^2, and so does y=1/ty = 1/t, since 2/t+(−1/t)−1/t=02/t + (-1/t) - 1/t = 0 after the same step.

Continuous coefficients and existence

Definition Let pp, qq and gg be continuous on an open interval II, and let t0t_0 be in II. For any real numbers y0y_0 and y1y_1, the initial value problem

y′′+p(t) y′+q(t) y=g(t),y(t0)=y0,y′(t0)=y1\begin{gathered} y'' + p(t)\,y' + q(t)\,y = g(t), \\[0.4em] y(t_0) = y_0, \qquad y'(t_0) = y_1 \end{gathered}

has exactly one solution, and that solution exists on the whole interval II.

In plain words As long as the normalized coefficients are continuous, a starting position and a starting slope pin down one solution, and it lives as long as the coefficients stay continuous.

Example For (t−3) y′′+t y′+(ln⁡t) y=0(t - 3)\,y'' + t\,y' + (\ln t)\,y = 0 with y(1)=2y(1) = 2 and y′(1)=1y'(1) = 1, the normalized coefficients are

p(t)=tt−3,q(t)=ln⁡tt−3,p(t) = \frac{t}{t - 3}, \qquad q(t) = \frac{\ln t}{t - 3},

which are continuous on (0,3)(0, 3), the largest open interval around t0=1t_0 = 1 that avoids t=0t = 0 and t=3t = 3. The theorem guarantees a unique solution on all of (0,3)(0, 3), even though we cannot write it in closed form.

Superposition

Definition If y1y_1 and y2y_2 solve the homogeneous equation y′′+p(t) y′+q(t) y=0y'' + p(t)\,y' + q(t)\,y = 0 on an interval II, then so does c1y1+c2y2c_1 y_1 + c_2 y_2 for all constants c1c_1 and c2c_2. This follows from linearity of differentiation. The principle fails for nonhomogeneous equations, since there the right side would be multiplied by c1+c2c_1 + c_2.

In plain words Add solutions, or scale them, and you still have a solution. The solutions form a vector space.

Example Both ete^t and e−te^{-t} solve y′′−y=0y'' - y = 0, so cosh⁡t=12et+12e−t\cosh t = \frac{1}{2}e^t + \frac{1}{2}e^{-t} and 3et−5e−t3e^t - 5e^{-t} do as well.

Independent solutions

Definition Two functions y1y_1 and y2y_2 are linearly independent on an interval II if the only constants with c1y1(t)+c2y2(t)=0c_1 y_1(t) + c_2 y_2(t) = 0 for every tt in II are c1=c2=0c_1 = c_2 = 0. Equivalently, neither is a constant multiple of the other on II. If pp and qq are continuous on II and y1y_1, y2y_2 are linearly independent solutions of the homogeneous equation on II, they form a fundamental set, and every solution on II has the form

y=c1y1+c2y2.y = c_1 y_1 + c_2 y_2.

In plain words Two genuinely different solutions are enough to build all the others. A solution and a multiple of it are not enough.

Example On t>0t > 0 the solutions tt and 1/t1/t of t2y′′+ty′−y=0t^2 y'' + t y' - y = 0 are independent, because t/(1/t)=t2t/(1/t) = t^2 is not constant. So every solution on t>0t > 0 is c1t+c2/tc_1 t + c_2/t. By contrast, ete^t and 4et4e^t both solve y′′−y=0y'' - y = 0 but are dependent, and their combinations miss e−te^{-t}.

The Wronskian

Definition The Wronskian of two differentiable functions is

W(y1,y2)(t)=∣y1y2y1′y2′∣=y1y2′−y1′y2.\begin{aligned} W(y_1, y_2)(t) &= \begin{vmatrix} y_1 & y_2 \\ y_1' & y_2' \end{vmatrix} \\[0.4em] &= y_1 y_2' - y_1' y_2. \end{aligned}

If y1y_1 and y2y_2 solve y′′+p(t) y′+q(t) y=0y'' + p(t)\,y' + q(t)\,y = 0 with pp and qq continuous on II, then Abel’s formula holds,

W(t)=W(t0) e−∫t0tp(s) ds,W(t) = W(t_0)\,e^{-\int_{t_0}^{t} p(s)\,ds},

so WW is either zero at every point of II or zero at no point of II. The two solutions are independent on II exactly when W≠0W \ne 0. For any two differentiable functions, W≠0W \ne 0 at a single point already forces independence. The converse needs the functions to be solutions. The functions t2t^2 and t∣t∣t \lvert t \rvert on R\mathbb{R} have W=0W = 0 everywhere, since the derivative of t∣t∣t \lvert t \rvert is 2∣t∣2\lvert t \rvert, yet they are independent on R\mathbb{R}, because t∣t∣t \lvert t \rvert equals t2t^2 for t>0t > 0 and −t2-t^2 for t<0t < 0. So they cannot both solve one equation of this form with pp and qq continuous on R\mathbb{R}.

In plain words The Wronskian is a single number at each tt that tells you whether the vectors (y1,y1′)(y_1, y_1') and (y2,y2′)(y_2, y_2') point in different directions. For solutions it is all or nothing.

Example For y1=ty_1 = t and y2=1/ty_2 = 1/t on t>0t > 0,

W=t⋅(−1t2)−1⋅1t=−2t,W = t \cdot \Bigl(-\frac{1}{t^2}\Bigr) - 1 \cdot \frac{1}{t} = -\frac{2}{t},

never zero, as expected. It also matches Abel’s formula, since p=1/tp = 1/t gives e−ln⁡t=1/te^{-\ln t} = 1/t.

The characteristic equation

Definition For the equation ay′′+by′+cy=0a y'' + b y' + c y = 0 with real constants aa, bb, cc and a≠0a \ne 0, the characteristic equation is

ar2+br+c=0.a r^2 + b r + c = 0.

Substituting y=erty = e^{rt} gives (ar2+br+c) ert=0(a r^2 + b r + c)\,e^{rt} = 0, and since ert≠0e^{rt} \ne 0, the exponential erte^{rt} is a solution exactly when rr is a root. The coefficients are constant, hence continuous on all of R\mathbb{R}, so every solution is defined for all real tt.

In plain words Guess an exponential, and the differential equation becomes a quadratic equation for the rate rr. The sign of b2−4acb^2 - 4ac decides which of three cases you are in.

Example For y′′−5y′+6y=0y'' - 5y' + 6y = 0 the characteristic equation r2−5r+6=0r^2 - 5r + 6 = 0 factors as (r−2)(r−3)=0(r - 2)(r - 3) = 0. So e2te^{2t} and e3te^{3t} are solutions.

Distinct and repeated real roots

Definition Let a≠0a \ne 0, bb, cc be real.

  • If b2−4ac>0b^2 - 4ac > 0, the roots r1≠r2r_1 \ne r_2 are real and the general solution is y=c1er1t+c2er2ty = c_1 e^{r_1 t} + c_2 e^{r_2 t}. The Wronskian is (r2−r1) e(r1+r2)t≠0(r_2 - r_1)\,e^{(r_1 + r_2)t} \ne 0.
  • If b2−4ac=0b^2 - 4ac = 0, there is one root r=−b/(2a)r = -b/(2a), and the general solution is y=(c1+c2t) erty = (c_1 + c_2 t)\,e^{rt}. The second solution tertt e^{rt} works because 2ar+b=02ar + b = 0, and the Wronskian is e2rt≠0e^{2rt} \ne 0.

In plain words Two different real rates give two exponentials. When the rates collide, the missing second solution is the same exponential times tt.

Example For y′′+4y′+4y=0y'' + 4y' + 4y = 0 with y(0)=1y(0) = 1 and y′(0)=0y'(0) = 0, the only root is r=−2r = -2, so y=(c1+c2t) e−2ty = (c_1 + c_2 t)\,e^{-2t}. Then y(0)=c1=1y(0) = c_1 = 1 and y′(0)=c2−2c1=0y'(0) = c_2 - 2c_1 = 0 give c2=2c_2 = 2, and

y=(1+2t) e−2t.y = (1 + 2t)\,e^{-2t}.

Its derivative is y′=−4t e−2ty' = -4t\,e^{-2t}, which is 00 at t=0t = 0 as required.

Conjugate complex roots

Definition If b2−4ac<0b^2 - 4ac < 0, the roots are r=α±iβr = \alpha \pm i\beta with

α=−b2a,β=4ac−b22∣a∣>0.\begin{gathered} \alpha = -\frac{b}{2a}, \\[0.4em] \beta = \frac{\sqrt{4ac - b^2}}{2\lvert a \rvert} > 0. \end{gathered}

Since the coefficients are real, the real and imaginary parts of the complex solution e(α+iβ)te^{(\alpha + i\beta)t} are real solutions, and the general real solution is

y=eαt(c1cos⁡βt+c2sin⁡βt),y = e^{\alpha t}(c_1 \cos \beta t + c_2 \sin \beta t),

with Wronskian βe2αt≠0\beta e^{2\alpha t} \ne 0.

In plain words Complex roots mean oscillation. The imaginary part β\beta sets the angular frequency and the real part α\alpha makes the swings grow or decay.

Example For y′′+2y′+10y=0y'' + 2y' + 10y = 0 the roots of r2+2r+10=0r^2 + 2r + 10 = 0 are r=−1±3ir = -1 \pm 3i. Every solution is

y=e−t(c1cos⁡3t+c2sin⁡3t),y = e^{-t}(c_1 \cos 3t + c_2 \sin 3t),

an oscillation with angular frequency 33 whose amplitude decays like e−te^{-t}.

Initial value and boundary value problems

Definition An initial value problem prescribes y(t0)y(t_0) and y′(t0)y'(t_0) at one point, and under the continuity hypotheses above it always has exactly one solution. A boundary value problem prescribes values at two different points, such as y(t0)=Ay(t_0) = A and y(t1)=By(t_1) = B with t0≠t1t_0 \ne t_1. A boundary value problem may have exactly one solution, no solution, or infinitely many, and the general solution must be checked against both conditions.

In plain words Fixing position and slope at one moment always works. Fixing positions at two moments is a different kind of question, and the answer depends on the numbers.

Example For y′′+y=0y'' + y = 0 every solution is c1cos⁡t+c2sin⁡tc_1 \cos t + c_2 \sin t, and y(0)=0y(0) = 0 forces c1=0c_1 = 0.

Second conditionResult
y(π/2)=1y(\pi/2) = 1one solution, y=sin⁡ty = \sin t
y(π)=0y(\pi) = 0infinitely many, y=c2sin⁡ty = c_2 \sin t
y(π)=1y(\pi) = 1none, since sin⁡π=0\sin \pi = 0

Check yourself

  1. Solve y′′−y′−6y=0y'' - y' - 6y = 0 with y(0)=3y(0) = 3 and y′(0)=−1y'(0) = -1.
Show the answer

The roots of r2−r−6=0r^2 - r - 6 = 0 are 33 and −2-2, so y=c1e3t+c2e−2ty = c_1 e^{3t} + c_2 e^{-2t}. The conditions give c1+c2=3c_1 + c_2 = 3 and 3c1−2c2=−13c_1 - 2c_2 = -1, so c1=1c_1 = 1 and c2=2c_2 = 2. The solution is

y=e3t+2e−2t.y = e^{3t} + 2e^{-2t}.
  1. Compute the Wronskian of e−tcos⁡3te^{-t}\cos 3t and e−tsin⁡3te^{-t}\sin 3t and decide whether they form a fundamental set for y′′+2y′+10y=0y'' + 2y' + 10y = 0.
Show the answer

A direct computation gives W=3e−2tW = 3e^{-2t}, which is never zero. Both functions are solutions, so they form a fundamental set on R\mathbb{R}. This agrees with the formula βe2αt\beta e^{2\alpha t} with α=−1\alpha = -1 and β=3\beta = 3.

  1. For y′′+4y=0y'' + 4y = 0 with y(0)=1y(0) = 1, decide how many solutions satisfy y(π/2)=2y(\pi/2) = 2, and how many satisfy y(π/2)=−1y(\pi/2) = -1.
Show the answer

Every solution is c1cos⁡2t+c2sin⁡2tc_1 \cos 2t + c_2 \sin 2t, and y(0)=1y(0) = 1 gives c1=1c_1 = 1. Then y(π/2)=cos⁡π+c2sin⁡π=−1y(\pi/2) = \cos \pi + c_2 \sin \pi = -1 whatever c2c_2 is. So y(π/2)=2y(\pi/2) = 2 has no solution, and y(π/2)=−1y(\pi/2) = -1 has infinitely many, namely cos⁡2t+c2sin⁡2t\cos 2t + c_2 \sin 2t.

How this lesson connects

It also uses Ordinary derivatives and the product rule, Exponential functions, sine and cosine and Complex conjugate roots and two by two determinants.