evan's notes
cal 3discreteodeenv systems

On this page

  • Homogeneous equations with constant coefficients
  • Why exponential trial functions work
  • Example: two distinct real roots
  • The principle of superposition
  • A complete initial-value problem
  • Existence and uniqueness on an interval
  • The Wronskian and fundamental solutions

Ordinary Differential Equations

§3 Second-Order Linear Equations

Evan Luo · Sep 23, 2026

Ordinary Differential Equations

§3 Second-Order Linear Equations

Evan LuoToday

6 min read

A second-order differential equation has y′′y''y′′ as its highest derivative; it may also contain y′y'y′ and yyy. For a linear equation, the normalized form is

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

Here ppp, qqq, and ggg are known functions of the independent variable ttt. The equation is homogeneous when g(t)=0g(t)=0g(t)=0 throughout the interval under consideration:

y′′+p(t)y′+q(t)y=0.y''+p(t)y'+q(t)y=0. y′′+p(t)y′+q(t)y=0.

If ggg is not identically zero, the equation is nonhomogeneous. The function ggg is the input or forcing term.

Because the equation is second order, an initial-value problem normally needs two pieces of information at the same time t0t_0t0​:

y(t0)=y0,y′(t0)=y0′.y(t_0)=y_0, \qquad y'(t_0)=y'_0. y(t0​)=y0​,y′(t0​)=y0′​.

The first condition fixes the starting value. The second fixes the starting slope.

Homogeneous equations with constant coefficients

Start with

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

where aaa, bbb, and ccc are constants. We assume a≠0a\ne0a=0 so that this is genuinely a second-order equation. If a=0a=0a=0 but b≠0b\ne0b=0, it reduces to the first-order equation

y′=−cby,y'=-\frac cb y, y′=−bc​y,

whose solutions are y(t)=Ce−(c/b)ty(t)=Ce^{-(c/b)t}y(t)=Ce−(c/b)t.

Why exponential trial functions work

For a≠0a\ne0a=0, we want a function whose derivatives keep the same shape, so all three terms can share a common factor. Exponentials do exactly that. Try

y(t)=ert,y(t)=e^{rt}, y(t)=ert,

where rrr is a constant that we still need to find. Exponentials are useful here because differentiating them only introduces powers of rrr:

y′(t)=rert,y′′(t)=r2ert.y'(t)=re^{rt}, \qquad y''(t)=r^2e^{rt}. y′(t)=rert,y′′(t)=r2ert.

Put these expressions into the differential equation:

ar2ert+brert+cert=0.ar^2e^{rt}+bre^{rt}+ce^{rt}=0. ar2ert+brert+cert=0.

Factor out the common exponential:

ert(ar2+br+c)=0.e^{rt}\bigl(ar^2+br+c\bigr)=0. ert(ar2+br+c)=0.

The exponential erte^{rt}ert is never zero, so the polynomial factor must vanish:

ar2+br+c=0.\boxed{ar^2+br+c=0}. ar2+br+c=0​.

The polynomial

P(r)=ar2+br+cP(r)=ar^2+br+c P(r)=ar2+br+c

is the characteristic polynomial, and P(r)=0P(r)=0P(r)=0 is the characteristic equation. Each real root gives a real exponential solution. Repeated or non-real roots require an extra step, so the first example uses two distinct real roots.

Example: two distinct real roots

Consider

y′′+5y′+6y=0.y''+5y'+6y=0. y′′+5y′+6y=0.

Its characteristic equation is

r2+5r+6=0.r^2+5r+6=0. r2+5r+6=0.

Factor it:

(r+2)(r+3)=0.(r+2)(r+3)=0. (r+2)(r+3)=0.

The roots are

r1=−2,r2=−3,r_1=-2, \qquad r_2=-3, r1​=−2,r2​=−3,

so we get two solutions:

y1(t)=e−2t,y2(t)=e−3t.y_1(t)=e^{-2t}, \qquad y_2(t)=e^{-3t}. y1​(t)=e−2t,y2​(t)=e−3t.

The principle of superposition

For the homogeneous equation

y′′+p(t)y′+q(t)y=0,y''+p(t)y'+q(t)y=0, y′′+p(t)y′+q(t)y=0,

any linear combination of solutions is another solution. This is the principle of superposition: if y1y_1y1​ and y2y_2y2​ solve the equation, then

y(t)=C1y1(t)+C2y2(t)y(t)=C_1y_1(t)+C_2y_2(t) y(t)=C1​y1​(t)+C2​y2​(t)

also solves it for every C1,C2∈RC_1,C_2\in\mathbb RC1​,C2​∈R.

