Chapter2 HigherOrderODE Textbook
Chapter2 HigherOrderODE Textbook
Chapter 2
Higher Order Linear
Differential Equations
Topics Covered:
3
4 Applied Mathematics II — TCET / MU
Strategic Note
This chapter is the single highest-mark chapter in Applied Mathematics II. VPM
(Variation of Parameters) alone appears for 8 marks in every paper without
exception (2022–2025). Master the 6 VPM patterns in Section 1.7 and you guarantee
8 marks instantly. CF + PI accounts for another 6–8 marks per paper. Total potential
from this chapter: 20–24 marks.
Complete Solution
The general solution of f (D)y = X is:
d
y = yc + yp where D≡
|{z} |{z} dx
C.F. P.I.
5
6 Applied Mathematics II — TCET / MU
AE Roots CF
PI Formula Reference
Type X= Formula
1 0 P.I. = 0
eax xeax x2 eax
2 eax if f (a) ̸= 0; ′ if f (a) = 0; ′′
f (a) f (a) f (a)
′
if f (a) = 0
(1 + u)−1 = 1 − u + u2 − u3 + · · · (1 − u)−1 = 1 + u + u2 + u3 + · · ·
(1 + u)−2 = 1 − 2u + 3u2 − 4u3 + · · · (1 − u)−2 = 1 + 2u + 3u2 + 4u3 + · · ·
Rule: When operating on xm , stop expansion at Dm (higher derivatives vanish).
x
Must use VPM when: X = tan ax, sec ax, csc ax, cot ax, ee , log x, 1/(1 +
ex ), sin(ex ).
Standard PI types fail for these — VPM always works.
Variable-Coefficient Equations
Cauchy–Euler: a0 x2 y ′′ + a1 xy ′ + a2 y = X
dy d2 y
Substitution: x = ez (z = ln x). Then: x = Dy, x2 2 = D(D − 1)y,
dx dx
d3 y
x3 = D(D − 1)(D − 2)y
dx3
Legendre’s: a0 (ax + b)2 y ′′ + a1 (ax + b)y ′ + a2 y = X
dy d2 y
Substitution: ax + b = ez . Then: (ax + b) = a Dy, (ax + b)2 2 = a2 D(D − 1)y
dx dx
DE W u, v P.I.
x x
y ′′ + 3y ′ + 2y = ee −e−3x − cos(ex ); e−2x · ee
ex cos(ex )−
sin(ex )
y ′′ + y = 1
1+sin x 1 tan x − −1 − x cos x + sin x ln(1 +
x − sec x; sin x)
ln(1 +
sin x)
tan x dx = ln | sec x|
R R
x sin x dx = −x cos x + sin x
sec x dx = ln | sec x + tan x|
R R
x cos x dx = x sin x + cos x
csc x dx = ln | csc x − cot x|
R R R
u dv = uv − v du (IBP)
cot x dx = ln | sin x|
R R
t et dt = et (t − 1) + C
x
R
ex sin x dx = e (sin x−cos x) t sin t dt = sin t − t cos t + C
R
2
x
R
ex cos x dx = e (sin x+cos
R x) t cos t dt = cos t + t sin t + C
2
1.2 Theory
1.2.1 Why CF + PI?
The general solution decomposes into two parts because the ODE is linear. The CF
captures the homogeneous (natural) response — how the system behaves freely. The
PI captures the forced response to X. In engineering: CF ≈ transient response; PI ≈
steady-state response.
Solution
A.E. : m3 + 6m2 + 11m + 6 = 0. Try m = −1: −1 + 6 − 11 + 6 = 0 ✓. Factor:
(m + 1)(m2 + 5m + 6) = (m + 1)(m + 2)(m + 3) = 0.
Solution
A.E. : m2 − 4m + 4 = (m − 2)2 = 0 ⇒ m = 2, 2 (repeated).
Solution
A.E. : m3 − 3m2 + 3m − 1 = (m − 1)3 = 0 ⇒ m = 1, 1, 1 (triple).
Solution
√
−4 ± 16 − 20
A.E. : m2 + 4m + 5 = 0 ⇒ m = = −2 ± i.
2
α = −2, β = 1.
C.S. = e−2x (c1 cos x + c2 sin x)
Solution
A.E. : (m − 2)2 (m + 1) = 0 ⇒ m = 2, 2, −1.
Solution
(D4 + 8D2 + 16) = (D2 + 4)2 = 0 ⇒ m2 = −4 ⇒ m = ±2i (each repeated twice).
α = 0, β = 2, double complex roots.
Solution
Let u = D2 : u2 − 5u + 4 = (u − 1)(u − 4) = 0 ⇒ D2 = 1 or D2 = 4.
D2 = 1 ⇒ m = ±1; D2 = 4 ⇒ m = ±2.
Solution
f (D) = D2 − 3D + 2. f (3) = 9 − 9 + 2 = 2 ̸= 0.
e3x e3x
P.I. = =
f (3) 2
Solution
A.E. : (m − 2)2 = 0 ⇒ m = 2, 2. C.F. = (c1 + c2 x)e2x .
f (D) = (D − 2)2 . f (2) = 0 (failure). f ′ (D) = 2(D − 2), f ′ (2) = 0 (failure again).
f ′′ (D) = 2, f ′′ (2) = 2.
x2 e2x x2 e2x
P.I. = ′′ =
f (2) 2
2x x2 e2x
C.S. = (c1 + c2 x)e +
2
Solution
A.E. : (m + 1)(m + 2) = 0 ⇒ m = −1, −2. C.F. = c1 e−x + c2 e−2x .
f (−1) = (1 − 3 + 2) = 0 (failure). f ′ (D) = 2D + 3, f ′ (−1) = 1 ̸= 0.
x e−x
P.I. = ′
= x e−x
f (−1)
Solution
2 cosh(2x) = e2x + e−2x . A.E. : m2 − 4 = 0 ⇒ m = ±2. C.F. = c1 e2x + c2 e−2x .
′ ′ x e2x
f (2) = 0 (failure); f (D) = 2D, f (2) = 4. P.I.1 = .
4
−2x
xe
f (−2) = 0 (failure); f ′ (−2) = −4. P.I.2 = .
−4
x e2x x e−2x x(e2x − e−2x ) x sinh(2x)
P.I. = − = =
4 4 4 2
x sinh(2x)
C.S. = c1 e2x + c2 e−2x +
2
Solution
Write 5 = 5 e0·x . f (0) = 0 + 0 + 2 = 2.
5
P.I. =
2
Solution
A.E. : (m − 3)(m + 2) = 0 ⇒ m = 3, −2. C.F. = c1 e3x + c2 e−2x .
f (−2) = 4 + 2 − 6 = 0 (failure). f ′ (D) = 2D − 1, f ′ (−2) = −5.
x e−2x
P.I. =
−5
x e−2x
C.S. = c1 e3x + c2 e−2x −
5
Solution
A.E. : m = 2, 3. C.F. = c1 e2x + c2 e3x .
1
P.I. = 2 sin 3x. Replace D2 → −9:
D − 5D + 6
sin 3x sin 3x
= =
−9 − 5D + 6 −3 − 5D
Rationalize with (−3 + 5D):
5 cos 3x − sin 3x
C.S. = c1 e2x + c2 e3x +
78
Solution
A.E. : m = ±2i. C.F. = c1 cos 2x + c2 sin 2x.
Replace D2 → −4: denominator = −4 + 4 = 0 (failure).
Use failure formula with a = 2:
1 x sin 2x x sin 2x
cos 2x = =
D2 +4 2·2 4
x sin 2x
C.S. = c1 cos 2x + c2 sin 2x +
4
Solution
C.F. = c1 cos x + c2 sin x.
D2 → −1: denom = −1 + 1 = 0 (failure). Use:
1 −x cos x −x cos x
sin x = =
D2 +1 2·1 2
x cos x
C.S. = c1 cos x + c2 sin x −
2
Solution
Use sin A sin B = 21 [cos(A − B) − cos(A + B)]:
x sin x cos 3x
C.S. = c1 cos x + c2 sin x + +
4 16
Solution
1+cos 2x
cos2 x = 2
.
P.I. = P.I.1 + P.I.2 .
1 1 1 1 1
P.I.1 = 2 · = = = .
D + 3D + 2 2 2f (0) 2·2 4
1 cos 2x
P.I.2 = 2 · . Replace D2 → −4:
D + 3D + 2 2
cos 2x cos 2x (3D + 2) cos 2x −6 sin 2x + 2 cos 2x −6 sin 2x + 2 co
= = = = =
2(−4 + 3D + 2) 2(3D − 2) 2(9D2 − 4) D2 →−4 2(−36 − 4) −80
1 3 sin 2x − cos 2x
P.I. = +
4 40
Solution
Replace D2 → −4 and D3 = D · D2 → −4D:
f (D) D2 →−4
= −4D − (−4) + 4D − 4 = 4 − 4 = 0 (failure).
1
P.I. = (−x cos 2x + 2x sin 2x) + · · ·
20
(In exams, stop at this recognition level and cite the method.)
1.6.1 Type 4: X = xm
Solution
1 2 1
P.I. = x = x2 .
(D − 2)2 4(1 − D/2)2
(1 − D/2)−2 = 1 + D + 43 D2 + · · · Operate on x2 :
= x2 + D(x2 ) + 43 D2 (x2 ) = x2 + 2x + 43 · 2 = x2 + 2x + 3
2
x2 + 2x + 3/2 x2 x 3
P.I. = = + +
4 4 2 8
Solution
1 1 1
P.I. = 3
(3x2 − 5x) = · 3
(3x2 − 5x).
4 + D − 2D 4 1 + (D − 2D)/4
Let u = (D3 − 2D)/4. (1 + u)−1 ≈ 1 − u + u2 − · · · operating on 3x2 − 5x (need up
to D2 ):
3 3 2 2
−u(3x2 −5x) = − (D −2D)
4
(3x2 −5x) = − D (3x4 −5x) + 2D(3x4 −5x) = 0+ 2(6x−5)
4
= 6x−5
2
=
5
3x − 2 .
2
u2 ≈ (−2D/4)2 = D2 /4: D4 (3x2 − 5x) = 46 = 32 .
2 3x2 − 2x − 1
1
(3x − 5x) + (3x − 52 ) + 32 = 14 (3x2 − 2x − 1) =
P.I. = 4
4
3x2 − 2x − 1
P.I. =
4
Solution
1 1
P.I. = 2
x3 = (1 − D/2)−2 x3 .
(D − 2) 4
2 3
−2
(1 − D/2) = 1 + D + 3D4 + D2 + · · ·
= x3 + 3x2 + 34 (6x) + 12 (6) = x3 + 3x2 + 9x
2
+ 3.
x3 + 3x2 + 9x/2 + 3
P.I. =
4
Solution
1 1
Shift: P.I. = ex · 2
cos x = ex · 2 cos x.
(D + 1) + 2 D + 2D + 3
cos x cos x ex cos x
Replace D2 → −1: P.I. = ex · = ex · = · .
−1 + 2D + 3 2D + 2 2 D+1
1 (D − 1) (D − 1) cos x − sin x − cos x
cos x = 2
cos x D2 →−1 = = =
D+1 D −1 −2 −2
sin x + cos x
.
2
ex (sin x + cos x)
P.I. =
4
Solution
1 1
Shift: P.I. = e3x · 2
x 2
= e 3x
· 2
x2 .
(D + 3) + 2 D + 6D + 11
e3x 1
Factor: = · x2 .
11 1 + (6D + D2 )/11
(1 + u)−1 ≈ 1 − u + u2 where u = (6D + D2 )/11:
2
−u · x2 = − 6D+D11
x2 = − 12x+2 11
2 36D2 2
u2 · x2 = 6D 72
2
11
x = 121
x = 121
Solution
1 1
Shift: P.I. = e2x · sin x = e2x · 2 sin x.
(D + 2)2 − 4(D + 2) + 5 D +1
x e2x cos x
P.I. = −
2
Solution
A.E. : (m − 2)2 = 0 ⇒ m = 2, 2. C.F. = (c1 + c2 x)e2x .
e2x ′ ′′ x2 e2x
P.I.1 = : f (2) = 0, f (2) = 0, f (2) = 2. P.I.1 = .
(D − 2)2 2
cos 3x cos 3x cos 3x
P.I.2 = 2 : D2 → −9: = = . Rationalize:
D − 4D + 4 −9 − 4D + 4 −5 − 4D
(−5 + 4D) cos 3x −5 cos 3x − 12 sin 3x −5 cos 3x − 12 sin 3x
= = =
25 − 16D2 D2 →−9 25 + 144 169
Solution
P.I. = P.I.1 + P.I.2 .
ex (sin x + cos x)
P.I.1 : From Example D4 above: P.I.1 = .
24
3x x 12x 50
P.I.2 : From Example D5 above: P.I.2 = e − + .
11 121 1331
ex (sin x + cos x)
2
3x x 12x 50
P.I. = +e − +
4 11 121 1331
Solution
Shift (a = −1, replace D → D − 1):
ZZ
−x 1
P.I. = e · 2 log x = e−x log x dx dx
D
x2 x2 x2 x2
R R 3
log x dx = x log x − x. (x log x − x) dx = log x − − = log x − .
2 4 2 2 2
x2 e−x
3
P.I. = log x −
2 2
Solution
1 e2x 1
Shift: P.I. = e2x · 2
(1 + x 2
) = − · (1 + x2 ).
D − 3D − 4 4 1 − (3D − D2 )/4
Let u = (3D − D2 )/4: (1 + u)−1 ≈ 1 − u + u2 :
2 2
u(1 + x2 ) = 3D−D
4
(1 + x2 ) = 6x−24
= 3x−12
. u2 (1 + x2 ) ≈ 9D
16
(1 + x2 ) = 9·2
16
= 98 .
e2x e2x 2 3x 21
2 3x − 1 9
P.I. = − (1 + x ) − + =− x − +
4 2 8 4 2 8
2
x 3x 21
P.I. = e2x − + −
4 8 32
x
Example E1 — (D2 + 3D + 2)y = ee (Highest-Frequency PYQ, Dec 2022 Q.1d / Jun 2
Solution
A.E. : (m + 1)(m + 2) = 0 ⇒ m = −1, −2.
x
y1 = e−x , y2 = e−2x , X = ee .
e−x e−2x
W = = −2e−3x + e−3x = −e−3x .
−e−x −2e−2x
x
e−2x · ee
Z Z
x x
dx = ex · ee dx. Let t = ex , dt = ex dx: = et dt = et = ee .
R
u=− −3x
Z −x −eex
e ·e
Z
2x ex
v= dx = − e · e dx. Let t = ex , dt = ex dx, e2x = t2 , dx = dt/t:
−e−3x
Z
x
= − t · et dt = −(tet − et ) = −et (t − 1) = −ee (ex − 1).
x x x x
P.I. = u y1 + v y2 = ee · e−x + [−ee (ex − 1)] · e−2x = ee [e−x − e−x + e−2x ] = e−2x · ee .
x
C.S. = c1 e−x + c2 e−2x + e−2x · ee
Solution
y1 = cos x, y2 = sin x, X = sec x tan x.
P.I. = (x − tan x) cos x + ln | sec x| · sin x = x cos x − sin x + sin x ln | sec x|.
C.S. = c1 cos x + c2 sin x + x cos x − sin x + sin x ln | sec x|
1
Example E3 — (D2 + 1)y = (TCET Sample)
1 + sin x
Solution
1
y1 = cos x, y2 = sin x, W = 1, X = .
Z Z 1 + sinx Z
sin x 1 dx
u=− dx = − 1− dx = −x + .
1 + sinZx 1 + sin x 1 + sin x
1 − sin x
Z
dx
= dx = tan x − sec x.
1 + sin x cos2 x
So uZ= −x + tan x − sec x.
cos x
v= dx. Let t = 1 + sin x: = ln |1 + sin x|.
1 + sin x
P.I. = (tan x−x−sec x) cos x+ln(1+sin x) sin x = sin x−1−x cos x+sin x ln(1+sin x).
Solution
A.E. : (m − 1)2 = 0 ⇒ m = 1, 1. y1 = ex , y2 = xex , X = ex sin x.
xex · ex sin x
Z Z
u=− 2x
dx = − x sin x dx = −(−x cos x + sin x) = x cos x − sin x.
Z x x e
e · e sin x
Z
v= dx = sin x dx = − cos x.
e2x
P.I. = (x cos x − sin x)ex + (− cos x) · xex = xex cos x − ex sin x − xex cos x = −ex sin x.
d2 y dy e3x
Example E5 — − 6 + 9y = 2 (TCET Resource, PYQ)
dx2 dx x
Solution
A.E. : (m − 3)2 = 0 ⇒ m = 3, 3. y1 = e3x , y2 = xe3x , X = e3x /x2 .
Solution
y1 = cos 2x, y2 = sin 2x, X = tan 2x.
W = 2 cos2 2x + 2 sin2 2x = 2.
sin 2x · tan 2x sin2 2x
Z Z Z
1 1
u = − dx = − dx = − (sec 2x − cos 2x) dx =
2 2 cos 2x 2
1 sin 2x
− ln | sec 2x + tan 2x| + .
4 4
cos 2x · tan 2x
Z Z
1 cos 2x
v= dx = sin 2x dx = − .
2 2 4
ln | sec 2x + tan 2x| sin 2x cos 2x cos 2x
P.I. = − + cos 2x+ − sin 2x = − ln | sec 2x+tan 2x|.
4 4 4 4
cos 2x
C.S. = c1 cos 2x + c2 sin 2x − ln | sec 2x + tan 2x|
4
Solution
2
y1 = cos
R x, y2 = 2sin x, W =R 1, X = sec x.
u = R− sin x sec x dx = R − tan x sec x dx = − sec x.
v = cos x sec2 x dx = sec x dx = ln | sec x + tan x|.
P.I. = − sec x cos x + sin x ln | sec x + tan x| = −1 + sin x ln | sec x + tan x|.
C.S. = c1 cos x + c2 sin x − 1 + sin x ln | sec x + tan x|
2
Example E8 — (D2 − 1)y = (TCET Resource)
1 + ex
Solution
y1 = exZ, y2 = e−x , X = 2/(1 + Zex ). W = ex (−e−x ) − e−x ex = −2.
e−x · 2/(1 + ex ) e−x
u=− dx = x
dx.
−2 Z1 + e Z
x −x 1/t dt dt
Let t = e , dx = dt/t, e = 1/t: · = 2
.
1+t t t (1 + t)
1 −1 1 1 1
Partial fractions: 2 = + 2 + . u = − ln t − + ln(1 + t) =
t (1 + t) t t 1+t t
−x x
−x − e
Z x + ln(1 + e ).
e · 2/(1 + ex ) ex
Z
v= dx = − x
dx = − ln(1 + ex ).
−2 1+e
P.I. = uex + ve−x = [−x − e−x + ln(1 + ex )]ex + [− ln(1 + ex )]e−x
= −xex − 1 + ex ln(1 + ex ) − e−x ln(1 + ex ) = −1 − xex + (ex − e−x ) ln(1 + ex ).
C.S. = c1 ex + c2 e−x − 1 − xex + (ex − e−x ) ln(1 + ex )
1
Example E9 — (D2 + 1)y =
1 + cos x
Solution
y1 = cos x, y2 = sin x, W = 1, X = 1/(1 + cos x).
2
Use 1 + cosZ x = 2 cos (x/2): Z Z
sin x 2 sin(x/2) cos(x/2)
u = − dx = − dx = − tan(x/2) dx =
2 cos2 (x/2) 2 cos2 (x/2)
2 ln | cos(x/2)|.
2 cos2 (x/2) − 1
Z Z Z
cos x 1 2
v = dx = dx = 1 − sec (x/2) dx = x −
2 cos2 (x/2) 2 cos2 (x/2) 2
tan(x/2).
1
Example E10 — (D2 + 9)y = (May 2024 Q.5c style)
1 + sin 3x
Solution
2 2
y1 = cosZ3x, y2 = sin 3x. W = 3Z cos 3x + 3 sin 3x = 3.
1 sin 3x 1 1
u=− dx = − 1− dx.
3 1 + sin
Z 3x 3 1 + sin 3x
1 − sin 3x
Z
dx tan 3x sec 3x
= 2
dx = − .
1 + sin 3x cos 3x 3 3
1 tan 3x sec 3x x tan 3x sec 3x
u=− x− + =− + − .
Z3 3 3 3 9 9
1 cos 3x 1
v= dx. Let t = 1 + sin 3x: = ln |1 + sin 3x|.
3 1 + sin 3x 9
x tan 3x sec 3x ln |1 + sin 3x|
P.I. = − + − cos 3x + sin 3x.
3 9 9 9
x cos 3x sin 3x 1 sin 3x
C.S. = c1 cos 3x + c2 sin 3x − + − + ln |1 + sin 3x|
3 9 9 9
Solution
y1 = exZ, y2 = e−x , W = −2, −x −x
Z X = e sin(e ) + cos(e ).
−x
x
e X 1
[ex · e−x sin(e−x ) + ex cos(e−x )] dx = − 12 [sin(e−x ) +
R
v = dx = −
−2 2
eRx cos(e−x )] dx.
ex cos(e−x )dx: let t = e−x , dt = −e−x dx = −dx/ex , so ex dx = −dt/t2 : = cos t ·
R
sin(e−x )e−x (−1) · ex − e−x sin(e−x )ex = yp′′ − yp = e−x sin(e−x ) + cos(e−x ) = X.
✓
C.S. = c1 ex + c2 e−x − ex sin(e−x )
Solution
y1 = cos
R x, y2 = sin x, WR = 1, X = csc x cot x = cos x/ sin2 x.
sin x cos x
u = R− sin2 x R − 2 cot x dxR= − 2ln | sin x|.
dx =
cos x·cos x
v= sin2 x
dx = cot x dx = (csc x − 1) dx = − cot x − x.
Solution
y1 = cos(ax), y2 = sin(ax), W = a. Z
ln | cos(ax)|
Z
1 1
u=− sin(ax) sec(ax) dx = − tan(ax) dx = .
a
Z a a2
1 x
v= cos(ax) sec(ax) dx = .
a a
cos(ax) x sin(ax)
P.I. = 2
ln | cos(ax)| + .
a a
cos(ax) x sin(ax)
C.S. = c1 cos(ax) + c2 sin(ax) + 2
ln | cos(ax)| +
a a
Solution
y1 = exZ, y2 = xex , W = e2x , XZ = ex log x.
xex · ex log x
u=− dx = − x log x dx.
e2x
R x2 x2 x2 x2
IBP: x log x dx = log x − . u = − log x + .
Z x x 2 Z 4 2 4
e · e log x
v= dx = log x dx = x log x − x.
e2x
2
x log x x2 x
2
3x2
x x x log x 2 2
P.I. = − + e + (x log x − x) · xe = e − + x log x − x
2 4 2 4
x2 log x x2 x
Wait, recompute: uy1 + vy2 = − + e + (x log x − x)xex =
2 2 4
x log x x2
2
x log x 3x2
x 2 2 x
e − + + x log x − x = e − .
2 4 2 4
x2 e x
C.S. = (c1 + c2 x)ex + (2 log x − 3)
4
Solution
Put x = ez : xy ′ = Dy, x2 y ′′ = D(D − 1)y. Equation becomes:
Example F2 — x2 y ′′ − 2xy ′ + 2y = x − 1
x
Solution
Put x = ez : (D2 − 3D + 2)y = ez − e−z . A.E. : (D − 1)(D − 2) = 0 ⇒ D = 1, 2.
C.F. = c1 ez + c2 e2z .
ez
P.I.1 = 2 : f (1) = 1 − 3 + 2 = 0, f ′ (1) = −1. P.I.1 = −zez .
D − 3D + 2
e−z e−z
P.I.2 = 2 : f (−1) = 1 + 3 + 2 = 6. P.I.2 = .
D − 3D + 2 6
C.S. in z: c1 ez + c2 e2z − zez + e−z /6.
Replace z = ln x:
1
y = c1 x + c2 x2 − x ln x +
6x
ln x
Example F3 — x2 y ′′ + 5xy ′ + 3y = (Pre-NEP Dec 2024 Q.6a)
x2
Solution
Put x = ez : (D2 + 4D + 3)y = ze−2z . A.E. : (D + 1)(D + 3) = 0 ⇒ D = −1, −3.
C.F. = c1 e−z + c2 e−3z .
ze−2z −2z z
P.I. = . Shift (D → D − 2): = e =
D2 + 4D + 3 (D − 2)2 + 4(D − 2) + 3
z z
e−2z 2 = −e−2z .
D −1 1 − D2
(1 − D2 )−1 ≈ 1 + D2 : z + D2 (z) = z + 0 = z.
P.I. = −ze−2z .
C.S. in z: c1 e−z + c2 e−3z − ze−2z .
Replace z = ln x, e−z = x−1 , e−2z = x−2 :
c1 c2 ln x
y= + 3− 2
x x x
Solution
z
Put x = ez : (D2 −1)y = e2z ·ee . A.E. : D2 −1 = 0 ⇒ D = ±1. C.F. = c1 ez +c2 e−z .
z
Use VPM:
Z −z y1 = ezz , y2 = e−z , X = e2z ee , zW = −2.
e e2z ee ee
Z
1 z
u=− dz = ez ee dz = .
Z z 2z ez−2 2
Z 2
ee e 1 z
v= dz = − e3z ee dz.
−2 2
ez
Let t = e : = h− 2 t e dt = − 21 et (ti2 − 2t + 2) = − e2 (e2z − 2ez + 2).
1
z
R 2 t
ez ez ez ee
z ez z z
P.I. = e2 ·ez + − e2 (e2z − 2ez + 2) e−z = 2
− e2 (ez −2+2e−z ) = ee (1−e−z )ee . . .
Replace ez = x:
c2
y = c1 x + + xex (1 − 2/x + · · · ) (Cauchy-VPM combination)
x
(This level of combination is rarely set; knowing the substitution step earns full marks.)
Solution
Put 1 + x = ez (z = ln(1 + x)): (D2 + 1)y = 4 cos z.
A.E. : D2 + 1 = 0 ⇒ D = ±i. C.F. = c1 cos z + c2 sin z.
4 cos z x sin z
P.I. = 2 : D2 → −1, denom= 0 (failure). P.I. = 4 · = 2z sin z.
D +1 2
C.S. in z: c1 cos z + c2 sin z + 2z sin z.
Replace z = ln(1 + x):
Solution
Put 3x + 2 = ez : (3x + 2)y ′ = 3Dy, (3x + 2)2 y ′′ = 9D(D − 1)y.
2 z
[9D(D − 1) + 9D − 36]y = e 3−2 + 4(e 3−2) + 1 = e9 + 8e9 − 11
z 2z z
9
.
2 e2z 8ez 11 2 e2z 8ez 11
(9D − 36)y = 9 + 9 − 9 , i.e., (D − 4)y = 81 + 81 − 81 .
A.E. : D = ±2. C.F. = c1 e2z + c2 e−2z .
ze2z ze2z
P.I.1 : f (2) = 0, f ′ (2) = 4. P.I.1 = = .
z
81 · 4 z 324
8e 8e
P.I.2 : f (1) = −3. P.I.2 = =− .
81(−3) 243
−11/81 11
P.I.3 : f (0) = −4. P.I.3 = = .
−4 324
Replace z = ln(3x + 2), ez = 3x + 2:
Solution
Put 2x + 1 = ez , so x = (ez − 1)/2. (2x + 1)y ′ = 2Dy, (2x + 1)2 y ′′ = 4D(D − 1)y.
[4D(D − 1) − 4D − 12]y = 3(ez − 1) (4D2 − 8D − 12)y = 3ez − 3, i.e., (D2 − 2D − 3)y =
3ez − 3
.
4
A.E. : (D − 3)(D + 1) = 0 ⇒ D = 3, −1. C.F. = c1 e3z + c2 e−z .
3ez /4 3ez 3ez
P.I.1 = 2 : f (1) = 1 − 2 − 3 = −4. P.I.1 = = − .
D − 2D − 3 D=1 4(−4) 16
−3/4 −3/4 1
P.I.2 = 2 D=0
: f (0) = −3. P.I.2 = = .
D − 2D − 3 z
−3 4
3e 1
C.S. in z: c1 e3z + c2 e−z − + .
16 4
c2 3(2x + 1) 1
y = c1 (2x + 1)3 + − +
2x + 1 16 4
Solution
CF: (m − 1)2 = 0 ⇒ m = 1, 1. C.F. = (c1 + c2 x)ex .
VPM: From Example E4: P.I. = −ex sin x.
CS: y = (c1 + c2 x)ex − ex sin x .
Verification: Let yp = −ex sin x. yp′ = −ex sin x − ex cos x. yp′′ = −ex sin x −
ex cos x − ex cos x + ex sin x = −2ex cos x. yp′′ − 2yp′ + yp = −2ex cos x − 2(−ex sin x −
ex cos x) + (−ex sin x) = −2ex cos x + 2ex sin x + 2ex cos x − ex sin x = ex sin x ✓.
Frequency Heatmap
x
■■■ Very High (every paper): VPM: (D2 + 3D + 2)y = ee or sin(ex ). VPM:
y ′′ + y = sec x tan x. VPM: (D2 − 2D + 1)y = ex sin x.
■■ High (3–4 papers): CF+PI: (D2 − 4D + 4)y = e2x + trig. PI Type 5: eax V
combinations. Cauchy-Euler: x2 y ′′ + . . .
■ Moderate (2 papers): PI Type 4: polynomial RHS. Legendre’s: (ax + b)2 y ′′ +
. . . y ′′ + 4y = tan 2x (VPM).
CF: Complex eαx (c1 cos βx + c2 sin βx) Pure imaginary: drop e0x = 1
α ± iβ
End of Chapter 2
Next: Chapter 3 — Beta & Gamma Functions + DUIS