Differential Equations: Learning Objectives
Differential Equations: Learning Objectives
15
only see what remains to be done.”
….Buddha
CHAPTER
Differential Equations
Learning Objectives
After reading this chapter, you will know:
1. Order of Differential Equation , First Order Equations (Linear and Nonlinear)
2. Higher Order Linear Differential Equations with Constant Coefficients
3. Cauchys-Euler's Equations, Legendre's Linear Equation
4. Partial Differential Equation
Example:
i) ex dx + ey dy = 0 → 1st order, 1st degree
dy x
ii) y = xdx + dy/dx → 1st order, 2nd degree (power of highest derivative)
3/2
dy 2 d2 y
iii) [1 + (dx) ] = C dx2 → 2nd order, 2nd degree
d2 x
iv) + n2 x = 5 → 2nd order, 1st degree
dt2
d2 y dy
Example: The differential equation dx2
+ sin x dx
+ yex = sin hx is
(A) 1st order and linear (C) 2nd order and linear
(B) 1st order and non-linear (D) 2nd order and non-linear
Solution: Option (C)
1. Variable Separable
f(x)dx + g(y)dy = 0
∫ f(x)dx + ∫ g(y)dy = c is the solution
2. Homogenous Equation
dy f(x, y)
=
dx g(x, y)
To solve a homogeneous equation, substitute y = Vx
dy dv
= V + x
dx dx
Separate the variable V and x and integrate.
dy (2x + 3y) + 4
𝐄𝐱𝐚𝐦𝐩𝐥𝐞: =
dx 2(2x + 3y) + 5
a b
𝐒𝐨𝐥𝐮𝐭𝐢𝐨𝐧: Here = , put 2x + 3y = t
a1 b1
dy dt dy 1 dt
2 + 3 = ⇒ = ( − 2)
dx dx dx 3 dx
1 dt t+4
( − 2) =
3 dx 2t + 5
dt 3(t + 4) 7t + 22
=2+ =
dx 2t + 5 2t + 5
2t + 5
dt = dx
7t + 22
3. Linear Equations
The standard form of a linear equation of first order:
dy
+ P(x) y = Q(x) , where P and Q are functions of x
dx
d2 y dy
Second order linear eqn: 2
+ P(x) + Q(x)y = R(x)
dx dx
Commonly known as “Leibnitz’s linear equations”
Integrating factor, I.F. = e∫ Pdx
dy
e∫ Pdx + y(e∫ Pdx P) = Q e∫ Pdx
dx
d
(y e∫ Pdx ) = Q e∫ Pdx
dx
ye∫ Pdx = ∫ Q (I. F)dx + C ⇒ y(I. F) = ∫ Q(I. F)dx + C
Bernoulli’s Equation
dy⁄dx + P y = Q y n where, P & Q are functions of x only .
Divide by y n
y −n dy⁄dx + Py1−n = Q
Substitute, y1−n = z
(1 − n)y −n dy⁄dx = dz⁄dx
y −n dy/dx = 1/(1 − n) dz⁄dx
1
dz⁄dx + Pz = Q
1−n
dz⁄dx + (1 − n)Pz = Q (1 − n) → This is linear equation and can be solved easily
z x −5 = ∫ x −5 (−5x 2 ) dx
z x −5 = −5 ∫ x −3 dx
−5x −2
y −5 x −5 = +C
−2
5
( + Cx 2 ) x 3 y 5 = 1
2
1
Rule 1: when M dx + N dy = 0 is homogenous in x and y and M x + N y ≠ 0 then I.F. =
Mx+Ny
1
Rule 2: If the equation f1 (xy)ydx + f2 (xy)xdy = 0 and M x – N y ≠ 0 then I.F. =
Mx−Ny
1 ∂M ∂N
Rule 3: If the M dx + N dy = 0 and (
N
− ) = f(x), then I.F. = e∫ f(x)dx
∂y ∂x
1 ∂N ∂M
Rule 4: If the equation M dx + N dy = 0 and M
( − ) = f(y), then I.F. = e∫ f(y)dy
∂x ∂y
dy dx x y
Example: Solve the differential equation − = −
dx dy y x
Solution: Let P = dy/dx
1 x y
P− = −
P y x
y x
P 2 + P ( − ) −1 = 0
x y
y x
(P + ) (P − ) = 0
x y
y x
P+ =0 P− =0
x y
d(xy) = 0 y dy − x dx = 0
xy = C x2 − y2 = C
(x y – C) (x − y 2 − C) = 0
2
Example: The family of straight lines passing through origin is represented by differential equation
(A) ydx + xdy = 0 (C) xdx + ydy = 0
(B) xdy − ydx = 0 (D) ydy − xdx = 0
Solution: [Ans. B]
(A) d(yx) = 0 (C) y 2 − x 2 = c
yx = c
(B) d(y /x) = 0 (D) x 2 + y 2 = c
y
= c → y = cx
x
Case IV: If two or more than two pairs of root are imaginary
i.e., Repeated imaginary root α ± iβ, α ± iβ
y = eαx [(C1 x + C2 ) cos βx + (C3 x + C4 ) sin βx ]
Illustration: D4 + 4 = 0
(D2 + 2D + 2) (D2 − 2D + 2) = 0
D = −1 ± i, 1 ± i
y = e−x (C1 cos x + C2 sin x) + ex (C3 cos x + C4 sin x )
Complete Solution
dn y k1 (dn−1 y)
Given differential equation: + + . . . . . . . kn = X
dx n dx n−1
The complete solution is: y = C.F + P.
d2 y dy
Example: Particular integral of dx2
+ dx = x 2 + 2x + 4
x2 x3
(A) 3
+ 4x (C) 3
+ 4x
x3 x3
(B) 3
+4 (D) 3
+ 4x 2
1
Solution: D(D+1)
{ x 2 + 2x + 4}
1
= (1−D + D2 - - - -) { x 2 + 2x + 4}
D
1
= D[x 2 + 2x + 4 − 2x − 2 + 2]
1 x3
= D { x 2 + 4 } ⇒ 3 + 4x , Option (C)
C.F. = C1 ex + C2 e3x
1
P.I. = (D2 −4D+3)
(2xe3x + 3ex cos 2x)
1 1
= 2e3x (D+3)2 −4(D+3)+3 x + 3ex (D+1)2 −4(D+1)+3 cos 2x
1 1
= 2e3x D2 +2D x + 3ex D2 −2D cos 2x
1 1
= 2e3x x
D x +3e −4−2D cos 2x
2D(1+ )
2
1 D 1 3
= e3x D (1+ 2 )−1 x − 2 ex 2+D
cos 2x
1 D D 2 3 2−D
= e3x D (1 − 2 + 4 ) x − 2 ex 4−D2 cos 2x
1 1 D 3 (2−D)
= e3x (D − 2 + 4 )x − 2 ex 4+4 cos 2x
x2 1 1 3
= e3x [ 2 − 2 x + 4 ] − 16 ex (2cos 2x + 2 sin 2x)
x2 x 1 3
C.S. = C.F. + P.I. = C1 ex + C2 e3x + e3x ( 2 − 2 + 4 ) − 8 e3x (cos 2x + sin 2x)
e3x
Note: is omitted from P.I. since C2 e3x is already in C.F.
4
In the above problem, try to find P.I. using different approach,
1 1 x2 x
P.I. = e3x D (x – 2
+ 0) = e3x ( 2 − 2) → Note the difference, (overall C.S. is unchanged)
d2 y dy
Example: Find the complete solution of dx2
+ dx = x 2 + 2x + 4
Solution: D (D+1) = 0
C.F. = C1 e0.x + C2 e−x = C1 + C2 e−x
Other approach
1
P. I. = (1− D +D2 − D3 ) (x 2 + 2x + 4 )
D
1
= (D − 1 + D − D2 ) (x 2 + 2x + 4 )
3
x x3
+ x 2 + 4x) − (x 2 + 2x + 4 ) + (2x + 2) − 2 =
=( + 4x – 4
3 3
Note: This is because a constant C1 is present in C.F. (overall C.S is unchanged) C.S. = C1 + C2 e−x +
x3
3
+ 4x
Cauchy’s Euler Equation: (Homogenous Linear Equation)
dn y dn−1 y dy
xn n + k1 x n−1 n−1 + - - - - - - - k n−1 + kny = x
dx dx dx
d
Put x = et ⟹ t = ln x, then if D = dt
dy dy dt 1 dy
= = Dy ⇒x = Dy
dx dt dx x dx
d2 y dy 1 dy dy 1 d dt
= ( ) = − x12 +
dy
( dx )
1 1
= − x2 Dy + x D2 y
1
dx2 dx x dt dx x dt dx x
d2 y 1
= (D2 − D)y
dx2 x2
2
d y
x 2 2 = D (D − 1) y
dx
3
3 d y
x = D (D − 1)(D − 2) y
dx3
After substituting these differentials, the Cauchy – Euler equation results in a linear equation with
constant coefficients.
d2 y dy
Example: Solve the given differential equation x 2 −x dx + y = ln x
dx2
Solution: Substitute x = et ⇒ t = ln x
dy d2 y
x = Dy x 2 = D(D − 1)
dx dx2
2
(D − D − D + 1) y = t
(D − 1)2 = 0, y = (C1 +C2 t) et
1
P. = t = [1 − D]−2 t = t + 2
(D−1)2
y = (C1 +C2 ln x) x + ln x + 2
d2 y dy
Example: (2x + 3)2 dx2
− (2x + 3) dx − 20y = 0
t
Solution: 2x +3 = e , t = ln (2x+3)
4D (D −1)y – 2Dy – 20y = 0
(4D2 − 6D − 20)y = 0
2D2 − 3D − 10 = 0
D = 3.1, −1 .6
y = C1 e3.1t + C2 e−1.6 t = C1 (2x +5)3.1 + C2 (2x +5)−1.6
Partial Differential Equation
z = f(x, y)
∂z ∂z ∂2 z ∂2 z ∂2 z
= p, = q , 2 = r, = s, 2 = t
∂x ∂y ∂x ∂x ∂y ∂y
Working Rule
dx dy dz
1. Form subsidiary equation = =
P Q R
2. Solve these simultaneous equations giving u = a and v = b as its solution.
3. Write the complete solution as f (u ,v) = 0 , u = f(v)
Working Rule:
dz dz
1. Assume, u = x + ay , substitute p = , and q = a
du du
2. Solve resulting ordinary differential equation (O.D.E.) in z and u
3. Replace u by x + ay
Form III: f (x, p) = F (y, q). Terms not containing z and terms containing x and p can be separated
from those containing y and q
Assume f (x, p) = F (y, q) = a
p = ф(x), q = Ψ(y)
∂z ∂z
dz = dx + dy = pdx + qdy
∂x ∂y
dz = ф(x) dx + Ψ(y) dy
z = ∫ ф(x)dx + ∫ Ψ(y)dy + b
Charpit’s Method: It is a general method for finding the complete integral of a nonlinear PDF of first
order of the form
f(x, y, z, p, q ) = 0 ----------------------- ①
Since z depends on x and y
∂z ∂z
dz = ∂x dx + ∂y dy = pdx + qdy -------------- ②
dx dy dz dp dq dф
= = = = =
∂f ∂f ∂f ∂f ∂f ∂f ∂f ∂f 0
− − −p −q +p +q
∂p ∂q ∂p ∂q ∂x ∂z ∂y ∂z
m2 − 4m + 4 = 0 ⇒ m = 2, 2
C.F = f1 (y +2x) + x f2 (y +2x)
1
P.I. = (D−2D′ )2
e2x+y
D = 2, D = 1, ⇒ D – 2D′ = 0
1
P.I. = D−2D′
∫ e(2x+c−2x) dx
1
= D−2D′
xey+2x
c−2x+2x
x 2 ey+2x
= ∫ xe dx =
2
1
Hence, C.S. is, Z = f1 (y +2x) + x f2 (y +2x) + 2 x 2 e2x+y
Example: Solve the given differential equation y (2xy + ex ) dx = ex dy
Solution: (y ex dx − exdy ) + 2x y 2 dx = 0
yex dx − ex dy
+ 2x dx = 0
y2
ex
d ( ) + 2x dx = 0
y
e x ⁄y + x 2 = 0
2 ∂2 z 5 ∂2 z 2 ∂2 z
Example: Solve the given differential equation ∂x2
+ ∂x ∂y
+ ∂y2
=0
Solution: (2D2 + 5DD′ + 2D′2 ) z = 0
2m2 + 5m + 2 = 0
−5±√25−16 −5±3 1
m= = = −2, −
2×2 4 2
1
Solution is z = f1 (y −2x) + f2 (y − x)
2
∂2 z ∂2 z
Example: Solve the given differential equation ∂x2
− ∂x ∂y
= cos x cos 2y
Solution: (D2 − DD′ ) z = cos x cos 2y
m2 − m = 0; ⇒ m = 0, 1
C.F. = f1 (y) + f2 (y + x)
1
P.I. = D2 −DD′ cos x cos 2y
1 1
=2 D2 −DD′
[cos(x + 2y) + cos(x − 2y)]
1 1 1
= [ cos(x + 2y) + 2 cos(x − 2y)]
2 D2 −DD′ D −DD′
1 cos(x+2y) 1
= [ −1+2 + −1−(+2) cos(x − 2y)]
2
1 1
= 2
cos(x + 2y) − 6 cos(x − 2y)
∂3 z 2 ∂3 z
Example: Solve the given differential equation − = 2e2x + 3x 2 y
∂x3 ∂x2 dy
Solution: D3 − 2D2 D′ = 2e2x + 3x 2 y
m3 − 2m2 = 0
m2 (m − 2) = 0, ⇒ m = 0, 0, 2
C.F. = f1 (y) + x f2 (y) + f3 (y + 2x)
1
P.I. = (2e2x + 3x 2 y)
D3 −2D2 D′
1 1
= 2 e2x + 3 1−2D′
x2y
D3 −2D2 D′ D3 ( )
D
2 e2x 3 −1
1−2D′
= + ( ) x2y
23 −2 22 0 D3 D
2 e2x 3 2D′ (−1)(−1−1) 4D′
2
= + (1 + + ) x2y
8 D3 D 22 D2
1 2x 3 2x3
= 4
e + D3 (x 2 y + 3
+ 0)
2x 5
e x y 2 x6
= + 3( + )
4 60 3 4 × 5 × 6
e2x x 5 y x 6
= + +
4 20 60
Example: Solve the given differential equation
∂2 z ∂2 z 6 ∂2 z
+ − = cos (2x + y)
∂x 2 ∂x dy ∂y 2
Solution: (D2 + DD′ − 6D′2 )z = cos(2x +y)
m2 + m − 6 = 0, ⇒m = −3, 2
C.F. = f1 (y−3x) + f2 (y+2x)
cos(2x + y) cos(2x + y) cos(2x + y)
P. I. = 2 ′ ′2
= 2 2
=
D + DD − 6D −2 + (−2 × 1) − 6(−1 ) 0
1 1
= cos(2x + y) = cos(2x + y)
D2 + DD′ − 6D′2 (D + 3D′ )(D − 2D′ )
1
= ∫ cos(2x + ̅̅̅̅̅̅̅̅
c − 2x) dx
D + 3D′
1
= ∫ cos c dx
D + 3D′
1
= x cos c
D + 3D′
1
= x cos(2x + y)
D + 3D′
= ∫ x cos(2x + ̅̅̅̅̅̅̅̅
c + 3x) dx
= ∫ x cos(5x + c) dx
sin (5x + c) cos (5x + c)
= x +
5 25
x cos (5x + ̅̅̅̅̅̅̅̅
y − 3x)
= sin (5x + ̅̅̅̅̅̅̅̅
y − 3x) +
5 25
x 1
= sin (2x + y) + cos(2x + y)
5 25
Approach 2
1
P. I = cos(2x + y)
(D −2D′ )(D + 3D′ )
1
= ∫ cos(2x + ̅̅̅̅̅̅̅̅
c + 3x) dx
(D − 2D′ )
1
= ∫ cos(c + 5x) dx
(D − 2D′ )
1 sin(c + 5x)
= ′
(D − 2D ) 5
1 sin(y − 3x + 5x)
=
(D − 2D′ ) 5
1 1 1
= ′
sin(y + 2x) = ∫ sin(c̅̅̅̅̅̅̅̅
− 2x + 2x) dx
(D − 2D ) 5 5
1 x x
= ∫ sin c dx = cos c = cos(y + 2x)
5 5 5
Approach 3
1
P. I. = cos (2x + y)
(D − 2D′ )(D + 3D′ )
1
= cos (2x + y)
D2
+ DD′ − 6D′2
1
= x cos (2x + y)
2D + D′
D
= x cos (2x + y)
2D + DD′
2
D cos(2x + y) x
= x = D cos(2x + y)
2(−4) − 2 −10
x
= sin(2x + y) × 2
−10
x
= sin(2x + y)
5
Note: The difference in answer of approach 1 & approach 2/3 (why?)
To find order and degree of differential equation simplify all the decimal and fractional power to
3/2 or 1/2 to integer.
3/2
dy 2 ∂2 y
𝐄𝐱𝐚𝐦𝐩𝐥𝐞: [1 + ( ) ]
dx ∂x 2
3/2 2
dy 2 2
∂2 y
Square both side [1 + ( ) ] C ( 2)
dx ∂x