The reason is the linearity of differentiation. Define the differential operator

L[y]=y′′+p(t)y′+q(t)y.L[y]=y''+p(t)y'+q(t)y. L[y]=y′′+p(t)y′+q(t)y.

Then

L[C1y1+C2y2]=C1L[y1]+C2L[y2]=C1(0)+C2(0)=0.\begin{aligned} L[C_1y_1+C_2y_2] &=C_1L[y_1]+C_2L[y_2]\\ &=C_1(0)+C_2(0)\\ &=0. \end{aligned} L[C1​y1​+C2​y2​]​=C1​L[y1​]+C2​L[y2​]=C1​(0)+C2​(0)=0.​

This argument depends on the right side being zero. If L[y1]=gL[y_1]=gL[y1​]=g and L[y2]=gL[y_2]=gL[y2​]=g, then

L[C1y1+C2y2]=(C1+C2)g,L[C_1y_1+C_2y_2]=(C_1+C_2)g, L[C1​y1​+C2​y2​]=(C1​+C2​)g,

which is not generally equal to ggg. Arbitrary superposition is therefore a property of the associated homogeneous equation, not of a nonhomogeneous equation with fixed forcing.

A complete initial-value problem

Return to

y′′+5y′+6y=0.y''+5y'+6y=0. y′′+5y′+6y=0.

Superposition gives the family

y(t)=C1e−2t+C2e−3t.y(t)=C_1e^{-2t}+C_2e^{-3t}. y(t)=C1​e−2t+C2​e−3t.

Now impose

y(0)=2,y′(0)=3.y(0)=2, \qquad y'(0)=3. y(0)=2,y′(0)=3.

Evaluating the solution and its derivative at t=0t=0t=0 gives

C1+C2=2,C_1+C_2=2, C1​+C2​=2,

and

−2C1−3C2=3.-2C_1-3C_2=3. −2C1​−3C2​=3.

The same system can be written as

(11−2−3)(C1C2)=(23).\begin{pmatrix} 1&1\\ -2&-3 \end{pmatrix} \begin{pmatrix} C_1\\ C_2 \end{pmatrix} = \begin{pmatrix} 2\\ 3 \end{pmatrix}. (1−2​1−3​)(C1​C2​​)=(23​).

Its determinant is −1-1−1, so the system has one solution:

C1=9,C2=−7.C_1=9, \qquad C_2=-7. C1​=9,C2​=−7.

Therefore the unique solution of the IVP is

y(t)=9e−2t−7e−3t.\boxed{y(t)=9e^{-2t}-7e^{-3t}}. y(t)=9e−2t−7e−3t​.

A quick check confirms both initial conditions:

y(0)=9−7=2,y′(t)=−18e−2t+21e−3t,y′(0)=−18+21=3.\begin{aligned} y(0)&=9-7=2,\\ y'(t)&=-18e^{-2t}+21e^{-3t},\\ y'(0)&=-18+21=3. \end{aligned} y(0)y′(t)y′(0)​=9−7=2,=−18e−2t+21e−3t,=−18+21=3.​

The initial slope lifts the solution before both modes decay

Solution of the second-order initial-value problemThe solution starts at height two with initial slope three, rises to one maximum, and then approaches zero from above as time increases.y(0)=2y′(0)=3y(t)=9e−2t−7e−3tty
The curve y(t)=9e−2t−7e−3t starts at y(0)=2 and initially rises because y′(0)=3. Both exponential modes decay, so the solution eventually approaches zero from above.

The positive initial slope explains why the curve rises at first. Both exponential terms decay as ttt increases, so the solution eventually approaches zero.

Existence and uniqueness on an interval

Consider the linear IVP

y′′+p(t)y′+q(t)y=g(t),y(t0)=y0,y′(t0)=y0′.\begin{aligned} y''+p(t)y'+q(t)y&=g(t),\\ y(t_0)&=y_0,\\ y'(t_0)&=y'_0. \end{aligned} y′′+p(t)y′+q(t)yy(t0​)y′(t0​)​=g(t),=y0​,=y0′​.​

Let III be an open interval with t0∈It_0\in It0​∈I. If ppp, qqq, and ggg are continuous on the interval III, then the IVP has exactly one solution on the entire interval III:

y:I→R.y:I\to\mathbb R. y:I→R.

This is stronger than a merely local statement. The solution does not stop at an interior point where the normalized coefficients remain continuous. The interval ends only where the chosen coefficient interval ends; it need not be all of R\mathbb RR.

If an equation is initially written as

