2024 - Lecture Notes - Differential Equations
2024 - Lecture Notes - Differential Equations
Maritime Transport
Department of Basic & Applied Sciences
Di¤erential Equations
Lecture & Exercise Notes
2024
Contents
1
CONTENTS 2
B Assignments 140
B.1 Assignment 1 on Separable and
Homogeneous Di¤erential Equations . . . . . . . . . . . . . . . 141
CONTENTS 4
a0 = an = bn =
1
R 1
R 1
R
general case f (x) dx f (x) cos nx dx f (x) sin nx dx
2
R
f (x) is odd 0 0 0
f (x) sin nx dx
2
R 2
R
f (x) is even 0
f (x) dx 0
f (x) cos nx dx 0
1
CHAPTER 1. THE FOURIER SERIES 2
) Ff (x) = 2
4
cos x + 19 cos 3x + 1
25
cos 5x + 1
49
cos 7x + .
1 x < 0;
Problem 1.2.2 f (x) =
1 0 x< :
This function
R is odd. ) a0 = R0; an = 0 for all n.
bn = 2 0 f (x) sin nx dx = 2 0 sin nx dx = n 2 cos nx 0
= n
2
cos n + n
2
= n
2
( 1)n + n
2
) bn = n
2
( ( 1)n + 1) :
) b1 = ; 4
b2 = 0; b3 = 34 ; b4 = 0; b5 = 54 ;
b6 = 0; b7 = 74 ; b8 = 0; b9 = 94 :
P1 X1
a0
Ff (x) = 2
+ (an cos nx + bn sin nx) = 4 1
2j+1
sin ((2j + 1) x) :
n=1 j=0
) Ff (x) = 4
sin x + 13 sin 3x + 15 sin 5x + 71 sin 7x + 91 sin 9x + .
P
1 8 8 8 8 8 8 8
Hence, = = + + + + + + :
n=0 (4n + 1) (4n + 3) 3 35 99 195 323 483
P
n=100 8
For instance, = 3:1366
n=0 (4n + 1) (4n + 3)
) Ff (x) = 2 sin x 2
2
sin 2x + 23 sin 3x 2
4
sin 4x + 25 sin 5x .
) Ff (x) =
2 4 4 4 4 4
3 1
cos x + cos 2x
4 9
cos 3x + 16
cos 4x 25
cos 5x + .
First-order Di¤erential
Equations
We learn how to recognize equations of each one of these …ve types, and
how to obtain their analytic solutions.
dy
x = 3y:
dx
7
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 8
C1 x3 ; x 0
y (x) = ;
C2 x3 ; x<0
y (x) = C1 x3 :
dy A (x)
= , (2.2)
dx B (y (x))
dy
which is rewritten as B (y (x)) dx = A (x). To solve this separable D.E., we
integrate both sides, and get
Z Z
dy
B (y (x)) dx = A (x) dx + C:
dx
dy y (x)
=f .
dx x
u + x du
dx
= f (u); that is,
du
x = f (u) u:
dx
dy 2 R ydy
R dx
y (4 + x2 ) dx = (y 2 + 25) Separable D.E. (y 2 +25)2
= 4+x2
R R dx
) 2y (y 2 + 25) dy = 2 4+x
2
2
) ) y 2 + 25 =
1
(y 2 + 25) = 22 tan 1 x
C 1
.
2 C tan 1
( x2 )
dy
Problem 2.1.2 (4xy 3 + 2x3 y) dx = 3y 4 + x2 y 2 + 5x4 :
y 3 y 4 y 2
Divide by x4 : ) 4 x
+2 y
x
dy
dx
=3 x
+ x
+5
Homogeneous D.E., Put u = y
x
.) y = xu ) dy
dx
= u + x du
dx
) (4u3 + 2u) u + x du
dx
= 3u4 + u2 + 5
) 4u4 + 2u2 + (4u3 + 2u) x du
dx
= 3u4 + u2 + 5
) (4u3 + 2u) x dx = (u4 + u2 5)
du
Separable D.E. in u (x).
R 4u3 +2u R dx
) u4 +u2 5 du = x
) ln (u + u2 5) = ln (C1 x) = ln
4 C
x
4 2
) u4 + u2 5 = Cx ) xy 4 + xy 2 5 = Cx
) y 4 + x2 y 2 5x4 = Cx3 .
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 10
dy
Problem 2.1.4 dx
5x4 = 5x4 y 2 ; y (0) = 1:
R dy R 4
dy
dx
= 5x4 (1 + y 2 ) Separable D.E. ) 1+y 2 = 5x dx
) tan y = C + x ; i.e. y = tan (C + x ) :
1 5 5
At x = 0; y = 1: ) tan 1 (1) = C + 0 ) 4 =C
) y = tan 4
+ x5 .
2y dy 2y 2y
Problem 2.1.5 2x3 cos x
+ 15xy 2 dx
= 2x2 y cos x
+x3 sin x
+20y 3 :
Divide by x3 :
y 2 y 3
) 2 cos 2yx
+ 15 x
dy
dx
= 2 xy cos 2y
x
+ sin 2y
x
+ 20 x
dy
Problem 2.1.6 (x7 8x4 y 3 + 7xy 6 ) dx = 4y 7 2x3 y 4 2x6 y;
y (1) = 1:
Divide by x7 :
3 6 dy y 7 y 4
1 8 xy + 7 xy dx
=4 x
2 x
2 y
x
:
Homogeneous D.E.
Put u = xy . ) y = xu ) dy
dx
= u + x du
dx
) (1 8u3 + 7u6 ) u + x du
dx
= 4u7 2u4 2u
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 11
) (1 8u3 + 7u6 ) x du
dx
+u 8u4 + 7u7 = 4u7 2u4 2u
) (1 8u3 + 7u6 ) x dudx
= 3u7 + 6u4 3u: Separable D.E. in u (x).
R R
) 1u7 8u2u+7u
3 6
4 +u du = 3 dx x
) y7 2x3 y 4 + x6 y = Kx4 :
At x = 1, y = 1, )K= 1 2 1= 4:
) y7 2x3 y 4 + x6 y = 4x4 :
p q
dy y+ x2 y 2 y y 2
Problem 2.1.7 dx
= x
= x
+ 1 x
; y (1) = 0
) p1duu2 = dx x
) sin 1
u = ln Cx
) x = u = sin (ln Cx) :
y
At x = 1, y = 0. ) 0 = sin (ln C)
) 0 = ln C )C=1 ) y = x sin (ln x) :
dy y
Problem 2.1.8 x dx = y + x cot x
:
) dy
dx
= y
x
+ cot y
x
Homogeneous D.E. Put u = y
x
.
) y = xu ) =u+ dy
dx
)u+ x du
= u + cot u
dx
x du
dx
)x du
Separable D.E. in u (x).
R dx = cot u R dx
) tan u du = x ) ln (sec u) = ln (Cx)
) sec u = Cx ) y
x
= u = sec 1
(Cx) ) y = x sec 1
(Cx) .
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 12
dy 2y 2 +7xy+5x2
Problem 2.1.10 dx
= xy+x2
; y (1) = 3:
y 2
2( x ) +7 xy +5
) dy
dx
= y
+1
: Homogeneous D.E.
x
Put u = xy . ) y = xu ) dx
dy
= u + x du
dx
) u + x dx = u+1
du 2u2 +7u+5
) x du
2 +7u+5 2
= 2u u+1 u = u +6u+5 = (u+1)(u+5) =u +5
R dxdu R dx u+1 u+1
) u+5 = x ) ln (u + 5) = ln (Cx) :
) x + 5 = u + 5 = Cx:
y
that is, y = Cx2 5x .
At x = 1; y = 3: ) C = 8: 2
) y = 8x 5x .
dy
Problem 2.1.11 (5xy 4 3x3 y 2 + 2x5 ) dx = 2y 5 4x4 y:
y 4 y 2 dy y 5
Divide by x5 : 5 x
3 x
+2 dx
=2 x
4 xy :
) (5u4 3u2 + 2) u + x du
dx
= 2u5 4u
) 5u 5 3
3u + 2u + (5u 4
3u + 2) x du
2
dx
= 2u5 4u
) (5u4 3u2 + 2) x du
dx
= 3u5 + 3u3 6u = 3 (u5 u3 + 2u)
(5u4 3u2 +2)du
Separable D.E. in u (x). ) u5 u3 +2u = 3dx
x
R 4 3u2 +2 R dx
) 5u5 3
u u +2u
du = 3 x
= 3 ln Cx
) ln (u5 u3 + 2u) = ln C1
x3
y5 y3
) u5 u3 + 2u = C1
x3
: ) x5
+ 2 xy = Cx31 :
x3
Multiply by x5 : ) y5 x y + 2x4 y = C1 x2 :
2 3
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 13
dy
Problem 2.1.12 (x2 2x + 1) dx 2y = 0:
p dy
Problem 2.1.16 2 1 x2 dx sec y sin3 y = 0:
dy
Problem 2.1.17 (3xy 2 + 10x2 y) dx + 6x3 5y 3 20xy 2 = 0:
y dy ex
Problem 2.1.18 x dx
= ln y
:
y
Problem 2.1.19 xy 0 = y + x cot x
:
1
Final Answer: Homogeneous D.E. y = x sec (Cx) :
dy
Problem 2.1.20 x2 + 3y 2 2xy dx = 0:
p
Final Answer: Homogeneous D.E. y= x Cx 1:
p dy
p
Problem 2.1.21 xy x2 + y 2 dx = x3 + y 2 x2 + y 2 :
2=3
Final Answer: Homogeneous D.E. y 2 = x2 [ln (Cx3 )] x2 :
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 14
d
) = p: (2.7)
dx
So, given any di¤erentiable function y (x), we …nd that
dy
d
dx
( y) = dx + ddx y = dx
dy
+ py = dy
dx
+ py :
Therefore, a function y (x) becomes a solution of the di¤erential equation
dy
+ py = q
dx
d dy
if and only if dx ( y) = dx
+ py = q:
When we integrate both sides w.r.t. x, we obtain the following formula
for all solutionsR of equation (2.4):
y=C+ q dx;
as claimed in (2.6).
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 15
dy
Linear D.E., dx
+ p (x) y = q (x) :
5x
R R
p = tan x; q = 10e cos x; p dx = tan x dx = ln (sec x) :
R
p dx
Integrating factor: =e = eln(sec x) = sec x; q = 10e 5x
:
R R
)y= 1 C+ q dx = cos x C + 10e 5x dx = cos x (C 2e 5x
):
At x = 0, y = 0, )0=C 2 ) C = 2:
5x
) y = 2 cos x (1 e ):
dy
Problem 2.2.3 dx
+ (1 + tan x) y = cos x:
dy
Linear D.E., dx
+ p (x) y = q (x) :
p = 1 + tan x q = cos x:
R R
p dx = (1 + tan x) dx = x + ln (sec x) :
R
p dx
Integrating factor: =e = ex+ln(sec x) = ex eln(sec x) = ex sec x:
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 16
R R R
q dx = ex sec x cos x dx = ex dx = ex :
R
)y= 1
C+ 1
q dx = ex sec x
(C + ex )
) y = Ce x
cos x + cos x :
dy
Problem 2.2.4 dx
+ (2 3 cot 3x) y = 8 sin 3x:
dy
Linear D.E., dx
+ p (x) y = q (x) :
p = 2 3 cot 3x q = 8 sin 3x:
R R
p dx = (2 3 cot 3x) dx = 2x ln (sin 3x) :
R
e2x e2x
Integrating factor: = e p dx = e2x ln(sin 3x) = eln(sin 3x)
= sin 3x
:
R R e2x R 2x
q dx = 8 sin 3x
sin 3x dx = 8 e dx = 4e2x :
R
)y= 1 C+ q dx = e 2x sin 3x (C + 4e2x )
) y = Ce 2x
sin 3x + 4 sin 3x :
dy
Problem 2.2.6 x2 dx + 3xy = 3ex :
dy 3ex 3ex
Linear D.E. dx
+ x3 y = x2
:, p= 3
x
q= x2
R R 3 3
p dx dx
Integrating factor: =e =e x = e3 ln x = eln x = x3 :
C+
R
q dx C 3
R
x ex dx 1
R x
y= = x3
= x3
C 3 x ex e dx
)y= 1
x3
(C 3 (x ex ex )) : )y= 1
x3
(C 3x ex + 3ex ) :
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 17
dy 1 1
Problem 2.2.7 dx
+ x
cot x y = x
cos x csc x:
1 1
Linear D.E. p= x
cot x q= x
cos x csc x
Integrating
R factor:
R 1
=e p dx
= e x cot x dx = eln x ln(sin x) = eln(x= sin x) = sinx x :
R R R R
q dx = cot x x csc2 x dx = cot x dx+x cot x cot x dx = x cot x
R
)y= 1 C+ q dx = sinx x (C + x cot x)
) y = C sinx x + cos x:
dy 2
Problem 2.2.8 dx
+ 1 x
y = x2 :
2
Linear D.E. p=1 x
q = x2
Integrating
R
factor:
R
= e 1 x dx = ex 2 ln x = ex eln(x ) =
2 2
p dx ex
=e x2
:
R R x R
q dx = xe 2 :x2 dx = ex dx = ex
R
)y= 1 C+ q dx = x2 e x (C + ex )
) y = Cx2 e x
+ x2 :
dy
Problem 2.2.9 dx
+ 3 (tan x) y = ex cos3 x; y (0) = 5:
dy
Linear D.E., dx
+ p (x) y = q (x) :
p = 3 tan x q = ex cos3 x:
R R
p dx 3 tan x dx
Integrating factor: =e =e = e3 ln(sec x) = sec3 x:
R R
)y= 1 C+ q dx = sec13 x C + sec3 x ex cos3 xdx
R x 3
) y = cos3 x C + ex dx ) y = (C + e ) cos x:
At x = 0, y = 5, ) 5 = 1 (C + 1) ) C = 4:
x 3
)y = (4 + e ) cos x.
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 18
dy 2x 2
Problem 2.2.10 dx
+ x2 16
y = x 4
:
dy
Linear D.E., dx
+ p (x) y = q (x) :
p = x22x16 q = x2 4 :
R 2x
R dx
= eln(x 16)
p dx x2 16 2
Integrating factor: =e =e = x2 16:
R R 2(x2 16) R
y= 1 C+ q dx = x2 1 16 C+ x 4
dx = 1
x2 16
C+ 2x + 8dx
) y= x2 +8x+C
x2 16
:
)y= 1
C+ q dx = cos1 x (C + x3 cos x) = C sec x + x3 :
dy
Problem 2.2.12 dx
+ (1 + tan x) y = sin x + sec x:
dy
Linear D.E., dx
+ p (x) y = q (x) :
p = 1 + tan x q = sin x + sec x:
R R
p dx = (1 + tan x) dx = x + ln (sec x) :
R
Integrating factor: = e p dx = ex+ln(sec x) = ex eln(sec x) = ex sec x:
R R x R R
R x q dx = e sec x (sin xR + sec x) dx = ex tan x dx + ex sec2 x dx
= e tan x dx + ex tan x ex tan x dx = ex tan x:
R
)y= 1 C+ q dx = ex sec 1
x
(C + ex tan x)
) y = Ce x
cos x + sin x :
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 19
dy
Problem 2.2.13 dx
+ (1 cot x) y = 0:
dy
Linear D.E., dx + p (x) y =R q (x) ; R p = 1 cot x q = 0:
p dx 1 cot x dx x ln(sin x) ex
Integrating factor: =e =e =e = sin x
R
) y= 1
C+ x
q dx = C e sin x:
dy
Problem 2.2.14 dx
+ (1 cot x) y = 2e x
sin2 x
dy
Linear D.E., dx + p (x) y =R q (x) ; p = 1 cot x q = 2e x sin2 x
ex
RIntegrating factor:
R = e p dx = sin x
(see Problem 2.2.13)
q dx = 2 sin x dx = 2 cos x: Final Answer:
R
y= 1 C+ q dx = C e x sin x e x sin 2x = e x (C sin x sin 2x) :
dy
Problem 2.2.15 dx
+ (1 cot x) y = 2 sin x cos x 2 sin2 x
dy
Linear D.E. dx +p (x) y = qR(x) ; p = 1 cot x q = 2 sin x cos x 2 sin2 x
ex
Integrating
R factor:
R x = e p dx = sin x
(see Problem 2.2.13)
x
q dx = 2 e (cos x sin x) dx = 2e cos x
R
Final Answer y = 1 C + q dx = C e x sin x + sin 2x:
dy
Problem 2.2.16 dx
+ (1 cot x) y = sin x
dy
Linear D.E., dx + p (x) y =R q (x) ; p = 1 cot x q = sin x
p dx ex
RIntegratingRfactor:
x x
=e = sin x (see Problem 2.2.13)
q dx = e dx = e
R
Final Answer y = 1 C + q dx = (C e x + 1) sin x:
dy
Problem 2.2.17 dx
+ (1 cot x) y = cos x csc x
dy
Linear D.E. dx + p (x) y = Rq (x) ; p = 1 cot x; q = cos x csc x
ex
Integrating
R R factor: = e p dx =Rsin x
(see Problem 2.2.13)
R x
q dx = e cot x e csc x dx = e cot x dx+ex cot x
x x 2 x
e cot x dx =
ex cot x R
Final Answer y = 1 C + q dx = C e x sin x + cos x:
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 20
2y 2 +7xy+5x2
Problem 2.2.18 y 0 = xy+x2
; y (1) = 3: (See Problem 2.1.10.)
2y 2 +7xy+5x2 (2y+5x)(y+x)
y0 = xy+x2
= x(y+x)
:
) y0 = 2y+5x
x
= x2 y + 5 ) y0 2
x
y = 5:
Final Answer: Linear D.E., y = Cx2 5x:
At x = 1; y = 3: ) C = 8: ) y = 8x2 5x .
dy 2
Problem 2.2.19 dx x
y = 4x2 sec2 x:
dy
p
Problem 2.2.20 dx
+ (tan x) y = sec x; y (0) = 3:
dy
Linear D.E., dx
+ p (x) y = q (x) ; p = tan x q = sec x:
R R
Integrating factor: = e p dx = e tan x dx = eln(sec x) = sec x:
R R
y= 1 C+ q dx = cos x C + sec2 x dx = cos x (C + tan x)
) yG:S: = C cos x + sin x:
p p
When x = 0; y = 3: ) 3 = C cos 0 + sin 0 = C:
p
) y = 3 cos x + sin x = 2 sin x + 3 :
y 2
0
-5 0 5 10
-1
-2
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 21
and get
n dy
(1 n) y + (1 n) g (x) y 1 n
= (1 n) h (x) :
dx
The substitution
z = y1 n
(2.10)
dz dy
(and so, dx = (1 n) y n dx ) turns this Bernoulli D.E. into the following
linear di¤erential equation in the unknown function z (x) :
dz
+ (1 n) g (x) z = (1 n) h (x) :
dx
We obtain its solution z (x). Then we get y (x) from (2.10).
dy
Problem 2.3.3 dx
(cos x) y = x3 e4 sin x y 3 : Bernoulli D.E., n = 3:
n
Multiply by (1 n) y = 4y 3 .
) 4y 3 dx
dy
4 (cos x) y 4 = 4x3 e4 sin x : Put z = y 4 . ) dz
dx
dy
= 4y 3 dx :
dz
) D.E. becomes Linear in z (x) : dx
4 (cos x) z = 4x3 e4 sin x :
p p
) y = 4 e4 sin x (C + x4 ) = esin x 4 C + x4 :
dy p
Problem 2.3.4 dx
+ 2y = 6e2x y:
n 1
Bernoulli D.E, n = 1=2 Multiply by (1 n) y = p
2 y
.
p p
) 2p1 y dx
dy
+ y = 3e2x Put z = y: ) dx
dz
= 1 dy
p
2 y dx
:
) D.E. becomes dz
dx
+ z = 3e2x :
Linear D.E. in z, p = 1, q = 3e2x :
R
p dx
Integrating factor: =e = ex :
p R R
y=z= 1 C+ q dx = e x
C+ 3e3x dx :
p
) ) y = (Ce
x 2
y=e (C + e3x ) = Ce x
+ e2x x
+ e2x ) :
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 23
dy
Problem 2.3.6 4 cos x dx + (sin x) y = y 3 :
(1 n)y n =4y 3
Bernoulli D.E., n = 3. Multiply by 4 cos x
:
) 4y 3 dx
dy
+ (tan x) y 4 = sec x: Put z=y : ) dx
dz 4 dy
= 4y 3 dx :
) D.E. becomes dz
dx
+ (tan x) z = sec x
R LinearR D.E. in z, p = tan x, q = sec x: Integrating factor: =
e p dx = e tan x Rdx = eln(sec x) = sec x: R
z= 1 C+ q dx = cos x C + sec2 x dx = cos x (C + tan x)
p
y 4 = z = C cos x + sin x: ) y = 4 C cos x + sin x:
dy
Problem 2.3.7 dx
+ 14 (tan x) y = 2 (e x
cos x) y 3 :
n
Bernoulli D.E., n = 3: Multiply by (1 n) y = 4y 3 .
) 4y 3 dx
dy
+ (tan x) y 4 = 8e x
cos x: Put z = y 4 . ) dz
dx
dy
= 4y 3 dx :
dz x
) D.E. becomes Linear in z : dx
+ (tan x) z = 8e cos x: p = tan x
q = 8e x cos x::
R R
p dx tan x dx
Integrating factor: =e =e = eln sec x = sec x:
R R
) y4 = z = 1
C+ q dx = cos x C + 8e x
dx = cos x (C 8e x ) :
p
) ) y= 4
cos x (C 8e x ):
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 24
dy
Problem 2.3.8 x2 dx xy = y 2 :
2
) dx
dy y
x
= xy Bernoulli D.E., n = 2
Multiply by (1 n) y n = y 2 . ) y 2 dxdy
+ x1 y 1
= 1
x2
Put z = y 1: ) dx
dz
= y 2 dxdy
: ) D.E. becomes dz
dx
+ x1 z =
1
x2
: R
p dx
Linear D.E. in zR (x), Integrating factor:
R = e = eln x = x:
1 1 1 1 1
y
=z= C+ q dx = x C x
dx = x (C ln x)
) y= x
C ln x
= x
ln K+ln x
= x
ln(Kx)
:
dy 1
Problem 2.3.9 (sec x) dx 3
(sec2 x) y = y4; y (0) = 1:
) dy
dx
1
3
(sec x) y = (cos x) y 4 :
n
Bernoulli D.E, n = 4 Multiply by (1 n) y = 3y 4 .
) 3y 4 dy
dx
+ (sec x) y 3
= 3 cos x:
Put z = y 3
. ) dz
dx
= 3y 4 dy
dx
:
) D.E. becomes dz
dx
+ (sec x) z = 3 cos x:
R R
Linear D.E. in z, p dx = sec x dx = ln (sec x + tan x) :
R
p dx
Integrating factor: =e
= sec x + tan x:
R R
y 3=z= 1 C+ q dx = sec x+tan1
x
C + 3 cos x (sec x + tan x) dx
R
) y 3 = sec x+tan
1
x
C + 3 + 3 sin x dx = C+3x 3 cos x
sec x+tan x
:
q
) y = 3 C+3x sec x+tan x
3 cos x
:
At x = 0, y = 1. )1= C+0 3
1+0
) C = 4.
q
) y = 3 4+3x
sec x+tan x
3 cos x
:
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 25
dy
Problem 2.3.10 dx
+ x3 y = 2x9 y 3 ; y (1) = 1:
dy
Problem 2.3.11 (2x3 9xy 2 ) dx = 2x3 + 6x2 y 15y 3 ; y (1) = 1:
dy
Problem 2.3.12 (6xy 2 3x2 y) dx = 10y 3 6xy 2 + 2x3 ; y (1) = 2:
dy
Problem 2.3.13 dx
(cot x) y = 2 (cos x) y 2 :
dy
Problem 2.3.14 x dx + 3y = (x7 tan x) y 3 :
dy 4
Problem 2.3.15 5 cos x dx + (sin x) y = y :
dy 1
Problem 2.3.16 dx 3
(sec x) y = (cos x) y 4 ; y (0) = 1:
dy
Problem 2.3.17 dx
+ 34 y = (e 3x
cos 2x) y 3 :
dy
Problem 2.3.18 dx
+ 14 (tan x) y = 3 (e3x cos x) y 3 :
dy
Problem 2.3.19 dx
+ x2 y = 2x4 y 3 :
dy
Problem 2.3.20 dx
+ x1 y = 1 2x
x3
e y 2:
dy
Problem 2.3.21 3 dx (tan x) y = 2 (sin x) y 2 :
dy 1
Problem 2.3.22 dx x
y = x2 y 2 :
dy
Problem 2.3.23 x2 dx xy = y 2 :
dy
Problem 2.3.24 dx
+ 2 (cot x) y = 8 (cos x) y 1=2 :
Final Answer: Bernoulli D.E., n = 1=2:
2
y = csc2 x (C 2 cos2 x) :
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 26
@ (x; y) @ (x; y)
= M (x; y) and = N (x; y) : (2.14)
@x @y
That function (x; y) is one and only one (up to additive constants).
Proof. First, suppose that M and N satisfy the exactness condition (2.13).
Choose a point (x ; y ) in their common region of de…nition. De…ne a function
(x; y) as follows:
Z x Z y
(x; y) = M (t; y) dt + N (x ; u) du: (2.15)
x y
@y
@
= @y @x
= @ @y@x = @ @x@y @
= @x @ (x;y)
@y
= @N@x
(x;y)
;
that is, M and N must satisfy the exactness condition (2.13). This com-
pletes the proof.
Our task in Step 2 is to determine that function (x; y). We can do that
through computing the following two integrals, then comparing them with
each other:
Z
(x; y) = M (x; y) dx + A (y) ;
Z
(x; y) = N (x; y) dy + B (x) :
(Alternatively, you can compute from formula (2.15). But this could be a
bit more laborious!)
Step 3: We recall the following basic theorem of Di¤erential Calculus:.
Theorem 7 (Chain Rule in Several Variables) Let z (x), y (x) and (z; y)
be continuously di¤erentiable functions around points x and (z (x) ; y (x)), re-
spectively. Then at that point x we have:
dz
In the special case that z (x) = x; we have dx
= 1. Therefore, the Chain
Rule becomes:
d (x; y (x)) @ (x; y) @ (x; y) dy
= + ;
dx @x jat (x;y(x)) @y jat (x;y(x))
dx
d (x; y (x))
= 0:
dx
But, the only functions of one variable whose derivatives are identically 0 are
the constant functions. We conclude that the general solution y (x) of the
di¤erential equation (2.11) is implicit in the equation:
M (x; y) dx + N (x; y) dy = 0
@M (x;y) @ (9 cos x+y 2 cos x)
@y
= @y
= 2y cos x;
@N (x;y)
= @2y@xsin x = 2y cos x. ) @M@y(x;y)
= @N@x
(x;y)
: ) Exact D.E.
@x R R
(x; y) = M (x; y) dx+A (y) = (9 + y ) cos x dx+A (y) = (9 + y 2 ) sin x+
2
A (y) ; R R
(x; y) = N (x; y) dy + B (x) = 2y sin x dy + B (x) = y 2 sin x + B (x) :
p
(x; y (x)) = (9 + y 2 ) sin x = C: So by initial condition y = 25 csc x 9 .
2
Another
R cos Solution:
R 2ySeparable D.E. (9 + y ) cos x dx = 2y sin x dy
) x
sin x
dx = 9+y2 dy
) ln (C sin x) = ln (9 + y 2 ) ) 9 + y 2 = C sin
1
x
:
p
At x = 2 , y = 4. ) C = 1
25
and y = 25 csc x 9 .
) (x; y) = x4
4
ln y sin x sin4 y + sec y = C :
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 30
dy
Problem 2.4.4 (x3 sec2 y + x sec y tan y 2e2y ) dx = (3x2 tan y + cos x + sec y) :
M (x; y) dx + N (x; y) dy = 0
@M (x;y)
@y
= @(cos x@yy sin x) = sin x; @N (x;y)
@x
= @(2y+cos x)
@x
= sin x:
) @M (x;y)
@y
= @N (x;y)
@x
:) Exact D.E.
R R
(x; y) = M (x; y) dx + A (y) = (cos x y sin x) dx + A (y)
= sin x + y cos x + A (y) ;
R R
(x; y) = N (x; y) dy + B (x) = (2y + cos x) dy + B (x)
= y 2 + y cos x + B (x) :
) (x; y (x)) = sin x + y cos x + y 2 = C:
At x = 0; y = 2: ) sin 0 + 2 cos 0 + 4 = C. ) C = 6
) sin x + y cos x + y = 6 .
2
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 31
M (x; y) dx + N (x; y) dy = 0
@M (x;y) @ (2e2x sin y+4x3 tan 1 y ) 4x3
@y
= @y
= 2e2x cos y + 1+y 2
;
4
@ e2x cos y+ x 2 +sec2 y
@N (x;y) 4x3
@x
= @x
1+y
= 2e2x cos y + 1+y 2
:
) @M@y (x;y)
= @N@x(x;y)
: ) Exact D.E.
R R
(x; y) = M (x; y) dx + A (y) = (2e sin y + 4x3 tan 1 y) dx + A (y)
2x
3
Problem 2.4.7 (2xey + cos x) dx + x2 ey y
dy = 0:
M (x; y) dx + N (x; y) dy = 0
3
@M (x;y) y @ (x2 ey )
@y
= @(2xe @y+cos x) = 2xey ; @N (x;y)
@x
= @x
y
= 2xey :
) @M@y(x;y)
= @N@x
(x;y)
: ) Exact D.E.
R R
(x; y) = M (x; y) dx + A (y) = (2xey + cos x) dx + A (y)
= x2 ey + sin x + A (y) ;
R R
(x; y) = N (x; y) dy + B (x) = x2 ey y3 dy + B (x)
= x2 ey 3 ln y + B (x) :
) (x; y (x)) = x2 ey + sin x 3 ln y = C:
dy
Problem 2.4.8 (x5 sec2 y + 5y 4 ) dx = 5x4 tan y sec2 x:
) (5x4 tan y + sec2 x) dx + (x5 sec2 y + 5y 4 ) dy = 0:
@ (5x4 tan y+sec2 x)
M (x; y) dx + N (x; y) dy = 0: @M@y (x;y)
= @y
= 5x4 sec2 y;
@ (x5 sec2 y+5y 4 )
@N (x;y)
= = 5x4 sec2 y: ) @M@y (x;y)
= @N@x(x;y)
: ) Exact D.E.
@x R @x R 4 2
(x; y) = M (x; y) dx + A (y) = (5x tan y + sec x) dx + A (y)
= x5 tan
R y + tan x + A (y) ; R
(x; y) = N (x; y) dy + B (x) = (x5 sec2 y + 5y 4 ) dy + B (x)
= x5 tan y + y 5 + B (x) :
) (x; y (x)) = x5 tan y + tan x + y 5 = C:
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 32
1
Problem 2.4.9 (sec2 x + yexy ) dx + y
+ xexy dy = 0; y (0) = 1
M (x; y) dx + N (x; y) dy = 0
@M (x;y) @ (sec2 x+yexy )
@y
= @y
= exy + xyexy ;
@N (x;y) @ ( y1 +xexy )
@x
= @x
= exy + xyexy :
) @M@y (x;y)
= @N@x (x;y)
: ) Exact D.E.
R R
(x; y) = M (x; y) dx + A (y) = (sec2 x + yexy ) dx + A (y)
= tan x + exy + A (y) ;
R R 1
(x; y) = N (x; y) dy + B (x) = y
+ xexy dy + B (x)
= ln y + exy + B (x) :
) (x; y (x)) = tan x + exy + ln y = C:
At x = 0; y = 1: ) 0 + 1 + 0 = C: ) tan x + exy + ln y = 1 .
2 dy
Problem 2.4.10 4y 3 tan x 5 xy + 6e2y dx
= y 4 sec2 x + 10x ln y:
2
4y 3 tan x 5 xy + 6e2y dy + (y 4 sec2 x 10x ln y) dx = 0
N (x; y) dy + M (x; y) dx = 0
@M (x;y) @ (y 4 sec2 x 10x ln y )
@y
= @y
= 4y 3 sec2 x 10 xy
2
@N (x;y) @ 4y 3 tan x 5 xy +6e2y
@x
= @x
= 4y 3 sec2 x 10 xy + 0:
) @M (x;y)
@y
= @N@x
(x;y)
: ) Exact D.E.
R R
(x; y) = M (x; y) dx = (y 4 sec2 x 10x ln y) dx + A (y)
= y 4 tan x 5x2 ln y + A (y)
R R 2
(x; y) = N (x; y) dy = 4y 3 tan x 5 xy + 6e2y dy + B (x)
= y 4 tan x 5x2 ln y + 3e2y + B (x)
(x; y (x)) = y 4 tan x 5x2 ln y + 3e2y = C .
(by taking A (y) = 3e2y and B (x) = 0)
) 1
y3
=z= 1
x3
(C 3x ex + 3ex ) ) y= p
3
x
C 3x ex +3ex
:
dy
p
Problem 2.4.14 y dx = 2x 9 + y2, y (0) = 4:
Final Answer:
q Separable d.e.
p
y= (x2 + 5)2 9 = x4 + 10x2 + 16
dy y y
Problem 2.4.15 dx
= x
+ sin2 x
, y (1) = 4 :
1 e
Final Answer: Homogeneous d.e. y = x cot ln x
+ (x + y)2 = 0,
2 dy
Problem 2.4.16 2xy + x2 2yey dx
y (3) = 0:
2
Final Answer: Exact d.e. x3 + 3x2 y + 3xy 2 3ey = 24
1
Problem 2.4.17 (2y 3 e2x + cos x) dx + 3y 2 e2x + y
dy = 0; y (0) = 1:
dy
Problem 2.4.18 (4y 3 + 4ex cos y) dx = 5ex 4ex sin y:
Final Answer: Exact D.E. y 4 + 4ex sin y 5ex = C:
1
Problem 2.4.19 y + x2 + x
dx + (x sec2 y) dy = 0:
x3
Final Answer: Exact D.E. xy + 3
+ ln x tan y = C:
dy
Problem 2.4.21 y 2 + 2x sin y ex + (2y + 2xy + x2 cos y) dx = 0; y (0) = 0:
dy
Problem 2.4.23 dx
+ 3y = 7e4x , y (0) = 3:
dy
Problem 2.4.25 x2 dx 2xy = 3y 4 :
x2
Final Answer: Bernoulli D.E., y = p
3
C (9x5 =5)
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 35
dy
2. (xy + 3x2 ) dx = 2y 2 + 9xy 7x2 ; y (1) = 2:
dy
3. (3xy 2 + 6x2 y) dx = 2y 3 + 3xy 2 2x3 :
dy
4. (4xy 3 + 2x3 y) dx = 3y 4 + x2 y 2 + 5x4 :
dy
5. (4xy 3 15x2 y 2 ) dx = 2y 4 5xy 3 + 2x4 :
dy y
6. x dx = 3y ln x
+ y; y (1) = 1:
dy
7. (3xy 2 4x2 y + 3x3 ) dx = 5y 3 8xy 2 + 9x2 y:
dy
8. (3xy 2 + 3x3 ) dx = 2y 3 + 4x3 ; y (1) = 0:
dy
9. (9x2 y 2 + 4xy 3 ) dx = 2y 4 + 3xy 3 2x4 :
dy
10. (12x2 y 3 10xy 4 ) dx = 3x5 + 3xy 4 4y 5 :
dy
11. (3xy + x2 ) dx = 6y 2 + 3xy; y (1) = 2:
dy
12. x2 dx xy = y 2 :
1 dy
13. (3e2y cos 3x sin x) + 2e2y sin 3x + 1+y 2 dx
=0; y 3
= 1:
y2 1
14. 3x2 sin y x
+ 1+x2
dx + (x3 cos y 2y ln x) dy = 0; y (1) = :
3
15. (2e2x + 3x2 tan y) dx + y
+ x3 sec2 y dy = 0; y (0) = 1:
1
16. (3x2 y 2 + ey sec x tan x cos x) dx + 2x3 y + ey sec x + y
dy = 0 :
4
17. (2e2x sin 3y + sec2 x) dx + 3e2x cos 3y + y
dy = 0 :
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 36
y2 dy
18. x
sin ( x) + 2x tan y + (3y 2 + 2y ln x + x2 sec2 y) dx =0;
y (1) = 0:
2y
21. (2x cos y + sec2 x) dx + 1+y 2
x2 sin y dy = 0; y (0) = 0:
dy
22. (2 cos 2x tan y + 3e3x ) + (sin 2x sec2 y + 5y 4 ) dx =0; y (0) = 1:
3
23. y 2 ex + y x
dx + (2yex + x + cos y) dy = 0:
dy
26. dx
+ y tan x = 5 cos2 x; y (0) = 7:
dy
27. dx
+ x5 y = 6 sin x
x5
; y (1) = 4
dy 2x 2
28. dx
+ x2 16
y = x 4
; y (0) = 1:
dy
29. e3x dx + 3e3x y = 4 cos 2x:
dy
30. (csc x) dx (sec x) y = 4:
dy x
31. dx
+y =e cos x; y (1) = 2:
dy 4
32. dx
4x3 y = 3ex cos 3x; y (0) = 2:
dy
33. x dx 2y = x3 sec2 x; y( ) = 2
:
x2
34. y 0 + 2xy = 2x: Answer y = 1 + Ce :
Problem 2.5.2 RIn a series RC-circuit, the current I (t) is governed by the
equation RI + C1 I dt = E (t), where E (t) is the electric potential applied to
the circuit. When we di¤erentiate both sides, we get the di¤erential equation
dI 1 dE
R + I= :
dt C dt
Given that E (t) = E e kt and the initial current is 0, …nd I (t).
8
kE C
< kRC 1
e kt e t=RC ; if kRC 6= 1;
Final Answer: I (t) =
: kER t e t=RC ; if kRC = 1:
y 0 + p (x) y = 0: (2.23)
C
And the functions Cy1 = are all the solutions of equation (2.23).
R
(x)q(x) dx
Also, yp (x) = (x)
is one of the solutions of equation (2.20) itself;
which we obtain from formula (2.21) when we choose the value 0 for the
constant C. R
1 (x)q(x) dx
It is easy to see that these y1 = (x) and yp = (x)
can be replaced
in (2.22) by any other nonzero solutions of equations (2.23) and (2.20),
respectively.
CHAPTER 2. FIRST-ORDER DIFFERENTIAL EQUATIONS 39
equivalently,
y0 yp0 y10
= 0;
y yp y1
and so,
y1 y 0 yp0 y10 (y yp ) = 0;
that is
0 y10 y10
y y = yp0 yp :
y1 y1
We show that this new equation is exact, and we derive its solution:
@M (x;y)
@y
= @( (x)p(x)y@y (x)q(x)) = p; @N@x
(x;y)
= ddx = p:
@M (x;y) @N (x;y)
Hence, @y
= @x
; ensuring that equation (2.25) is exact.
Next, we obtain the function (x; y) such that
@ (x;y)
@x
= M (x; y) and @ @y (x;y)
= N (x; y).
R R
R (x;
d
y) = R M (x; y) dx = ( (x) p (x) y (x) q (x)) dx
= y dx dx (x) Rq (x) dx
) (x; y) = y R (x) q (x) dxR+ A (y)
Also, (x; y) = N (x; y) dy = (x) dy = y + B (x) :
By comparing these two Rbits of information on the function ,
we have to take B (x) = (x) q (x) dx and A (y) = 0:
Hence, the solution of the exact R equation (2.25) is given by
(x; y (x)) = (x) y (x) q (x) dx = C; that is,
Z
1
y= C+ (x) q (x) dx :
(x)
This is the formula, obtained in Section 2.2, for the general solution of
the linear di¤erential equation (2.4).
Chapter 3
Second-order Linear
Di¤erential Equations
d2 y dy
a (x) 2
+ b (x) + c (x) y = Q (x) : (3.1)
dx dx
The collection of all solutions of this D.E. is called its general solution yG:S: (x).
The corresponding homogeneous equation is
d2 y dy
a (x) 2
+ b (x) + c (x) y = 0: (3.2)
dx dx
The general solution of this shorter equation is called the homogeneous
solution of (3.1) (also called the complementary function of (3.1)). It is
denoted by yh (x).
We say that two functions y1 (x) and y2 (x) are independent if their ratio
y1 (x)
y2 (x)
is not a constant function.
41
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS42
Q 2y 0 b R
whereby q (x) = ; p (x) = 1 + ; (x) = e pdx , and C1 , C2 are
ay1 y1 a
arbitrary constants. Then we get V (x) by integrating V 0 (x):
R R 1R Q R 1
V (x) = V 0 (x) dx + C1 = dx dx + C2 dx + C1 :
ay1
We compute y (x) from V (x) and the substitution (3.5), and get the
following general solution of equation (3.1):
Z Z Z
1 Q 1
y (x) = V y1 = y1 (x) dx dx + C2 y1 (x) dx + C1 y1 (x) :
ay1
R 1R Q
Using the notations yp (x) = y1 (x) dx dx and
ay1
R 1
y2 (x) = y1 (x) dx, this general solution of (3.1) becomes:
When the two constants C2 ; C1 are taken to be 0; we …nd that this yp (x)
is a particular solution for (3.1). By replacing Q (x) by the constant function
0; we …nd that C2 y2 (x) + C1 y1 (x) is the required homogeneous solution yh .
y2 (x) R 1
Also, = dx; which is not constant, so y2 (x) ; y1 (x) are indeed two
y1 (x)
independent solutions of the corresponding homogeneous equation (3.2).
Finally, elementary Linear Algebra ensures that we can replace the
function yp (x), above, by any other solution of equation (3.1), and we can
replace y1 (x) and y2 (x) by any two independent solutions of equation (3.2).
This completes the proof of all claims in the theorem.
) a dx
d2 d x x
2 y1 + b dx y1 + cy1 = 0e sin x + 0e cos x = 0:
This shows that y1 = e x sin x is a solution of the homogeneous equation
(3.7). Similarly, y2 = e x cos x is a solution of (3.7). These two solutions
y1
are independent, because = tan x, which is a nonconstant function.
y2
This completes the proof.
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS46
Auxiliary Equation: m2 3m + 2 = 0:
) (m 1) (m 2) = 0 ) m1 = 1; m2 = 2:
) General solution is: yh = C1 ex + C2 e2x :
Problem 3.3.3 y 00 6y 0 + 9y = 0:
Auxiliary Equation: m2 6m + 9 = 0:
) (m 2
3) = 0 ) m1 = m2 = 3:
) General solution is: yh = (C1 + C2 x) e3x :
Auxiliary Equation: m2 + 8m + 25 = 0;
p p p
) m1;2 = b b2 4ac
2a
= 8 64 100
2
= 8
2
36
= 8 6i
2
:
) m1;2 = 4 3i = i ) = 4 and = 3:
x
and yh = e (C1 sin ( x) + C2 cos ( x)) :
) General solution is: yh = e 4x
(C1 sin (3x) + C2 cos (3x)) :
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS47
in which V1 (x) and V2 (x) are two new unknown functions, that we need to
determine.
We are going to substitute this yp (x) for y (x) in equation (3.12). This
will result in one di¤erential equation only in the two unknown functions
V1 (x) and V2 (x). This is consistent with our prior knowledge that equation
(3.12) has an in…nite number of solutions, any one of them can serve as the
required yp (x) : So, in order to come out with one solution only for yp (x) ;
we impose another di¤erential equation on the two unknowns V1 and V2 .
We choose it in such a manner that facilitates our solution greatly, without
causing any contradictions. That expedient equation is:
0 y2 y1 0
y2 Q (x) Q (x) y20 y1 Q (x) 0
y1 Q (x)
V10 = = ; V20 = = : (3.23)
W y1 y2 W y1 y2
y10 y20 y10 y20
Step 6: Obtain the two functions V1 (x) and V2 (x), by integrating V10 (x)
and V20 (x), respectively.
Step 7: Obtain a particular solution yp (x) ; as:
yp (x) = V1 y1 + V2 y2 . (3.24)
Step 8: Then write the general solution yG:S: (x) of (3.18); as:
yG:S: = yp + yh .
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS50
y1 V10 + y2 V20 = 0;
y10 V10 + y20 V20 = Q (x) :
y2 Q
V10 = W
= 12 tan x sin x cos x
sec x
= 12 sin2 x cos x:
Z Z
) V1 = V10 dx = 12 sin2 x cos x dx = 4 sin3 x:
y1 Q sec xQ
V20 = W
= = Q = 12 sin x cos x:
sec x
Z Z
) V2 = 0
V2 dx = 12 sin x cos x dx = 6 cos2 x:
Consequently,
yp (x) = V1 y1 + V2 y2 = 4 sin3 x sec x + 6 cos2 x tan x
= 4 sin2 x tan x + 4 cos2 x tan x + 2 cos2 x tan x:
) yp (x) = 4 tan x + 2 sin x cos x: Finally,
The Wronskian of y1 ; y2 is
y1 y2 e 4x e2x 2x
W = = = 6e :
y10 y20 4e 4x
2e2x
y2 Q 36xe2x e2x
V10 = W
= 6e 2x
= 6xe6x :
Z Z Z
V1 = V10 dx = 6xe 6x
dx = xe 6x
+ e6x dx:
) V1 = xe6x + 16 e6x
y1 Q 36xe 4x e2x
V20 = W
= 6e 2x = 6x
Z Z
) V2 = V20 dx = 6x dx = 3x2
) yp (x) = V1 e 4x
+ V2 e2x = xe2x + 61 e2x + 3x2 e2x
) yp (x) = 3x2 x+ 1
6
e2x : Hence;
d2 y dy
(x + 1) 2
+x 1y = (x + 1)2 ;
dx dx
x
knowing that its homogeneous solution is: yh = C1 e + C2 x:
y1 V10 + y2 V20 = 0;
y10 V10 + y20 V20 = Q (x) = x + 1:
y2 Q x(x+1)
By Cramer’s Rule, V10 = = x (x+1)) = xex
Z Z W e Z
) V1 = V1 dx =
0
xex dx = xex + ex dx = xex + ex ;
Z Z
x (x+1)
V20 = y1 Q
W
= e
e x (x+1) = 1: ) V2 = V2 dx = 1 dx = x:
0
Consequently, yp (x) = y1 V1 + y2 V2 = e x ex (1 x) + xx = 1 x + x2 :
x
yG:S: = yh + yp = C1 e + C2 x + 1 x + x2 = C1 e x
+ C3 x + 1 + x2 .
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS53
Problem 3.4.4 Verify that the function yh (x) = C1 sec2 x + C2 sec x tan x
is the homogeneous solution of the following linear di¤erential equation:
(cos x) y 00 3 (sin x) y 0 2 (cos x) y = sec2 x.
Then use this yh (x) to …nd the its general solution, by means of the
Method of Variation of Parameters.
y1 V10 + y2 V20 = 0;
y10 V10 + y20 V20 = Q (x) = sec3 x:
y2 Q sec3 x
V10 = W
= sec3 x
sec x tan x = sec x tan x
Z Z
V1 = V10 dx = sec x tan x dx = sec x:
y1 Q sec3 x
V20 = W
sec2 x = sec2 x:
= sec3 x
Z Z
) V2 = V2 dx =
0
sec2 x dx = tan x
Consequently,
yp (x) = V1 e x + V2 e 2x = sec3 x sec x tan2 x
) yp (x) = sec x (sec2 x tan2 x) = sec x: Finally,
Problem 3.4.5 Find the general solution of the following linear di¤erential
equation, using the Method of Variation of Parameters:
d2 y dy
6 + 10y = 2e3x cos x:
dx2 dx
y1 V10 + y2 V20 = 0;
y10 V10 + y20 V20 = Q (x)
y2 Q 3x 3x
V10 = Z W = e cos x2e cos x
Z e6x = 2 cos2 x = 1 + cos 2x
V1 = V10 dx = (1 + cos 2x) dx = x + 21 sin 2x = x + sin x cos x:
y1 Q e3x sin x2e3x cos x
V20 = W
= e2 x
= 2 sin x cos x
Z Z
V2 = V20 dx = 2 sin x cos x dx = sin2 x:
Consequently,
yp (x) = V1 e3x sin x+V2 e3x cos x = (x + sin x cos x) e3x sin x sin2 xe3x cos x =
xe3x sin x
Finally, yG:S: = yp + yh = (C1 + x) e3x sin x + C2 e3x cos x:
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS55
y1 V10 + y2 V20 = 0:
1
y10 V10 + y20 V20 = Q (x) = :
1 + ex
y1 y2 e x e 2x 3x
The Wronskian of y1 ; y2 is = = e :
y10 y20 e x 2e 2x
y2 Q e 2x ex
V10 == W
= e 3x (1+ex ) = 1+e x:
Z Z
) V1 = V10 dx = ex
1+ex
dx = ln (1 + ex )
y1 Q e x e2x
V20 == W
= e 3x (1+ex ) = 1+ex
Z Z Z
) V2 = V20 dx = e2x
1+ex
dx = ex
1+ex
ex dx:
We use Zthe.substitution
Z z = ex ) dz = ex dz :
) V2 = 1+z
z
dz = 1
1 + 1+z dz = z + ln (1 + z) = ex + ln (1 + ex ) :
Consequently,
yp (x) = V1 e x + V2 e 2x = e x ln (1 + ex ) + e 2x
( ex + ln (1 + ex ))
) yp (x) = (e x + e 2x ) ln (1 + ex ) e x : Finally,
) yG:S: = yp + yh = (e x
+e 2x
) ln (1 + ex ) + C3 e x
+ C2 e 2x
:
) e yp + 3e yp + 2e yp = 1+ex :
x 00 x 0 x ex
1
We divide by ex , and get yp00 + 3yp0 + 2yp = 1+e x ; as requested.
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS56
y1 V10 + y2 V20 = 0:
y10 V10 + y20 V20 = Q (x) = sin (ex ) :
The Wronskian of y1 ; y2 is
y1 y2 sin x cos x
W = 0 0 = = sin2 x cos2 x = 1:
y1 y2 cos x sin x
y2 Q
V10 = W = y2 Q = cos x sec x tan x = tan x:
Z Z
0
V1 = V1 dx = tan x dx = ln (sec x) :
y1 Q
V20 = W
= y1 Q = sin x sec x tan x
) V20 = Z tan2 x = 1Z sec2 x
) V2 = V20 dx = 1 sec2 x dx = x tan x
Consequently,
yp (x) = V1 e x + V2 e 2x = sin x ln (sec x) + x cos x cos x tan x
) yp (x) = sin x ln (sec x) + x cos x sin x: Finally,
y1 V10 + y2 V20 = 0:
y10 V10 + y20 V20 = Q (x) = xex :
The Wronskian of y1 ; y2 is
y1 y2 ex sin x ex cos x
W = =
y10 y20 ex sin x + ex cos x ex cos x ex sin x
) W = e2x sin x sin x sin2 x sin x sin x cos2 x :
) W = e2x sin2 x + cos2 x = e2x
y2 Q ex cos x (xex )
V10 = W
= 2x
= x cos x
Z Ze Z
V1 = V10 dx = x cos x dx = x sin x sin x dx = x sin x + cos x:
ex sin x (xex )
y1 Q
V20 = W
= = x sin x
e2x
Z Z Z
) V2 = V2 dx =
0
x sin x dx = x cos x cos x dx = x cos x sin x:
Problem 3.4.9 Find the general solution of the following second-order lin-
ear di¤erential equation, using the Method of Variation of Parameters:
y 00 + 3y 0 + 2y = sin (ex ) :
) yG:S: = yp + yh = e 2x
sin (ex ) + C1 e x
+ C2 e 2x
:
d2 y dy
sin2 x 2 (sin x cos x) + sin2 x + 2 cos2 x y = 2 sin3 x;
dx2 dx
y1 V10 + y2 V20 = 0;
y10 V10 + y20 V20 = Q (x) = 2 sin x:
y2 Q
By Cramer’s Rule, V10 = W
= 2 x sin x sin x
sin2 x
= 2x
Z Z
) V1 = V1 dx =
0
2xdx = x2 ;
y1 Q
V20 = = 2 sinsinx 2sinx x = 2
W
Z Z
) V2 = V2 dx = 2 dx = 2x:
0
Consequently,
yp (x) = y1 V1 + y2 V2 = x2 sin x + 2xx sin x = x2 sin x:
Finally, yG:S: = yh + yp = (C1 + C2 x + x2 ) sin x:
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS60
y1 V10 + y2 V20 = 0:
y10 V10 + y20 V20 = Q (x) = 8e3x :
The Wronskian of y1 ; y2 is
y1 y2 e3x e 5x
W = = = 5e 2x 3e 2x = 8e 2x :
y10 y20 3e3x 5e 5x
Z Z
0 y2 Q 8e 2x
V1 = W = 8e 2x = 1: ) V1 = V1 dx = 1 dx = x:
0
Z Z
0 y1 Q 8e6x 8x
V2 = W == 8e 2x = e : ) V2 = V2 dx =
0
e8x dx = 1 8x
8
e
1 3x
Consequently, yp (x) = V1 e3x + V2 e 5x
= xe3x 8
e
General Solution is yG:S: = yp + yh = xe3x + C1 e 3x
+ C2 e 5x
.
2. y 00 9y = 6e 3x
.
3x 3x
Final Answer: yG:S: = xe + C1 e + C2 e3x
3. y 00 + 9y = 6e 3x
.
1
Final Answer: yG:S: = 3
e 3x + C1 sin 3x + C2 cos 3x
4. y 00 + 6y 0 + 5y = 8ex :
Final Answer: yG:S: = 32 ex + C1 e x
+ C2 e 5x
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS61
d2 y dy
a + b + cy = Q (x) : (3.25)
dx2 dx
whose three coe¢ cients a; b; c are constants, Q (x) 6= 0 and a 6= 0, we may
obtain a particular solution yp (x) by means of the (easier) Method of
Undetermined Coe¢ cients. This consists of the following …ve steps:
Step 1: Obtain the homgeneous solution yh (x) = C1 y1 + C2 y2 of (3.25);
that is, the solution of its corresponding homogeneous equation:
d2 y dy
a 2
+ b + cy = 0: (3.26)
dx dx
Step 2: Try to put Q (x) (the right-hand side of (3.25)) in a lineaer form
2. yp (x) = An xn + An 1 xn 1
+ An 2 x n 2
+ + A2 x 2 + A1 x + A0
(polynomials of degree n), because
yp0 (x) = 0xn + nAn xn 1
+ (n 1) An 1 xn 2
+ (n 2) An 2 xn 3
+ +
2A2 x + A1 :
Particular cases are:
yp (x) = A (constants),
yp (x) = Ax + B (polynomials of degree 1),
yp (x) = Ax2 + Bx + C (polynomials of degree 2),
yp (x) = Ax3 + Bx2 + Cx + D (polynomials of degree 3).
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS62
Auxiliary Equation: m2 + 4m 12 = 0
) (m 2) (m + 6) = 0 ) m1 = 2; m2 = 6:
) Homogeneous solution is: yh = C1 e2x + C2 e 6x
:
Auxiliary Equation: m2 + 4m + 3 = 0
) (m + 3) (m + 1) = 0 ) m1 = 3; m2 = 1:
) Homogeneous solution is: yh = C1 e 3x
+ C2 e x :
Auxiliary Equation:
p mp2 + 4 = 0 p
) m1;2 = b b2 4ac
2a
= 0 0 16
2
= 6 2 144 = 24i = 2i = i
) =0 and =2
) yh = e (C1 sin ( x) + C2 cos ( x)) = C1 sin (2x) + C2 cos (2x) :
x
Auxiliary Equation:
p mp2 4m + 13 p = 0
) m1;2 = b 2
b 4ac 4 16 52 4 36
2a
= 2
= 2
= 4 26i
) m1;2 = 2 3i = i ) Homogeneous solution is:
) yh = e (C1 sin ( x) + C2 cos ( x)) = e2x (C1 sin (3x) + C2 cos (3x)) :
x
yp00 = 4Ax sin (2x) 4Bx cos (2x)+2A cos (2x) 2B sin (2x)+2A cos (2x)
2B sin (2x)
= ( 4Ax 4B) sin (2x) + ( 4Bx + 4A) cos (2x)
) yp +4yp = ( 4Ax 4B) sin (2x)+( 4Bx + 4A) cos (2x)+4Ax sin (2x)+
00
Auxiliary Equation: m2 3m + 2 = 0
) (m 1) (m 2) = 0 ) m1 = 1; m2 = 2:
) Homogeneous solution is: yh = C1 ex + C2 e2x :
Problem 3.5.8 y 00 + 4y 0 + 4y = 6e 2x
:
Auxiliary Equation: m2 + 4m + 4 = 0
) (m + 2)2 = 0 ) m1 = m2 = 2:
) Homogeneous solution is: yh = C1 e 2x
+ C2 xe 2x
:
Auxiliary Equation:
p m2p+ 6m + 45 =p0
) m1;2 = b 2ab 4ac = 6 236 180 = 6 2 144 = 6 2 12i
2
) m1;2 = 3 6i = i ) = 3 and =6
) Homogeneous solution is:
yh = e x (C1 sin ( x) + C2 cos ( x)) = e 3x (C1 sin (6x) + C2 cos (6x)) :
= (9A 36B) sin (6x)+(36A + 9B) cos (6x) = 81 sin (6x)+18 cos (6x)
) 9A 36B = 81; 36A + 9B = 18 ) A = 1. and B = 2
) yp = sin (6x) 2 cos (6x) : General Solution is
yG:S: = yp + yh = sin (6x) 2 cos (6x) + e 3x (C1 sin (6x) + C2 cos (6x)) :
Auxiliary Equation: m2 6m + 9 = 0
) (m 3) = 0
2
) m1 = m2 = 3:
) Homogeneous solution is: yh = (C1 + C2 x) e3x :
Auxiliary Equation:
p mp2 8m + 25p= 0
) m1;2 = b b2 4ac
2a
= 8 64 100
2
= 8 2 36 = 8 6i
2
yp = Ax2 + Bx + C :
Auxiliary Equation: m2 2m = 0
) m (m 2) = 0 ) m1 = 0; m2 = 2:
) Homogeneous solution is: yh = C1 + C2 e2x :
) Particular solution has form:
Problem 3.5.14 y 00 + 5y 0 + 6y = 2e 3x
; y (0) = 1; y 0 (0) = 4:
Auxiliary Equation: m2 + 5m + 6 = 0
) (m + 3) (m + 2) = 0 ) m1 = 3; m2 = 2:
) Homogeneous solution is: yh = C1 e 3x
+ C2 e 2x :
Auxiliary Equation:
p m2 + 9 = 0 ) m2 = 9
) m1;2 = 9 = 3i = i
) =0 and =3
) yh = e (C1 sin x + C2 cos x) = C1 sin 3x + C2 cos 3x :
x
1. y 00 + y = 6 sin x + 2 cos x.
Final Answer: yG:S: = (x + C1 ) sin x + (3x + C2 ) cos x:
2. y 00 2y 0 + y = 12xe3x .
Final Answer: yG:S: = C1 ex + C2 xex 3e3x + 3xe3x :
3. y 00 2y 0 + y = 12xex .
Final Answer: yG:S: = C1 ex + C2 xex + 2x3 ex :
4. y 00 4y 0 + 4y = xe3x :
Final Answer: yG:S: = C1 e2x + C2 xe2x + 2e3x xe3x :
5. y 00 + 2y 0 15y = 14e 4x
.
Final Answer: yG:S: = C1 e3x + C2 e 5x
2e 4x
:
6. y 00 + y = 6 cos x.
Final Answer: yG:S: = C1 sin x + C2 cos x + 3x sin x:
7. y 00 y = 6 cos x.
Final Answer: yG:S: = C1 ex + C2 e x
3 cos x:
8. y 00 2y 0 + y = 12xe3x .
Final Answer: yG:S: = C1 ex + C2 xex 3e3x + 3xe3x :
9. y 00 2y 0 + y = 12xex .
Final Answer: yG:S: = C1 ex + C2 xex + 2x3 ex :
12. y 00 2y 0 + y = 2 cos x :
an xn y (n) + an 1 xn 1 y (n 1)
+ + a2 x2 y 00 + a1 xy 0 + a0 y = f (x) ; (3.29)
dy d2 y dk y
where y 0 = dx
; y 00 = dx2
and y (k) = dxk
.
x = et ; (3.30)
that is,
t = ln x; (3.31)
and we let be the following di¤erentail operator:
d
= ; (3.32)
dt
we …nd that, for all natural numbers n :
!
Y1
n
n (n)
x y = ( k) y = ( 1) ( 2) ( (n 1))y: (3.33)
| {z }
k=0 n multiples
dz
) z=x : (3.34)
dx
We use Mathematical Induction:
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS75
xy 0 = y;
2
x2 y 00 = ( 1) y = y;
3 000 3 2
xy = ( 1) ( 2) y = 3 +2 y;
4 3 2
x4 y (4) = ( 1) ( 2) ( 3) y = 6 + 11 6 y:
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS76
) Particular solution takes the form yp (t) = A sin (2t) + B cos (2t) :
) yp = 2B sin (2t) + 2A cos (2t)
) 2 yp = 4A sin (2t) 4B cos (2t)
) 2 5 + 4 yp = 10B sin (2t) 10A cos (2t) = 10 sin (2t)
) A = 0; B = 1: ) yp (t) = cos (2t) :
General Solution is yG:S: (t) = yp + yh = cos (2t) + C1 et + C2 e4t .
As a function in x; this general solution becomes
yG:S: (x) = cos (2 ln x) + C1 x + C2 x4 :
2 4t
+ 8 + 25 y = 18e :
Auxiliary Equation: m2 + 8m + 25 = 0
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS77
p p p
) m1;2 = b 2ab 4ac = 8 264 100 = 8 2
2 36
= 82 6i
) m1;2 = 4 3i = i ) = 4 and =3
and yh = e t (C1 sin ( t) + C2 cos ( t)) :
yh (t) = e 4t (C1 sin (3t) + C2 cos (3t)) :
) Particular solution takes the form: yp = Ae 4t
) yp = 4Ae 4t
) 2 yp = 16Ae 4t
) 2 + 8 + 25 yp = Ae 4t (16 32 + 25) = 9Ae 4t = 18e 4t :
) A = 2: ) yp (t) = 2e 4t :
General Solution is yG:S: (t) = yp + yh = 2e 4t + e 4t (C1 sin (3t) + C2 cos (3t)) .
As a function in x; this general solution becomes
yG:S: (x) = x14 (2 + C1 sin (3 ln x) + C2 cos (3 ln x)) :
Problem 3.6.8 Find the general solution of the following linear di¤erential
equation:
d2 y dy
x2 2 5x + 5y = 8x:
dx dx
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS81
x x
Ih = C1 sin p + C2 cos p : (3.41)
LC LC
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS82
d2 I dI 2 !E
2
+2 + + !2 I = sin (!x) ; (3.50)
dx dx L
whereby ; ! denote the following two constant values:
p
R 4LC R2 C 2
= ; ! = : (3.51)
2L 2LC
Auxiliary Equation: m2 + 2 m 2 2
p + ( + ! ) = 0:
2 4 2 4 ( 2 + !2)
Its two roots are m1;2 = = i! :
2
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS84
) Homogeneous solution is
x
Ih (x) = C1 e sin (! x) + C2 e x cos (! x) (3.52)
p p
4LC R2 C 2 4LC R2 C 2
= e Rx=2L C1 sin x + C2 cos (3.53)
x :
2LC 2LC
We assume the following form for a particular solution Ip (x) (remember
that 6= 0):
Ip = A sin (!x) + B cos (!x) : (3.54)
) Ip0 = !B sin (!x)+!A cos (!x) ; Ip00 = ! 2 A sin (!x) ! 2 B cos (!x) :
Substituting in equation (3.50) we get:
Ip00 + 2 Ip0 + ( 2 + ! 2 ) Ip
= ! 2 A sin (!x) ! 2 B cos (!x) 2 !B sin (!x) + 2 !A cos (!x)
+ ( 2 + ! 2 ) A sin (!x) + ( 2 + ! 2 ) B cos (!x)
= ( ! 2 A 2 !B + ( 2 + ! 2 ) A) sin (!x)+( ! 2 B + 2 !A + ( 2 + ! 2 ) B) cos (!x)
= ( 2 !B + ( 2 + ! 2 ! 2 ) A) sin (!x)+(2 !A + ( 2 + ! 2 ! 2 ) B) cos (!x)
!E
= sin (!x) + 0 cos (!x) :
L
!E
Coe¢ cients of sin (!x): ( 2 + ! 2 ! 2 ) A 2 !B = ;
L
Coe¢ cients of cos (!x): 2 !A + ( 2 + ! 2 ! 2 ) B = 0:
We solve these two equations in A; B by means of Cramer’s Rule. Let:
( 2 + !2 !2) 2 ! 2
= = ( 2 + !2 !2) + 4 2!2:
2 ! ( 2 + !2 !2)
!E
2 !
L
0 ( 2 + !2 !2) 2
+ !2 !2 !E
)A= = 2 ;
2
( +! 2 ! ) +4 ! L
2 2 2
!E
( 2 + !2 !2)
L
2 ! 0 2 ! !E
B= = 2 :
2
( +! 2 ! ) +4 ! L
2 2 2
!E
) Ip = 2 (( 2 + ! 2 ! 2 ) sin (!x) 2 ! cos (!x)) :
L ( +!2 2 2
! ) +4 ! 2 2
Let 0 6 < 2 be the unique angle that has the following two
trigonometric ratios:
2 ! RC!
sin = q =q (3.55)
;
( 2 + !2 ! 2 )2 +4 2!2 (1 LC! 2 )2 + R2 C 2 ! 2
2
+ !2 !2 1 LC!
cos = q =q (3.56)
:
( 2 + !2 ! 2 )2 +4 2!2 (1 LC! 2 )2 + R2 C 2 ! 2
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS85
!CE
= q sin (!x + ) : (3.58)
(1 LC! 2 )2 + R2 C 2 ! 2
yG:S: = C1 y1 + C2 y2 + yp
(in which C1 and C2 are arbitrary constants, and y1 ; y2 are two independent
functions) is valid for second-order linear di¤erential equations only.
That equation has to be the following one (in the unknown function y (x)):
y y1 y2 yp y1 y2
equivalently,
y 00 yp00 y100 y200
y yp y1 y2
Also rewritten as:
0 0
10 0 0
10
y yp
B y1 C B y1 C
@ 0A =@ 0A ;
y2 y2
y1 y1
equivalently, 0 10
0
y yp
B y1 C
@ 0 A = 0:
y2
y1
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS87
Proof. Multiply both sides of equation (3.59) by the non-zero function (x) :
We get the equivalent equation:
d2 y dy
(x) 2
+ (x) b (x) + (x) c (x) y = 0: (3.63)
dx dx
We have: R R R
d d b(x) dx b(x) dx d
dx
= dx
e = e : dx b (x) dx = (x) b (x) :
But, we showed in Subsection 3.8.1 that equation (3.63) must take the
following form (by taking yp = 0; because the equation is homogeneous):
y20 y10 y200 y100 y200 y100
00 0
(x) y (x) y + (x) y = 0;
0 0
y2 y1 y2 y1 y2 y1
CHAPTER 3. SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS88
for some nonzero function (x). This means that, for some such (x),
But, we have
y20 y10 y200 y100 y20 y10 y200 y100
d
dx
= + = + 0:
y2 y1 y2 y1 y20 y10 y2 y1
It follows that 0 1
y200 y100 y20 y10
(x) = (x) b (x) = d
= d @ (x) A
dx dx
y2 y1 y2 y1
0
y y10 y20 y10
d (x) 2 d
= + (x) dx
dx
y2 y1 y2 y1
y0 y10 y200 y100
d (x) 2
= + (x) :
dx
y2 y1 y2 y1
y20 y10
d (x)
Hence, = 0:
dx
y2 y1
y20 y10
d (x)
But 6= 0; because y1 ; y2 are independent. Therefore, =
dx
y2 y1
0; and we can take (x) any non-zero constant. We take (x) = 1, and we
y20 y10
…nd that = (x) :
y2 y1
This proves that y2 (x) is any solution of the following …rst-order linear
di¤erential equation in the unknown function y (x):
y0 y10
= (x) : (3.64)
y y1
d3 y d2 y dy
3
+ b (x) 2
+ c (x) + k (x) y = Q (x) ; (3.66)
dx dx dx
whereby we know the general solution
We solve this system of algebraic equations. Then V1 ; V2 ,V3 are obtained from
these derivatives V10 ; V20 ,V30 ; by integration, and yp is obtained from (3.69).
Proposition 12 Let f be a real function, and let T > 0 and s be two real
constants such that the following inequality holds for all t > TR : jf (t)j es t :
1
Then for all values s > s (of a constant s), the value 0 e st f (t) dt
is …nite (that is, it is not 1).
Proof. We have:
R 1 st R1 R 1 j f (t)j
T
e f (t) dt 6 T je st
f (t)j dt = T
dt
est
91
CHAPTER 4. THE LAPLACE TRANSFORM 92
R 1 es t R 1 (s s )t 1 t=1
6 T st
dt = T
e dt = e (s s )t t=T
e s s
1 (s s )1 (s s )T
= e e :
s s
1 1
But, e (s s )1 = (s s )1 = = 0 for all s > s : Hence,
e 1
R 1 st 1
T
e f (t) dt 6 e (s s )T , and so it is …nite.
s s
0 when t < 0
U (t) = :
1 when t 0
CHAPTER 4. THE LAPLACE TRANSFORM 94
n o
dF (s) df (t)
L ft f (t)g = ds
L dt
= sF (s) f (0)
n o
dn F d2 f (t)
L ftn f (t)g = ( 1)n dsn
L dt2
= s2 F (s) sf (0) f 0 (0)
R1 nR o
t
L 1t f (t) = s F (v) dv; L 0
f (z) dz = 1s F (s)
whenever one of them exists.
n o R
1
Rule 2: L dfdt(t) = 0 e st dfdt(t) dt
t=1 R 1 @e st
= [e st f (t)]t=0 f (t) dt (integration by parts)
R 1 0 st @t
= 0 f (0) + s 0 e f (t) dt (because s > s )
n o
) L dfdt(t) = sF (s) f (0) :
Applying this rule twice, we get
n 2 o n o n o
L d dtf 2(t) = L dtd dfdt(t) = sL dfdt(t) f 0 (0)
= s (sF f (0)) f 0 (0) = s2 F sf (0) f 0 (0) :
CHAPTER 4. THE LAPLACE TRANSFORM 95
1 !
1 = U (t) s
(s > 0) sin (!t) s2 +! 2
(s > 0)
1
t s2
(s > 0)
2 s
t2 s3
(s > 0) cos (!t) s2 +! 2
(s > 0)
6
t3 s4
(s > 0)
24 k
t4 s5
(s > 0) sinh (kt) s2 k 2
(s > jkj)
n! s
tn (n = 0; 1; 2; ) sn+1
(s > 0) cosh (kt) s2 k 2
(s > jkj)
1
ekt s k
(s > k) (t) 1
sin(3t)
Problem 4.3.2 t
sin(3t) R1 R1 3 3 1 v v=1
L t
= L (sin (3t)) = v 2 +9
dv = 3
tan 3 v=s
s s
1 1 s 1 s 1 s
= tan (1) tan 3
= 2
tan 3
= cot 3
:
d2 (t cos 5t)
Problem 4.3.3 h (t) = dt2
s2 +25 s(2s)
L (t cos 5t) = d (cot
ds
s5)
= d
ds
s
s2 +25
= (s2 +25)2
s(2s) (s +25)
2 2 25
= (s2 +25)2 = (ss2 +25) 2;
s s 3
L (e3t cos 4t) = s2 +16 js
= (s 3)2 +16
!s 3
d(e 3tcos 4t) s2 3s
L (k) = L dt
= sL (e3t cos 4t) (e3t cos 4t)jt=0 = s2 6s+25
1
s2 3s (s2 6s+25)
) L (k) = s2 6s+25
= 3s 25
s2 6s+25
.
CHAPTER 4. THE LAPLACE TRANSFORM 98
Rt
Problem 4.3.5 g (t) = z 10 e 4z
dz
0
L (t10 e 4t ) = s10!
11 j
10!
= (s+4) 11
s !s+4
Rt 10 4z
)L z e dz = s(s+4) 10!
11 :
0
L (t 5)10 e 4(t 5)
U (t 5) = L (t10 e 4t
)e 5s
= 10!
(s+4)11
e 5s :
2(t 3)
Problem 4.3.8 g (t) = (t 3) e sin 4 (t 3) U (t 3)
d d 4 8s
L (t sin (4t)) = ds
L (sin (4t)) = ds s2 +16
= (s2 +16)2
) L (e 2t
t sin (4t)) = 8(s)
2 = 8(s+2)
2
((s)2 +16) js !s+2
((s+2)2 +16)
) L (g) = L e 2(t 3)
(t 3) sin 4 (t 3) U (t 3) = e 3s
L (e 2t
t sin (4t))
) L (g) = 8(s+2)
2 e 3s
.
((s+2)2 +16)
4s 36
L (k) = 7 2e + s2 +16
e 5s
s(1+e s
)
Final Answer: s2 +1
:
CHAPTER 4. THE LAPLACE TRANSFORM 99
4.3.2 Convolution
De…nition 3 The convolution of two functions f and g is the function f ? g
de…ned as follows:
Rt
(f ? g) (t) = 0
f (z) g (t z) dz
2. Associativity: f ? (g ? h) = (f ? g) ? h
3. Linearity: f ? (g + h) = (f ? g) + (f ? h) ;
f ? (cg) = c (f ? g) = (cf ) ? g when c is a constant
Rt
4. 1 is not a neutral element for ?: (f ? 1) (t) = 0 f (z) dz
Convolution Theorem. For any two functions f (t) and g (t) which
have Laplace transforms:
L ff ? gg = L ff g.L fgg :
0 t < 1; 0 z t:
We use the substitution t = z+y, and get the following region in the yz-plane:
0 y < 1; 0 z < 1:
d 12(s 3)
L (k) = ds
L (e3t sin 6t) = 2 (Problem 4.3.23):
((s 3)2 +36)
1 6 6
L (g) = L (t) :L (e3t sin 6t) = s2 s2 +36 js
= s2 ((s 3)2 +36)
:
!s 3
6
L (h) = L (t) :L (e3t ) :L (sin 6t) = s2 (s 3)(s2 +36)
.
CHAPTER 4. THE LAPLACE TRANSFORM 101
4! 4!
L (z 4 e3z ) = L (z 4 )js !s 3
= s 5 js
= (s 3)5
!s 3
Rt
) L (g) = L z 4 e3z dz = 1s L (z 4 e3z ) = 4!
s(s 3)5
.
0
Rt 2z
Problem 4.3.13 e cosh (5z) dz;
0
2z (s) s+2
L (e cosh (5z)) = L (cosh (5z))js !s+2
= (s)2 25
= (s+2)2 25
js !s+2
Rt
)L e 2z
cosh (5z) dz = 1s L (e 2z
cosh (5z)) = s+2
s((s+2)2 25)
.
0
Rt
Problem 4.3.14 g (t) = z sin 3z dz
0
d d 3
L (z sin 3z) = ds L (sin 3z) = ds s2 +9
= (s26s
+9)2
.
Rt
) L (g) = L z sin 3z dz = 1s L (z sin 3z) = (s2 +9)
6
2 :
0
2t s s+2
L (e cos 5t) = s2 +25 js
= (s+2)2 +25
!s+2
d(e 2t cos 5t) 2t 2t s2 +2s
L (k) = L dt
= sL (e cos 5t) (e cos 5t)jt=0 = s2 +4s+29
1
) L (k) = s2 +2s s2 4s 29
s2 +4s+29
= 2s 29
s2 +4s+29
.
d
Problem 4.3.16 r (t) = dt
(t sin 6t)
d d 6 12s
L (t sin 6t) = ds
(L (sin 6t)) = ds s2 +36
= (s2 +36)2
1 cos 4(t 6)
Problem 4.3.17 k (t) = t 6
U (t 6) :
1 cos 4t
R1 R1 1 v 1 v=1
L t
= L (1 cos 4t) = v v 2 +16
dv = ln (v) 2
ln (v 2 + 16) v=s
s s
p v=1
h iv=1
= ln (v) ln v2 + 16 = ln p v = ln 1 ln p s :
v=s v 2 +16 v=s s2 +16
p
p s s2 +16
=0 ln s2 +16
= ln s
.
p
) L (k) = e 6s
ln s2 +16
s
:
) L (k) = L e3(t 7)
(t 7) cos 2 (t 7) U (t 7) = e 7s
L (e3t t cos (2t))
(s 3)2 4
) L (k) = 2 e 7s
= s2 6s+5
(s2 6s+13)2
e 7s .
((s 3)2 +4)
L sin2 t = L 1
2
1
2
cos 2x = 1
2s
s
2(s2 +4)
= 2
s(s2 +4)
.
) L (k) = L sin2 t L (sin t) L (sin t) = 2
s(s2 +4)
1
(s2 +1)2
.
Problem 4.3.22 Draw a neat sketch for the following function, and …nd its
Laplace transform:
g (t) = cos t cos t U t 2
) L (h) = 12(s 3)
2 + 12s
(s 3)(s2 +36)2
+ 3e 2s
.
((s 3)2 +36)
Rt Rt
Problem 4.3.24 g (t) = ze7z sin 4z dz + ((ze7z ) sin 4z) dz
0 0
1 8s 56
As in Problem 4.3.23, L (te7t ) = (s 7)2
; L (te7t sin 4t) = 2
((s 7)2 +16)
) L (q) = 1
2
7
s2 +49
3
s2 +9
+ 5s
(s2 +25)(s2 +4)
.
CHAPTER 4. THE LAPLACE TRANSFORM 104
20 20 20
L (h) = (s)2 25
= (s 3)2 25
= s2 6s 16
:
js !s 3
2(t 6)
Problem 4.3.27 k (t) = e cos 3 (t 6) U (t 6) :
2t s s+2
L (e cos 3t) = L (cos 3t)js !s+2
= s2 +9 js
= (s+2)2 +9
!s+2
)L e 2(t 6)
cos 3 (t 6) U (t 6) = e 6s
L (e 2t
cos 3t)
) L (k) = s+2
(s+2)2 +9
e 6s :
) L (g) = e 3s
L (e2t sin 4t) = 4e 3s
(s 2)2 +16
:
Rt
Problem 4.3.29 g (t) = z 6 e3z dz + (t 2)6 e3(t 2)
U (t 2)
0
6! 6!
L (t6 e3t ) = L (t6 )js !s 3
= s 7 js
= (s 3)7
:
!s 3
d d s s2 9
L (t cos 3t) = ds
L (cos 3t) = ds s2 +9
= (s2 +9)2
:
) L (p) = s2 9
(s2 +9)2
e s=2
.
CHAPTER 4. THE LAPLACE TRANSFORM 105
7(t 4)
Problem 4.3.31 k (t) = (t 4) e sin 6 (t 4) U (t 4) :
d 6 12s
L (t sin 6t) = ds s2 +6
= (s2 +36)2
:
) L (t e 7t
sin 6t) = 12s
(s2 +36)2
= 12(s+7)
2 :
js !s+7 ((s+7)2 +36)
) L (k) = e 4s
L (t e 7t
sin 6t) = 12(s+7)
2 e 4s
:
((s+7)2 +36)
4 4 4
= s2 +16 js
= (s 6)2 +16
= s2 12s+52
!s 6
) L (h) = L e6(t 5)
sin 4 (t 5) U (t 5) = e 5s
L (e6t sin 4t) = 4e 5s
s2 12s+52
.
2(t 6)
Problem 4.3.33 k (t) = e 3et 6
U (t 6)
6s 2t 1 3 ( 2s 7)
L (k) = e L (e 3et ) = s+2 s 1
e 6s
= (s 1)(s+2)
e 6s
) L (t 2) e9(t 2)
U (t 2) = e 2s
L (t e9t ) = e 2s
(s 9)2
:
1 e4t
Problem 4.3.35 t
Z1 Z1
e4t
L 1
t
= L (1 e4t ) = 1
v
1
v 4
dv = [ln (v) ln (v 4)]v=1
v=s
s s
v v=1 s s 4
= ln v 4 v=s
= ln 1 ln s 4
= ln s
:
CHAPTER 4. THE LAPLACE TRANSFORM 106
Rt Rt
Problem 4.3.36 g (t) = 0
e6z sin 4z dz 0
(e6z sin 4z) dz
4 4
L (e6t sin 4t) = s2 +16 js
= (s 6)2 +16
:
!s 6
Rt Rt Rt
Problem 4.3.37 g (t) = e3z cos (4z) dz+ z cos (4z) dz (z cos (4z)) dz:
0 0 0
1. t t = 16 t3 :
1
2. t tn = (n+1)(n+2)
tn+2 :
) L (h) = e 6s
L (t2 + 4t + 4) = e 6s 2
s3
+ 4
s2
+ 4
s
= 2+4s+4s2
s3
e 6s .
CHAPTER 4. THE LAPLACE TRANSFORM 107
Problem 4.3.43 te2t sin 4t + t (e2t sin 4t) + e2t (t sin 4t) :
8(s 2) 4 8s
Final Answer: 2 + s2 ((s 2)2 +16)
+ (s 2)(s2 +16)2
((s 2)2 +16)
2t
Problem 4.3.45 1 3e + 2e2t cos 3t + e2t sin 3t:
1 3 2s 1
Final Answer: s s+2
+ (s 2)2 +9
d
Rt
Problem 4.3.47 dt
(t9 e 5t
)+ v9 e 5v
dv + (t9 e 5t
) cos 7t:
0
9!(s) 9! 9!(s)
Final Answer: (s+5)10
+ s(s+5)10
+ (s+5)10 (s2 +49)
Problem 4.3.48 t9 (e 5t
cos 7t) + t9 e 5t
cos 7t:
9!(s+5) 9!(s)
Final Answer: s10 ((s+5)2 +49)
+ s10 (s+5)(s2 +49)
p p
) L (g) = e s=6
L 1
2
sin 2t + 2
3
cos 2t = -e s=6 2+ 3s
2(s2 +4)
.
h (t) = cos t+ 6 6
U t 6
= cos t + 6 jt
:
!t 6
p
3 1
But cos t + 6
= cos (t) cos 6
sin t sin 6
= 2
cos (t) 2
sin (t)
p
) L cos t + 6
= 1+ 3s
2(s2 +1)
) L (h) = L cos t + 6 jt
=e s=6
L cos t + 6
!t 6
p
) L (h) = 1+ 3s
2(s2 +1)
e s=6
.
Problem 4.3.53 g (t) = t e2t cos 6t + t (e2t cos 6t) + (t e2t ) cos 6t
+ (t cos 6t) e2t + t cos 6t e2t :
Final Answer:
(s 2)2 36 s 2 s s2 36 s
L (g) = 2 + s2 2 + (s 2)2 (s 2 +36) + (s2 +36)2 (s 2) + s2 (s2 +36)(s 2) :
((s 2
2) +36) ((s 2) +36)
CHAPTER 4. THE LAPLACE TRANSFORM 109
1
L (M ) = 4 (t) + 3 (t 2) 7 (t 6) + 5U (t 9) .
3s+41
Problem 4.4.3 H (s) = (s 3)(s+7)
e 9s
)L 1
(H) = (5e3t 2e 7t
)jt !t 9
= 5e3(t 9)
2e 7(t 9)
U (t 9)
16+s2
Problem 4.4.4 G (s) = ln s2
)L 1 s 3
s2 6s+13
=L 1 s
s2 +4 js !s 3
= e3t L 1 s
s2 +4
= e3t cos 2t
Also, L 1 s2s+4 e 7s = L 1 s2s+4 j = cos 2tjt !t 7
t !t 7
)L 1
(F ) = e3t cos 2t + cos 2 (t 7) U (t 7) .
s 1
Problem 4.4.6 F (s) = s2 2s+10
e 4s :
1 s 1 1 s 1 1 s
L s2 2s+10
=L (s 1)2 +9
=L s2 +9 js !s 1
t 1 s t
=eL s2 +9
= e cot 3t
)L 1
(F ) = et 4
cos 3 (t 4) U (t 4) .
s 3
Problem 4.4.7 K (s) = s(s2 6s+25)
C(s 3)
By partial fractions K (s) = s(s2 s 6s+25)
3
= As + B
s2 6s+25
+ s2 6s+25
) s 3 = A (s2 6s + 25) + Bs + C (s 3) s
at s = 0 : 3 = 25A ) A = 3=25;
at s = 3 : 0 = 16A + 3B ) B = 16=25;
coe¢ cients of s2 : 0=A+C = ) C = 3=25:
) K (s) = 1
25
3 1s + 4 (s 4
3)2 +16
+ 3 (s s 3
3)2 +16
1
= 25
3 1s + 4 s2 +16
4 s
+ 3 s2 +16 js !t 3
)L 1
(K) = 1
25
( 3 + 4e3t sin 4t + 3e3t cos 4t) .
(Compare with Problem 4.3.37.)
3
Problem 4.4.8 F (s) = s2 4s+13
e 7s
1 3 1 3 1 3
L s2 4s+13
=L (s 2)2 +9
=L s2 +9 js !s 2
3
= e2t L 1
s2 +9
= e2t sin 3t
)L 1
(F ) = (e2t sin 3t)jt!t 7
= e2(t 7)
sin 3 (t 7) U (t 7) .
CHAPTER 4. THE LAPLACE TRANSFORM 111
s2 +4s+1
Problem 4.4.9 H (s) = (s+3)2 (s+2)
e 7s :
) s2 + 4s + 1 = A (s + 2) + B (s + 3) (s + 2) + C (s + 3)2
at s = 3 : 2= A ) A = 2;
at s = 2 : 3 = C;
coe¢ cients of s2 : 1 = B + C = B 3; ) B = 4:
) s2 +4s+1
(s+3)2 (s+2)
= 2
(s+3)2
+ 4
s+3
3
s+2
)L 1 s2 +4s+1
(s+3)2 (s+2)
= 2te 3t
+ 4e 3t
3e 2t
)L 1
(H) = L 1 s2 +4s+1
(s+3)2 (s+2)
= (2te 3t
+ 4e 3t
3e 2t
)jt !t 7
jt !t 7
)L 1
(H) = 2 (t 7) e 3(t 7)
+ 4e 3(t 7)
3e 2(t 7)
U (t 7) .
1 s
Problem 4.4.10 F (s) = cot 7
1 s
Let L (f ) = F (s) = cot 7
) tf (t) = L 1 7
s2 +49
= sin 7t
)L 1
(F ) = f (t) = sin 7t
t
.
CHAPTER 4. THE LAPLACE TRANSFORM 112
)L 1
(F ) = (3 cos 4t 4 sin 4t) e3t . (Compare with Problem 4.3.4.)
2s 29
Problem 4.4.12 F (s) = s2 +4s+29
)L 1
(F ) = ( 2 cos 5t 5 sin 5t) e 2t
(Compare with Problem 4.3.15.)
4
Problem 4.4.13 Q (s) = (s 2)2 +16
e 3s
1 4 1 4
L (s 2)2 +16
=L s2 +16 js
= e2t sin 4t
!s 2
)L 1
(Q) = e2(t 3)
sin 4 (t 3) U (t 3) .
20
Problem 4.4.14 G (s) = (s 8)(s+2)
e 5s
20 A B
By partial fractions (s 8)(s+2)
= s 8
+ s+2
) 20 = A (s + 2) + B (s 8)
at s = 2 : 20 = 10B )B= 2;
at s = 8 : 20 = 10A ) A = 2:
) 20
(s 8)(s+2)
= 2s 18 1
2 s+2 :
)L 1 20
(s 8)(s+2)
= 2e8t 2e 2t
(Compare with Problem 4.3.26.)
)L 1
(G) = 2e8(t 5)
2e 2(t 5)
U (t 5) .
CHAPTER 4. THE LAPLACE TRANSFORM 113
4s2 24s+36
Problem 4.4.15 H (s) = (s 2)2 (s 4)
4s2 24s+36 A B C
(s 2)2 (s 4)
= s 2
+ (s 2)2
+ s 4
(partial fractions)
) 4s 2
24s + 36 = A (s 2) (s 4) + B (s 4) + C (s 2)2
at s = 2 : 4 = 2B ) B = 2;
at s = 4 : 4 = 4C ) C = 1;
coe¢ cients of s2 : 4 = A + C; ) A = 4 C = 3:
) H (s) = s 2 (s 2)2 + s 4
3 2 1
) L 1 (H) = (3 2t) e2t + e4t .
s+4
Problem 4.4.16 F (s) = s2 +8s+41
e 2s
1 s+4 s+4 s
L s2 +8s+41
= L 1 (s+4) 2
+25
= L 1 s2 +25 js !s+4
s
= e 4t L 1 s2 +25 = e 4t cot 5t
)L 1 s+4
(F ) = L 1 s2 +8s+41 jt !t 2
= e 4(t 2) cos 5 (t 2) U (t 2) .
2s 7
Problem 4.4.17 W (s) = (s 1)(s+2)
e 6s
2s 7 A B
We use partial fractions: (s 1)(s+2) = s 1
+ s+2
) 2s 7 = A (s + 2) + B (s 1)
at s = 1 : 9 = 3A ) A = 3;
at s = 2 : 3 = 3B ) B = 1;
) 2s 7
(s 1)(s+2)
= 3s 11 + 1
s+2
)L 1 2s 7
(s 1)(s+2)
= 3et + e 2t
)L 1
(W ) = L 1 2s 7
(s 1)(s+2)
= (e 2t
3et )jt !t 6
jt !t 6
)L 1
(W ) = e 2(t 6)
3et 6
U (t 6)
s+2
Problem 4.4.18 F (s) = (s+2)2 9
e 4s
1 s+2 1 s
L (s+2)2 9
=L s 2 9 js !s+2
2t 1 s 2t
=e L s2 9
=e cosh 3t
)L 1
(F ) = (e 2t
cosh 3t)jt!t 4
=e 2(t 4)
cosh 3 (t 4) U (t 4) .
CHAPTER 4. THE LAPLACE TRANSFORM 114
X
1
1 3ns
Problem 4.4.19 1 e 3s = e (in…nite geometric series)
n=0
3Ks
X1
K
1 e 3ns
1 e 3s = e (K is a positive integer, …nite geometric series),
n=0
X
1 X1
K
)L 1 1 e 3Ks
1 1
1 e 3s = (t 3n) , L 1 e 3s = (t 3n) .
n=0 n=0
13. T (s) = 7s e 9s
+ 5 + 4e 2s
:
CHAPTER 4. THE LAPLACE TRANSFORM 115
s2 +5s 6
14. V (s) = (s 2)2 (s+6)
:
6s2 +7s 2
15. W (s) = (s 2)(s+1)2
e 4s :
s 4
16. X (s) = s2 8s+25
e 7s :
13
17. Y (s) = (s+2)(s2 +9)
e 6s :
7!
18. Z (s) = (s 3)8
e 4s :
s+2
19. A (s) = s2 2s+2
:
6s2 +34s+76
20. C (s) = (s 3)(s2 +8s+25)
.
8+8s2
21. D (s) = (s+1)2 (s 3)
:
s 7
22. E (s) = s2 +s 2
e 5s +e 2s
:
4s2 24s+36
23. G (s) = (s 2)2 (s 4)
:
s 2
24. H (s) = s2 4s 5
e 6s +e 4s
:
5s+5
25. J (s) = (s 1)(s2 +4s+5)
:
6s s 4
26. K (s) = 1 2e + s2 8s+20
e 3s :
8!
27. N (s) = 1 + (s+4)9
e 2s :
1
28. P (s) = 2
tan (s).
s+7 1 e 7t
29. R (s) = ln s
: Answer: t
:
s2 +4 2e 3t 2
30. Q (s) = ln s2 +6s+13
: Answer: t
cos 2t:
s s+7
31. Y (s) = ln s+7
+ ln s
:
s2 9
32. T (s) = ln s2
:
)F s2
s2 1
= 2s2 +3s 2
s(s2 1)
)F = 2s2 +3s 2
s3
= 2 1s + 3 s12 2
s3
) f (t) = L 1
(F ) = 2 + 3t t2 :
Rt
Problem 4.5.3 f (t) = e3t e2t + 0
f (z) e2(t z)
dz
f (t) = e3t e2t + f ? e2t :
Denote L (f ) (s) by F (s). ) F = s 1 3 1
s 2
+ F: s 1 2
) F 1 s 1 2 = s 1 3 s 1 2 = (s(s 2)3)(s(s 2)3)
) F ss 32 = (s 3)(s
1
2)
)F = s 2
s 3
1
(s 3)(s 2)
= 1
(s 3)2
= 1
s 2 js !s 3
) f (t) = L 1
(F ) = te . 3t
CHAPTER 4. THE LAPLACE TRANSFORM 117
Rt
Problem 4.5.4 f (t) = 8te4t 4 0
f (z) e4(t z)
dz
Rt
Problem 4.5.6 f (t) = 3t2 + et 0
f (z) et z
dz:
f (t) = 3t2 + et (f ? et ) :
Denote L (f ) (s) by F (s), and apply the Laplace transform to both sides
of the equation. Then we get by the Convolution Theorem:
F = s63 + s 1 1 F: s 1 1
) F 1 + s 1 1 = s63 + s 1 1 i.e., F s s 1 = s63 + s 1 1
) F = s s 1 s63 + s 1 1 = 6(ss4 1) + s(ss 11) = s63 s64 + 1s
) f (t) = L 1
(F ) = 3t2 t3 + 1 .
CHAPTER 4. THE LAPLACE TRANSFORM 118
Rt
Problem 4.5.7 f (t) = et cos t + 3 0
f (z) e 2(t z)
dz: .
)F = s 1
(s 1)2 +1
1
+ 3F: s+2 : )F 1 3
s+2
= s 1
(s 1)2 +1
:
)F s 1
s+2
= s 1
(s 1)2 +1
: F = s+2
(s 1)2 +1
)F = s 1
(s 1)2 +1
+ 3
(s 1)2 +1
= s
s2 +1
+ 3
s2 +1 js
:
!s 1
) f (t) = L 1
(F ) = et (cos t + 3 sin t) .
Rt
Problem 4.5.8 3f (t) = e7t e4t 9 0
f (z) e7(t z)
dz:
)F 3(s 4)
s 7
= 3
(s 7)(s 4)
)F = s 7
s 4
1
(s 7)(s 4)
= 1
(s 4)2
)F = 1
s2 js !s 4
) f (t) = L 1 4t
(F ) = te .
Rt
Problem 4.5.9 f (t) = 3 sin 5t 8 0
f (z) cos 5 (t z) dz .
)F s2 +8s+25
s2 +25
= 15
s2 +25
)F = 15
s2 +8s+25
= 15
(s+4)2 +9
=5 3
s2 +9 js !s+4
) f (t) = L 1
(F ) = 5e 4t
sin 3t .
CHAPTER 4. THE LAPLACE TRANSFORM 119
t t 2t
Rt 2(t z)
Problem 4.5.10 f (t) = te e +e + 0
f (z) e dz
f (t) = te t e t + e 2t + f ? e 2t :
Denote L (f ) (s) by F (s).
) F = (s+1)
1
2
1
s+1
1
+ s+2 1
+ F: s+2
s+2 (s+1)(s+2)+(s+1)2
)F 1 1
s+2
= (s+1)2 (s+2)
1
)F s+1
s+2
= 2
(s + 1) (s + 2)
1 1
) F = s+2s+1 2 = 1
(s+1)3
= 2
s 3 js
(s + 1) (s + 2) 2 !s+1
1
) f (t) = L 1 (F ) = t2 e t .
2
Rt
Problem 4.5.11 f (t) = e5t sinh 3t 6 0
f (z) e8(t z)
dz:
Rt
Problem 4.5.12 f (t) = 4 sin 5t 6 f (z) cos 5 (t z) dz:
0
Rt
Problem 4.5.13 f (t) = e3t e2t + 0
f (z) e2(t z)
dz:
Rt
Problem 4.5.14 f (t) = 8e3t 2 0
f (z) cos (t z) dz:
Rt
Problem 4.5.15 f (t) = e2t + sin 2t 4 0
f (z) e2(t z)
dz:
Rt
Problem 4.5.16 f (t) = 4tet + 2 0
f (z) et z
dz:
Rt
Problem 4.5.17 f (t) = sinh 2t 4 0
f (z) e2(t z)
dz:
Rt
Problem 4.5.18 f (t) = sin 3t 6 0
f (z) cos 3 (t z) dz:
Rt
Problem 4.5.19 f (t) = cos 2t + sin 2t 4 0
f (z) e2(t z)
dz:
CHAPTER 4. THE LAPLACE TRANSFORM 120
Problem 4.6.2 y 00 + 4y 0 + 3y = 6e 3t
; y (0) = 1; y 0 (0) = 2
) s2 + 5s + 12 = A (s + 3) (s + 1) + B (s + 1) + C (s + 3)2
at s = 3 : 6 = 2B ) B = 3;
at s = 1 : 8 = 4C ) C = 2;
coe¢ cients of s2 : 1 = A + C; ) A = 1 C = 1:
) Y (s) = s+3 (s+3)2 + s+1
1 3 2
) y (t) = L 1
(Y ) = e 3t
3te 3t
+ 2e t
= (1 + 3t) e 3t
+ 2e t
.
CHAPTER 4. THE LAPLACE TRANSFORM 121
) Y (s) = s2
3
4s+5
= (s 2)3 2 +1 = 3 1
s2 +1 js !s 2
) y (t) = L 1
(Y ) = 3e2t sin t .
) y (t) = L 1
(Y ) = e3t + 12 t4 e3t = 1 + 21 t4 e3t .
) 3s2 + 4 = A (s 2) (s + 2) + B (s + 2) + C (s 2)2
at s = 2 : 16 = 4B ) B = 4;
at s = 2 : 16 = 16C ) C = 1;
coe¢ cients of s2 : 3 = A + C; ) A = 2:
) Y (s) = s 2 + (s 2)2 + s+2
2 4 1
) y (t) = L 1
(Y ) = 2e2t + 4te2t + e 2t
.
and so ) y (t) = L 1
(Y ) = (3 2t) e2t + e4t .
CHAPTER 4. THE LAPLACE TRANSFORM 123
d2 y
Problem 4.6.8 dt2
+ 2 dy
dt
8y = 36te2t ; y (0) = 0; y 0 (0) = 1:
) (s2 + 2s 8) Y = 36
(s 2)2
+1
36+(s2 4s+4)
) (s + 4) (s 2) Y = (s 2) 2 = s2 4s 32
(s 2)2
)Y = (s 8)(s+4)
(s 2)3 (s+4)
= s 8
(s 2)3
= s 2
(s 2)3
6
(s 2)3
= 1
(s 2)2
6
(s 2)3
)Y = 1
s2
3 s23 js !s 2
)y=L 1
(Y ) = (t 3t2 ) e2t :
Problem 4.6.9 y 00 + 5y 0 + 6y = 2e 3t
; y (0) = 1; y 0 (0) = 4:
) (s + 3) (s + 2) Y = (s+3)(s+1) 2
s+3
= s2 +4s+1
s+3
)Y = s2 +4s+1
(s+3)2 (s+2)
= A
(s+3)2
+ B
s+3
+ C
s+2
(partial fractions)
) s + 4s + 1 = A (s + 2) + B (s + 3) (s + 2) + C (s + 3)2
2
at s = 3 : 2= A ) A = 2;
at s = 2 : 3 = C;
coe¢ cients of s2 : 1 = B + C = B 3; ) B = 4:
)Y = s2 +4s+1
(s+3)2 (s+2)
= 2
(s+3)2
+ 4
s+3
3
s+2
) y (t) = L 1
(Y ) = 2te 3t
+ 4e 3t
3e 2t
.
(See also Problem 3.5.14.)
CHAPTER 4. THE LAPLACE TRANSFORM 124
Problem 4.6.20 y 00 + 5y 0 + 6y = 2e 3t
; y (0) = 1; y 0 (0) = 4:
Problem 4.6.26 y 00 + y 0 2y = 9e 2t
; y (0) = 0; y 0 (0) = 3.
Problem 4.6.28 y 00 9y = 6e 3t
; y (0) = 1; y 0 (0) = 4 .
d2 y
Problem 4.6.30 dt2
3 dy
dt
+ 2y = 4e2t ; y (0) = 3; y 0 (0) = 5:
Answer: y (t) = 3et + 6e2t 4te2t :
CHAPTER 4. THE LAPLACE TRANSFORM 126
6
) 6X + sY 3= : (4.2)
s+1
Also, (s + 1) X Y = 0; and when we multiply by s,
s2 + s X sY = 0; (4.3)
dy
L dx
dt
= sX x (0) = sX 1; L dt
= sY y (0) = sY + 1:
Apply the Laplace transform:
CHAPTER 4. THE LAPLACE TRANSFORM 127
) sX 1 6X 3Y = 0; 3X + sY + 1 = 2
(s 3)2
: This gives
dx dy
Problem 4.7.3 dt
+y = 2t; x 2y + dt
=t;
x (0) = 0; y (0) = 1 :
Problem 4.7.4
dx
2y = 6e3t sin t;
dt
dx dy
10x 3 +2 = 0;
dt dt
x (0) = 0; y (0) = 1:
6
sX 2Y = ; (4.8)
(s 3)2 + 1
and 10X 3sX + 2 (sY 1) = 0; that is,
)X= 2
s2 6s+10
= 2
(s 3)2 +1
=2 1
s2 +1 js!s 3
) x (t) = L 1
(X) = 2e3t sin t .
dx
As dt
2y = 6e3t sin t; then
dx
2y = dt
6e3t sin t = dtd (2e3t sin t) 6e3t sin t
= 6e3t sin t + 2e3t cos t 6e3t sin t = 2e3t cos t
dy dx
Problem 4.7.5 x dt
= 1; dt
2x + y = t 2,
x (0) = 2; y (0) = 0 .
1
)X sY = ; (4.11)
s
1 2
Also, sX 2 2X + Y = s2 s
; and so,
1 2s + 2s2
) s2 X 2sX + sY = : (4.12)
s
2s+2s2 2s(s 1)
We add (4.11) and (4.12): (s2 2s + 1) X = s
= s
) (s 1)2 X = 2 (s 1) ) X = s 2 1 ) x (t) = L (X) = 2e . 1 t
As dx
dt
2x + y = t 2; then y = t 2 dx dt
+ 2x
) y = t 2 2e + 4e
t t
) y (t) = t 2 + 2et .
Problem 4.7.6 4x 2 dy
dt
= 6e3t 8te3t ; dx
dt
3x + 2y = 8 ;
x (0) = 0; y (0) = 3 :
dy
L dx
dt
= sX x (0) = sX and L dt
= sY y (0) = sY 3:
We apply the Laplace transform: ) 4X 2 (sY 3) = s 6 3 8
(s 3)2
;
6s 26 42s 6s2 80
) 4X 2sY = 6+ = : (4.13)
(s 3)2 (s 3)2
s2 3s X + 2sY = 8: (4.14)
Problem 4.7.7
dx
+ 3x 2y = 3 sin 2t 6e t ;
dt
dy
2x + = 3e t ;
dt
x (0) = 0 , y (0) = 4:
)X= 2
s2 +4
) x (t) = L 1
(X) = sin 2t .
dx
As dt
+ 3x 2y = 3 sin 2t 6e t ; then
dx t
2y = dt
+ 3x 3 sin 2t + 6e = 2 cos 2t + 3 sin 2t 3 sin 2t + 6e t :
) y (t) = cos 2t + 3e t
.
CHAPTER 4. THE LAPLACE TRANSFORM 131
dy dx
Problem 4.7.8 x + dt
= 2et ; dt
y=0;
x (0) = 0; y (0) = 1 :
2 s 3
) X + sY = 1 = : (4.19)
s 1 s 1
Also, sX Y = 0. We multiply by s, and get
s2 X sY = 0: (4.20)
As dx
dt
y = 0; then y = dx
dt
) y (t) = et + 2 cos t sin t .
dx dy
+ y = 0; x+ = 2 cos t:
dt dt
(Answer: x = cos t, y = sin t.)
dx dx dy
2y = 6e3t sin t; 10x 3 + 2 = 0:
dt dt dt
(Answer: x = 2e3t sin t, y = e3t cos t.)
CHAPTER 4. THE LAPLACE TRANSFORM 132
2s2 1
Problem 4.8.2 Let F (s) = (s2 +4)2
. Find L (F ).
) L (g) = e bs
L (k (t + b)) e cs
L (k (t + c)) : (4.24)
Problem 4.8.3 Draw a neat sketch for the following function, and compute
its Laplace transform:
8
< 0; 1<t<1
g (t) = 4t3 22t2 + 38t 10; 1 t<3 :
:
0; 3 t<1
Problem 4.8.4 Draw a neat sketch for the following function, and …nd its
Laplace transform (b 8is a positive constant):
< 0 1<t<0
gb (t) = 1=b 0 t<b .
:
0 b t<1
Remark 4.8.1 In the preceding problem, the function g (t) is a positive pulse
in the shape of an inverted parabola, centred at t = 2, with width 2c and height
3
. (See its graph above, when c = 1=8.) It is easy to see (by integration)
4c
that the area under its curve equals 1. In the limit case, when c tends to
0, this pulse gets very narrow and very high, while its centre remains the
point t = 2 and its area remains 1. This means that g (t) tends to the Dirac
function (t 2) as c tends to 0. Meanwhile, its Laplace transform tends
that of (t 2), as we now show:
3 cs cosh (cs) sinh (cs)
lim L (g) (s) = 3 e 2s lim
c !0 s c !0 c3
2
3 s cosh (cs) + cs sinh (cs) s cosh (cs)
= 3 e 2s lim (by L’Hôpital’s Rule)
s c !0 3c2
3 2s s2 sinh (cs)
= 3 e lim
s c !0 3c
3 2s s3 cosh (cs)
= 3 e lim (by L’Hôpital’s Rule).
s c !0 3
) lim L (g) (s) = e 2s lim cosh (cs) = e 2s ;
c !0 c !0
which is the Laplace transform of the Dirac function (t 2).
When computations are performed within many engineearing applications,
and when the width 2c of this pulse g is narrow enough, it is often reasonable
and acceptable to replace g by the Dirac function (t 2), and to replace the
rather cumbersome formula (4.25), for the Laplace transform of g, by the
much simpler L ( (t 2)) = e 2s .
Similar conclusions hold for the function described in Problem 4.8.4.
Appendix A
R 1
R p R
xr dx = xr+1 (r 6= 1) p dx = ln x + x2 a2 sec2 x dx = tan x
r+1 x2 a 2
R 1
R p R
dx = ln (Cx) p dx = ln x + x 2 + a2 csc2 x dx = cot x
x x2 +a2
R R R
ex dx = ex + C sin x dx = cos x + C sec x tan x dx = sec x
R dx 1 1 x
R R
x2 +a2
= a
tan a
+C cos x dx = sin x + C csc x cot x dx = csc x
R dx 1 x a
R R
x2 a 2
= 2a
ln x+a
+C sec x dx = ln jsec x + tan xj tan x dx = ln jsec xj
R 1
R R
p dx = sin x
+C csc x dx = ln jcsc x + cot xj cot x dx = ln jsin xj
a 2 x2 a
Trigonometric Formulas
sin x 1 2 tan x
tan x = cos x
sec x = cos x
tan 2x = 1 tan2 x
cos x 1 1 cot2 x 1
cot x = sin x
= tan x
csc x = sin x
cot 2x = 2 cot x
d
dx
sin x = cos x sin2 x + cos2 x = 1 sin x cos x = 21 sin 2x
d
dx
tan x = sec2 x sec2 x = 1 + tan2 x sin2 x = 12 (1 cos 2x)
d
dx
sec x = sec x tan x csc2 x = 1 + cot2 x cos2 x = 21 (1 + cos 2x)
139
Appendix B
Assignments
140
APPENDIX B. ASSIGNMENTS 141
Specialization: Date: / 20
Specialization: Date: / 20
Specialization: Date: / 20
dy
Assignment Problem B.3.1 dx
+ x2 y = (x4 sec2 x) y 3 :
APPENDIX B. ASSIGNMENTS 152
Specialization: Date: / 20
Specialization: Date: / 20
Specialization: Date: / 20
Assignment Problem B.6.1 Find the general solution of the following second-
order linear di¤erential equation, using the Method of Undetermined
Coe¢ cients:
y 00 + 4y 0 + 5y = 4ex 10xex :
APPENDIX B. ASSIGNMENTS 167
Assignment Problem B.6.2 Find the general solution of the following second-
order linear di¤erential equation, using the Method of Undetermined
Coe¢ cients:
y 00 2y 0 8y = 6e 2x :
APPENDIX B. ASSIGNMENTS 168
Assignment Problem B.6.3 Find the general solution of the following second-
order linear di¤erential equation, using the Method of Undetermined
Coe¢ cients:
y 00 + 6y 0 + 45y = (37x 35) e 2x :
APPENDIX B. ASSIGNMENTS 169
Assignment Problem B.6.4 Find the general solution of the following second-
order linear di¤erential equation, using the Method of Undetermined
Coe¢ cients:
y 00 3y 0 18y = 36 sin 3x 108 cos 3x:
APPENDIX B. ASSIGNMENTS 170
Assignment Problem B.6.5 Find the general solution of the following second-
order linear di¤erential equation, using the Method of Undetermined
Coe¢ cients:
y 00 2y 0 + 5y = 15x2 12x 29:
APPENDIX B. ASSIGNMENTS 171
Specialization: Date: / 20
Specialization: Date: / 20
5. p (t) = (t 2) e6(t 2)
cos 3 (t 2) U (t 2)
d
6. h (t) = dt
(te6t cos 3t)
Rt
7. m (t) = 0
z e6z cos 3z dz
1
3. L (w) = s(s 6)(s2 +9)
:
s2 12s+27
4. L (k) = (s2 12s+45)2
:
s2 12s+27
5. L (p) = (s2 12s+45)2
e 2s :
s(s2 12s+27)
6. L (h) = (s2 12s+45)2
:
s2 12s+27
7. L (m) = s(s2 12s+45)2
:
s 4 s+4 s2
8. L (v) = 2((s 4)2 +9)
+ 2((s+4)2 +9)
+ (s2 16)(s2 +9)
:
APPENDIX B. ASSIGNMENTS 178
Specialization: Date: / 20
1. F (s) = 7e 3s
+ 3s e 4s
+ 10
s2 25
e 6s 12
(s+6)4
e 8s :
1 s
2. G (s) = cot 4
:
s 6
3. K (s) = s2 12s+40
e 7s :
APPENDIX B. ASSIGNMENTS 180
3. K (s) = 3e 2s
+ 4s e 5s
+ s
s2 4
e 3s :
APPENDIX B. ASSIGNMENTS 181
Specialization: Date: / 20
Specialization: Date: / 20
Specialization: Date: / 20
dx
3x + 2y = 0;
dt
dx dy
3 13x + 2 = 0;
dt dt
x (0) = 1 , y (0) = 0:
APPENDIX B. ASSIGNMENTS 193
dx
4x + 2y = 2te4t ;
dt
dy
x+3 12y = 0:
dt
x (0) = 3 , y (0) = 0 .
APPENDIX B. ASSIGNMENTS 194
dy 2t
8x = 8te ;
dt
dx 2t
2x y = 2te ;
dt
x (0) = 0 , y (0) = 0:
APPENDIX B. ASSIGNMENTS 195
dx dy
+ 4 + 4y = 3e t ;
dt dt
dy
x = te t ;
dt
x (0) = 1 , y (0) = 0 .
APPENDIX B. ASSIGNMENTS 196