a2(t)y′′+a1(t)y′+a0(t)y=h(t),a_2(t)y''+a_1(t)y'+a_0(t)y=h(t), a2​(t)y′′+a1​(t)y′+a0​(t)y=h(t),

first divide by a2(t)a_2(t)a2​(t). The normalized coefficients are

p(t)=a1(t)a2(t),q(t)=a0(t)a2(t),g(t)=h(t)a2(t).p(t)=\frac{a_1(t)}{a_2(t)}, \qquad q(t)=\frac{a_0(t)}{a_2(t)}, \qquad g(t)=\frac{h(t)}{a_2(t)}. p(t)=a2​(t)a1​(t)​,q(t)=a2​(t)a0​(t)​,g(t)=a2​(t)h(t)​.

The theorem applies on an interval where a2(t)≠0a_2(t)\ne0a2​(t)=0 and these three quotients are continuous. A zero of a2a_2a2​ can split the real line into separate coefficient intervals.

The Wronskian and fundamental solutions

The Wronskian answers a concrete question: do two solutions provide enough freedom to match any starting value and slope? Given two differentiable functions fff and ggg, their Wronskian is

W[f,g](t)=det⁡(f(t)g(t)f′(t)g′(t))=f(t)g′(t)−f′(t)g(t).\begin{aligned} W[f,g](t) &=\det \begin{pmatrix} f(t)&g(t)\\ f'(t)&g'(t) \end{pmatrix}\\ &=f(t)g'(t)-f'(t)g(t). \end{aligned} W[f,g](t)​=det(f(t)f′(t)​g(t)g′(t)​)=f(t)g′(t)−f′(t)g(t).​

The Wronskian tells us whether two solution functions provide independent initial data. Suppose y1y_1y1​ and y2y_2y2​ solve the same homogeneous second-order linear equation on III, with ppp and qqq continuous there. To make

y=C1y1+C2y2y=C_1y_1+C_2y_2 y=C1​y1​+C2​y2​

satisfy prescribed values at t0t_0t0​, the constants must solve

(y1(t0)y2(t0)y1′(t0)y2′(t0))(C1C2)=(y0y0′).\begin{pmatrix} y_1(t_0)&y_2(t_0)\\ y_1'(t_0)&y_2'(t_0) \end{pmatrix} \begin{pmatrix} C_1\\ C_2 \end{pmatrix} = \begin{pmatrix} y_0\\ y'_0 \end{pmatrix}. (y1​(t0​)y1′​(t0​)​y2​(t0​)y2′​(t0​)​)(C1​C2​​)=(y0​y0′​​).

The determinant of the coefficient matrix is exactly

W[y1,y2](t0).W[y_1,y_2](t_0). W[y1​,y2​](t0​).

The matrix can match every pair of initial values if and only if

W[y1,y2](t0)≠0.\boxed{W[y_1,y_2](t_0)\ne0}. W[y1​,y2​](t0​)=0​.

When this condition holds, the matrix is invertible, so every pair of initial values determines one pair of constants C1,C2C_1,C_2C1​,C2​. The existence-and-uniqueness theorem then shows that this linear combination is the only solution with those initial values.

In that case, {y1,y2}\{y_1,y_2\}{y1​,y2​} is called a fundamental set of solutions on III, and

y(t)=C1y1(t)+C2y2(t)\boxed{y(t)=C_1y_1(t)+C_2y_2(t)} y(t)=C1​y1​(t)+C2​y2​(t)​

is the general solution of the homogeneous equation: every solution on III appears exactly once for a choice of C1C_1C1​ and C2C_2C2​.

For the example y1=e−2ty_1=e^{-2t}y1​=e−2t and y2=e−3ty_2=e^{-3t}y2​=e−3t,

W[y1,y2](t)=e−2t(−3e−3t)−(−2e−2t)e−3t=−e−5t.\begin{aligned} W[y_1,y_2](t) &=e^{-2t}\bigl(-3e^{-3t}\bigr) -\bigl(-2e^{-2t}\bigr)e^{-3t}\\ &=-e^{-5t}. \end{aligned} W[y1​,y2​](t)​=e−2t(−3e−3t)−(−2e−2t)e−3t=−e−5t.​

This never vanishes, so the two exponentials form a fundamental set of solutions on R\mathbb RR.

Source: https://notes.ohevan.com/notes/ode/03-second-order-linear-equations

© 2026 Evan Luo. All rights reserved.

Back to Ordinary Differential Equations

0 reads

  • stay up to date

  • about me

  • coffee

© 2026 Evan Luo. All rights reserved.