First Order and Higher Order Differential Equations
First Order and Higher Order Differential Equations
Differential Equation:
Ordinary differential equations: Differential equations which involve only one independent variable
and differential co-efficients w.r.t it are called Ordinary differential equations.
Partial differential equations: Differential equations which involve two or more independent variable
and partial derivatives w.r.t it are called Partial differential equations.
Exact Differential Equation: A differential equation is called exact differential equation if it is obtained
only by differentiation of its primitive without applying any other operation like addition, subtraction,
elimination etc.
e.g. 𝑥𝑑𝑦 + 𝑦𝑑𝑥 = 0
Integrating Factor: A suitable factor by multiplying with which, the given equation becomes exact is
called integrating factor.
Rules to find Integrating Factor:
1. By Inspection:
i. 𝑥𝑑𝑦 + 𝑦𝑑𝑥 = 𝑑(𝑥𝑦)
𝑦𝑑𝑥−𝑥𝑑𝑦 𝑥
ii. 𝑦2
= 𝑑 (𝑦)
𝑦𝑑𝑥−𝑥𝑑𝑦 𝑥
iii. x2 +𝑦 2
= 𝑑 (tan−1 𝑦)
𝑦𝑑𝑥−𝑥𝑑𝑦 𝑥
iv. x2 +𝑦 2
= 𝑑 (log 𝑦)
𝑦𝑑𝑥−𝑥𝑑𝑦
v. xy
= 𝑑(log xy)
dy
e x ( A cos x B sin x) e x ( A sin x B cos x) ………………(2)
dx
d2y
2
e x ( A cos x B sin x) e x ( A sin x B cos x) e x ( A sin x B cos x) e x ( A cos x B sin x)
dx
d 2 y dy
e x ( A sin x B cos x) e x ( A cos x B sin x) ( fromequation(2))
dx 2 dx
d 2 y dy dy
2
( y) y
dx dx dx
d2y dy
2
2 2y 0
dx dx
Assignment:
(1) y cx c
2
x
(2) y Ae Be C
x
x
(3) xy Ae Be x2
x
(4) Find the differential equation of all parabolas whose axes are parallel to y-axis.
(5) Form the differential equation of all circles of radius a.
g ( y )dy f ( x)dx
f ( x)dx g ( y)dy
dy
Q2. Solve x cos x cos y sin y 0
dx
dy
Sol. x cos x cos y sin y 0
dx
dy
x cos x cos y sin y
dx
dy
x cos x tan y
dx
Homogeneous equation:-
If a differential equation of the first order and first degree
dy y
f( ) …………….(1)
dx x
Putting y =vx
dy dv
vx
dx dx
From equation(1)
dv vx
vx f( )
dx x
dv
vx f (v )
dx
dv
f (v ) v x
dx
dv dx
( f (v ) v ) x
dv dx
( f (v ) v ) x C
Q.3 Solve xdy ydx x 2 y 2 dx
dy y x y
2 2
Sol. ……………(1)
dx x
Putting y =vx
dy dv
vx
dx dx
From equation(1)
dv vx x 2 v 2 x 2
vx
dx x
dv
vx v 1 v2
dx
dv
x 1 v2
dx
dv dx
1 v2 x
dv dx
1 v2
x
log( v 1 v 2 ) log( x) log c
v 1 v 2 cx
y x 2 y 2 cx 2
Assignment:-
y y y
(1) ( x tan y sec 2 )dx x sec 2 dy 0
x x x
(2) ( x y )dx ( x y )dy 0
dy y y
(3) sin
dx x x
x x
x
(4) (1 e )dx e (1 )dy 0
y y
y
x x
(5) ye dx ( xe y )dy
y y
a b
When '
'
a b
dY a( X h) b(Y k ) c
'
dX a ( X h) b ' (Y k ) c '
dY aX bY ah bk c
' ………….(1)
dX a X b 'Y a ' h b ' k c '
ah bk c 0
a 'h b'k c' 0
h k 1
' '
bc b c ca c a ab a ' b
' ' '
dY aX bY
' which is hom ogenous in X , Y and can be solved by
dX a X b 'Y
putting Y vX
State and Prove Necessary and sufficient condition for differential eq. Mdx+Ndy =0 to be exact.
Solve (x2-ay)dx=(ax-y2)dy
Sol. Here 𝑀 = 𝑥 2 − 𝑎𝑦 and 𝑁 = 𝑦 2 − 𝑎𝑥
𝜕𝑀 𝜕𝑁
𝜕𝑦
= −𝑎 and 𝜕𝑥 = −𝑎
𝜕𝑀 𝜕𝑁
Therefore = Hence equation is exact.
𝜕𝑦 𝜕𝑥
Solution is given by ∫ 𝑀𝑑𝑥 + ∫ 𝑡𝑒𝑟𝑚𝑠 𝑜𝑓 𝑁 𝑓𝑟𝑒𝑒 𝑓𝑟𝑜𝑚 𝑥 𝑑𝑦 = 𝑐
𝑦𝑖𝑒𝑙𝑑𝑠 𝑦𝑖𝑒𝑙𝑑𝑠 𝑥 3 𝑦3
→ ∫(𝑥 2 − 𝑎𝑦)𝑑𝑥 + ∫(𝑦 2 )𝑑𝑦 = 𝑐 → − 𝑎𝑥𝑦 + =𝑐
3 3
Solve (x2-4xy-2y2)dx+(y2-4xy-2x2)dy=0
Sol. Here 𝑀 = 𝑥 2 − 4𝑥𝑦 − 2𝑦 2 and 𝑁 = 𝑦 2 − 4𝑥𝑦 − 2𝑥 2
𝜕𝑀 𝜕𝑁
𝜕𝑦
= −4𝑥 − 4𝑦 and 𝜕𝑥 = −4𝑦 − 4𝑥
𝜕𝑀 𝜕𝑁
Therefore 𝜕𝑦
= 𝜕𝑥
Hence equation is exact.
Solution is given by ∫ 𝑀𝑑𝑥 + ∫ 𝑡𝑒𝑟𝑚𝑠 𝑜𝑓 𝑁 𝑓𝑟𝑒𝑒 𝑓𝑟𝑜𝑚 𝑥 𝑑𝑦 = 𝑐
𝑦𝑖𝑒𝑙𝑑𝑠 𝑦𝑖𝑒𝑙𝑑𝑠 𝑥 3 𝑦3
→ ∫(𝑥 2 − 4𝑥𝑦 − 2𝑦 2 )𝑑𝑥 + ∫(𝑦 2 )𝑑𝑦 = 𝑐 → − 2𝑥 2 𝑦 − 2𝑦 2 𝑥 + =𝑐
3 3
Solve (3x2+6xy2)dx+(6x2y+4y3)dy=0
Sol. Here 𝑀 = 3𝑥 2 + 6𝑥𝑦 2 and 𝑁 = 6𝑥 2 𝑦 + 4𝑦 3
𝜕𝑀 𝜕𝑁
𝜕𝑦
= 12𝑥𝑦 and 𝜕𝑥 = 12𝑥𝑦
𝜕𝑀 𝜕𝑁
Therefore 𝜕𝑦
= 𝜕𝑥
Hence equation is exact.
Solution is given by ∫ 𝑀𝑑𝑥 + ∫ 𝑡𝑒𝑟𝑚𝑠 𝑜𝑓 𝑁 𝑓𝑟𝑒𝑒 𝑓𝑟𝑜𝑚 𝑥 𝑑𝑦 = 𝑐
𝑦𝑖𝑒𝑙𝑑𝑠 𝑦𝑖𝑒𝑙𝑑𝑠
→ ∫(3𝑥 2 + 6𝑥𝑦 2 )𝑑𝑥 + ∫(4𝑦 3 )𝑑𝑦 = 𝑐 → 𝑥 3 + 3𝑥 2 𝑦 2 + 𝑦 4 = 𝑐
𝑦 𝑥
Sol. Here 𝑀 = 𝑦 − and 𝑁=𝑥+
𝑥 2 +𝑦 2 𝑥 2 +𝑦 2
𝜕𝑀 𝑦 2 −𝑥 2 𝜕𝑁 𝑦 2 −𝑥 2
𝜕𝑦
= 1 + (𝑥 2 +𝑦2 )2and 𝜕𝑥
=1 + (𝑥 2 +𝑦2 )2
𝜕𝑀 𝜕𝑁
Therefore 𝜕𝑦
= 𝜕𝑥
Hence equation is exact.
Solution is given by ∫ 𝑀𝑑𝑥 + ∫ 𝑡𝑒𝑟𝑚𝑠 𝑜𝑓 𝑁 𝑓𝑟𝑒𝑒 𝑓𝑟𝑜𝑚 𝑥 𝑑𝑦 = 𝑐
𝑦𝑖𝑒𝑙𝑑𝑠 𝑦 𝑦𝑖𝑒𝑙𝑑𝑠 𝑥
→ ∫ (𝑦 − 2 2
) 𝑑𝑥 + ∫(0)𝑑𝑦 = 𝑐 → 𝑦𝑥 − 𝑡𝑎𝑛−1 = 𝑐
𝑥 +𝑦 𝑦
x3
Solve ydx-xdy+3x2y2 e dx =0
𝑦𝑑𝑥−𝑥𝑑𝑦 3 𝑦𝑖𝑒𝑙𝑑𝑠 𝑥 3
Sol. Dividing both sides by 𝑦 2 , we get 𝑦2
+ 3𝑥 2 𝑒 𝑥 𝑑𝑥 = 0 → 𝑑 (𝑦) + 𝑑(𝑒 𝑥 ) = 0
𝑥 𝑥3
Integrating both sides, we get 𝑦
+𝑒 =𝑐
Solve (x2y-2xy2)dx-(x3-3x2y)dy=0
Sol. Here 𝑀 = 𝑥 2 𝑦 − 2𝑥𝑦 2 and 𝑁 = 3𝑥 2 𝑦 − 𝑥 3
𝜕𝑀 𝜕𝑁
𝜕𝑦
= 𝑥 2 − 4𝑥𝑦 and 𝜕𝑥
= 6𝑥𝑦 − 3𝑥 2
𝜕𝑀 𝜕𝑁
Therefore ≠ Hence equation is not exact.
𝜕𝑦 𝜕𝑥
1 1 1
Since the given equation is homogenous, so I.F. = 𝑀𝑥+𝑁𝑦 = 𝑥 3 𝑦−2𝑥2 𝑦2 +3𝑥2 𝑦2 −𝑥3 𝑦 = 𝑥 2 𝑦2
1 1 2 𝑥 3
Multiplying the given equation by , we get ( − ) 𝑑𝑥 − ( − ) 𝑑𝑦 = 0
𝑥 2𝑦2 𝑦 𝑥 𝑦2 𝑦
1 2 𝑥 3
Now 𝑀 = (𝑦 − 𝑥)and 𝑁= − (𝑦2 − 𝑦)
𝜕𝑀 −1 𝜕𝑁 −1
𝜕𝑦
=
𝑦2
and 𝜕𝑥
= 𝑦2
𝜕𝑀 𝜕𝑁
Therefore 𝜕𝑦 = 𝜕𝑥 Hence equation is exact.
Solution is given by ∫ 𝑀𝑑𝑥 + ∫ 𝑡𝑒𝑟𝑚𝑠 𝑜𝑓 𝑁 𝑓𝑟𝑒𝑒 𝑓𝑟𝑜𝑚 𝑥 𝑑𝑦 = 𝑐
𝑦𝑖𝑒𝑙𝑑𝑠 1 2 3 𝑦𝑖𝑒𝑙𝑑𝑠 𝑥
→ ∫ ( − ) 𝑑𝑥 + ∫ ( ) 𝑑𝑦 = 𝑐 → − 2 𝑙𝑜𝑔 𝑥 + 3 𝑙𝑜𝑔 𝑦= 𝑐
𝑦 𝑥 𝑦 𝑦
Solve xdy-ydx=(x2+y2)dx
Solve xdy-ydx=xy2 dx
1 1 1 1 1
Multiplying the given equation by2𝑥2 𝑦2 , we get (2𝑥 2 𝑦 + 2𝑥) 𝑑𝑥 + (2𝑥𝑦2 − 2𝑦) 𝑑𝑦 = 0
1 1 1 1
Now 𝑀 = (2𝑥 2 𝑦 + 2𝑥)and 𝑁 = (2𝑥𝑦2 − 2𝑦)
𝜕𝑀 −1 𝜕𝑁 −1
= and =
𝜕𝑦 2𝑥 2 𝑦 2 𝜕𝑥 2𝑥 2 𝑦 2
𝜕𝑀 𝜕𝑁
Therefore 𝜕𝑦 = 𝜕𝑥 Hence equation is exact.
Solution is given by ∫ 𝑀𝑑𝑥 + ∫ 𝑡𝑒𝑟𝑚𝑠 𝑜𝑓 𝑁 𝑓𝑟𝑒𝑒 𝑓𝑟𝑜𝑚 𝑥 𝑑𝑦 = 𝑐
𝑦𝑖𝑒𝑙𝑑𝑠 1 1 1 𝑦𝑖𝑒𝑙𝑑𝑠 −1 1 1
→ ∫ ( 2 + ) 𝑑𝑥 + ∫ (− ) 𝑑𝑦 = 𝑐 → + 𝑙𝑜𝑔 𝑥 − 𝑙𝑜𝑔 𝑦= 𝑐
2𝑥 𝑦 2𝑥 2𝑦 2𝑥𝑦 2 2
Solve x2ydx-(x3+y3)dy=0
Sol. Here 𝑀 = 𝑥 2 𝑦 and 𝑁 = −𝑥 3 − 𝑦 3
𝜕𝑀 𝜕𝑁
= 𝑥 2 and = −3𝑥 2
𝜕𝑦 𝜕𝑥
𝜕𝑀 𝜕𝑁
Therefore 𝜕𝑦
≠ 𝜕𝑥 Hence equation is not exact.
1 1 −1
Since the given equation is homogenous, so I.F. = 𝑀𝑥+𝑁𝑦 = 𝑥 3 𝑦−𝑥3 𝑦−𝑦4 = 𝑦4
−1 −𝑥 2 𝑥3 1
Multiplying the given equation by 𝑦4 , we get ( 𝑦3 ) 𝑑𝑥 + (𝑦4 + 𝑦) 𝑑𝑦 = 0
−𝑥 2 𝑥3 1
Now 𝑀 = ( 𝑦3 )and 𝑁 = (𝑦4 + 𝑦)
𝜕𝑀 3𝑥 2 𝜕𝑁 3𝑥 2
𝜕𝑦
=
𝑦4
and 𝜕𝑥
= 𝑦4
𝜕𝑀 𝜕𝑁
Therefore 𝜕𝑦 = 𝜕𝑥 Hence equation is exact.
Solution is given by ∫ 𝑀𝑑𝑥 + ∫ 𝑡𝑒𝑟𝑚𝑠 𝑜𝑓 𝑁 𝑓𝑟𝑒𝑒 𝑓𝑟𝑜𝑚 𝑥 𝑑𝑦 = 𝑐
𝑦𝑖𝑒𝑙𝑑𝑠 −𝑥 2 1 𝑦𝑖𝑒𝑙𝑑𝑠 𝑥 3
→ ∫ ( 3 ) 𝑑𝑥 + ∫ ( ) 𝑑𝑦 = 𝑐 → + 𝑙𝑜𝑔 𝑦= 𝑐
𝑦 𝑦 3𝑦 3
Solve (2x2y-3y4)dx+(3x3+2xy3)dy=0
𝜕𝑦
=13
𝑥 13 𝑦 13 −
13
𝑥 13 𝑦 13 and
𝜕𝑥
= 13
𝑥 13 𝑦 13 − 13 𝑥 13 𝑦 13
𝜕𝑀 𝜕𝑁
Therefore = Hence equation is exact.
𝜕𝑦 𝜕𝑥
Solution is given by ∫ 𝑀𝑑𝑥 + ∫ 𝑡𝑒𝑟𝑚𝑠 𝑜𝑓 𝑁 𝑓𝑟𝑒𝑒 𝑓𝑟𝑜𝑚 𝑥 𝑑𝑦 = 𝑐
𝑦𝑖𝑒𝑙𝑑𝑠 −23 −15 −49 24 𝑦𝑖𝑒𝑙𝑑𝑠 −13 −10 −15 13 −36 24
→ ∫ (2𝑥 13 𝑦 13 − 3𝑥 13 𝑦 13 ) 𝑑𝑥 + ∫ 0𝑑𝑦 = 𝑐 → 𝑥 13 𝑦 13 + 𝑥 13 𝑦 13 = 𝑐
5 12
2 3 𝑦 3 2 )𝑑𝑥 3 3 𝑦
Solve (3𝑥 𝑦 𝑒 +𝑦 +𝑦 + (𝑥 𝑦 𝑒 − 𝑥𝑦)𝑑𝑦 = 0
Sol.Here 𝑀 = (3𝑥 2 𝑦 3 𝑒 𝑦 + 𝑦 3 + 𝑦 2 )and 𝑁 = (𝑥 3 𝑦 3 𝑒 𝑦 − 𝑥𝑦)
𝜕𝑀 𝜕𝑁
= (9𝑥 2 𝑦 2 𝑒 𝑦 + 3𝑥 2 𝑦 3 𝑒 𝑦 + 3𝑦 2 + 2𝑦)and = 3𝑥 2 𝑦 3 𝑒 𝑦 − 𝑦
𝜕𝑦 𝜕𝑥
𝜕𝑀 𝜕𝑁
Therefore 𝜕𝑦
≠ 𝜕𝑥
Hence equation is not exact.
𝜕𝑁 𝜕𝑀
− (3𝑥 2 𝑦 3 𝑒 𝑦 −𝑦)−(9𝑥 2 𝑦 2 𝑒 𝑦 +3𝑥 2 𝑦 3 𝑒 𝑦 +3𝑦 2 +2𝑦) −3
𝜕𝑥 𝜕𝑦
Now 𝑀
= 3𝑥 2 𝑦 3 𝑒 𝑦 +𝑦 3 +𝑦 2
= 𝑦
= 𝑔(𝑦), a function of y only.
−3
∫ 𝑦 𝑑𝑦 1
∴ I. F. = 𝑒 ∫ 𝑔(𝑦)𝑑𝑦 = 𝑒 = 𝑒 −3 ln 𝑦 = 3
𝑦
1
Multiplying the given equation by𝑦3 , we get
1 𝑥
(3𝑥 2 𝑒 𝑦 + 1 + ) 𝑑𝑥 + (𝑥 3 𝑒 𝑦 − 2 ) 𝑑𝑦 = 0
𝑦 𝑦
1
Now 𝑀 = (3𝑥 2 𝑒 𝑦 + 1 + 𝑦)
𝑥
and 𝑁 = (𝑥 3 𝑒 𝑦 − )
𝑦2
𝜕𝑀 1 𝜕𝑁 1
𝜕𝑦
= (3𝑥 2 𝑒 𝑦 − 𝑦2 )and 𝜕𝑥
= (3𝑥 2 𝑒 𝑦 − 𝑦2 )
𝜕𝑀 𝜕𝑁
Therefore 𝜕𝑦 = 𝜕𝑥 Hence equation is exact.
Solution is given by ∫ 𝑀𝑑𝑥 + ∫ 𝑡𝑒𝑟𝑚𝑠 𝑜𝑓 𝑁 𝑓𝑟𝑒𝑒 𝑓𝑟𝑜𝑚 𝑥 𝑑𝑦 = 𝑐
𝑦𝑖𝑒𝑙𝑑𝑠 1 𝑦𝑖𝑒𝑙𝑑𝑠 𝑥
→ ∫ (3𝑥 2 𝑒 𝑦 + 1 + ) 𝑑𝑥 + ∫ 0𝑑𝑦 = 𝑐 → 𝑥3𝑒𝑦 + 𝑥 + = 𝑐
𝑦 𝑦
and 𝑁 = 6𝑥 2 𝑦
𝜕𝑀 𝜕𝑁
= 12𝑥𝑦and = 12𝑥𝑦
𝜕𝑦 𝜕𝑥
𝜕𝑀 𝜕𝑁
Therefore 𝜕𝑦
= 𝜕𝑥
Hence equation is exact.
Solution is given by ∫ 𝑀𝑑𝑥 + ∫ 𝑡𝑒𝑟𝑚𝑠 𝑜𝑓 𝑁 𝑓𝑟𝑒𝑒 𝑓𝑟𝑜𝑚 𝑥 𝑑𝑦 = 𝑐
𝑦𝑖𝑒𝑙𝑑𝑠 𝑦𝑖𝑒𝑙𝑑𝑠
→ ∫(5𝑥 4 + 12𝑥 3 + 6𝑥𝑦 2 )𝑑𝑥 + ∫ 0𝑑𝑦 = 𝑐 → 𝑥 5 + 3𝑥 4 + 3𝑥 2 𝑦 2 = 𝑐
−11 −19
Multiplying the given equation by 𝑥 7 𝑦 7 , we get
3 −5 −11 −12 −4 −19 10 −12
(2𝑥 𝑦 7 7 +𝑥 7 𝑦 7 ) 𝑑𝑥 + (3𝑥 𝑦 7 7 −𝑥7𝑦 7 ) 𝑑𝑦 = 0
3 −5 −11 −12
Now 𝑀 = (2𝑥 𝑦 7 7 +𝑥 7 𝑦 7 )
−4 −19 10 −12
and 𝑁 = (3𝑥 7 𝑦 7 −𝑥7𝑦 7 )
3 −12 −11 −19
𝜕𝑀 −10 12 𝜕𝑁 −10 3 −12 12 −11 −19 𝜕𝑀 𝜕𝑁
𝜕𝑦
= 7
𝑥 7 𝑦 7 − 7
𝑥 7 𝑦 7 and 𝜕𝑥
= 7
𝑥 7𝑦 7 − 7
𝑥 7 𝑦 7 Therefore 𝜕𝑦
= 𝜕𝑥
Hence equation is exact.
Solution is given by ∫ 𝑀𝑑𝑥 + ∫ 𝑡𝑒𝑟𝑚𝑠 𝑜𝑓 𝑁 𝑓𝑟𝑒𝑒 𝑓𝑟𝑜𝑚 𝑥 𝑑𝑦 = 𝑐
𝑦𝑖𝑒𝑙𝑑𝑠 3 −5 −11 −12
→ ∫ (2𝑥 7 𝑦 7 + 𝑥 7 𝑦 7 ) 𝑑𝑥 + ∫ 0𝑑𝑦 = 𝑐
𝑦𝑖𝑒𝑙𝑑𝑠 7 10 −5
7 −4 −12
→ 5
𝑥 𝑦 7 7 − 4𝑥 𝑦 7 7 =𝑐
1 1 2 1 1
Multiplying the given equation by 3𝑥 3 𝑦 3
, we get (3𝑥2 𝑦 + 3𝑥) 𝑑𝑥 + (3𝑥𝑦2 − 3𝑦) 𝑑𝑦 = 0
1 2 1 1
Now 𝑀 = ( + )and 𝑁 = ( 2− )
3𝑥 2 𝑦 3𝑥 3𝑥𝑦 3𝑦
𝜕𝑀 −1 𝜕𝑁 −1
= and = 2 2
𝜕𝑦 3𝑥 2 𝑦 2 𝜕𝑥 3𝑥 𝑦
𝜕𝑀 𝜕𝑁
Therefore 𝜕𝑦
= 𝜕𝑥
Hence equation is exact.
Solution is given by ∫ 𝑀𝑑𝑥 + ∫ 𝑡𝑒𝑟𝑚𝑠 𝑜𝑓 𝑁 𝑓𝑟𝑒𝑒 𝑓𝑟𝑜𝑚 𝑥 𝑑𝑦 = 𝑐
𝑦𝑖𝑒𝑙𝑑𝑠 1 2 1 𝑦𝑖𝑒𝑙𝑑𝑠 −1 2 1
→ ∫ ( 2 + ) 𝑑𝑥 + ∫ (− ) 𝑑𝑦 = 𝑐 → + 𝑙𝑜𝑔 𝑥 − 𝑙𝑜𝑔 𝑦= 𝑐
3𝑥 𝑦 3𝑥 3𝑦 3𝑥𝑦 3 3
Solve (3𝑥𝑦 2 − 𝑦 3 )𝑑𝑥 − (2𝑥 2 𝑦 − 𝑥𝑦 2 )𝑑𝑦 = 0
Sol. Here 𝑀 = 3𝑥𝑦 2 − 𝑦 3 and 𝑁 = −2𝑥 2 𝑦 + 𝑥𝑦 2
𝜕𝑀 𝜕𝑁
𝜕𝑦
= 6𝑥𝑦 − 3𝑦 2 and 𝜕𝑥
= −4𝑥𝑦 + 𝑦 2
𝜕𝑀 𝜕𝑁
Therefore ≠ Hence equation is not exact.
𝜕𝑦 𝜕𝑥
Since the given equation is homogenous, so
1 1 1
I.F. = = 2 2 3 2 2 3= 2 2
𝑀𝑥+𝑁𝑦 3𝑥 𝑦 −𝑥𝑦 −2𝑥 𝑦 +𝑥𝑦 𝑥 𝑦
1
Multiplying the given equation by 𝑥 2𝑦2
, we get
3 𝑦 2 1
(𝑥 − 𝑥 2 ) 𝑑𝑥 − (y − 𝑥) 𝑑𝑦 = 0
3 𝑦 −2 1
Now 𝑀 = (𝑥 − 𝑥 2 ) and 𝑁 = ( y + 𝑥)
𝜕𝑀 −1 𝜕𝑁 −1
𝜕𝑦 𝑦2
=and 𝜕𝑥
= 𝑦2
𝜕𝑀 𝜕𝑁
Therefore 𝜕𝑦 = 𝜕𝑥 Hence equation is exact.
Solution is given by ∫ 𝑀𝑑𝑥 + ∫ 𝑡𝑒𝑟𝑚𝑠 𝑜𝑓 𝑁 𝑓𝑟𝑒𝑒 𝑓𝑟𝑜𝑚 𝑥 𝑑𝑦 = 𝑐
𝑦𝑖𝑒𝑙𝑑𝑠 3 𝑦 −2 𝑦𝑖𝑒𝑙𝑑𝑠 𝑦
→ ∫ ( − 2 ) 𝑑𝑥 + ∫ ( ) 𝑑𝑦 = 𝑐 → 3 𝑙𝑜𝑔 𝑥 + + 3 − 2 𝑙𝑜𝑔 𝑦= 𝑐
𝑥 𝑥 y x
Solve (𝑥 2 𝑦 2 + xy + 1)𝑦𝑑𝑥 + (𝑥 2 𝑦 2 − xy + 1)𝑥𝑑𝑦 = 0
Sol. Here 𝑀 = (𝑥 2 𝑦 2 + xy + 1)𝑦 and 𝑁 = (𝑥 2 𝑦 2 − xy + 1)𝑥
𝜕𝑀 𝜕𝑁
𝜕𝑦
= 3𝑥 2 𝑦 2 + 2xy + 1 and 𝜕𝑥
= 3𝑥 2 𝑦 2 − 2xy + 1
𝜕𝑀 𝜕𝑁
Therefore 𝜕𝑦
≠ 𝜕𝑥
Hence equation is not exact.
1 1
Since the given equation is of the form (𝑥𝑦)𝑦𝑑𝑥 + 𝑔(𝑥𝑦)𝑥𝑑𝑦 = 0 , I.F. = =
𝑀𝑥−𝑁𝑦 2𝑥 2 𝑦 2
1 1 2 1 1
Multiplying the given equation by , we get ( + ) 𝑑𝑥 +( − ) 𝑑𝑦 =0
2𝑥 2 𝑦 2 3𝑥 2 𝑦 3𝑥 3𝑥𝑦 2 3𝑦
1 2 1 1
Now 𝑀 = (3𝑥 2 𝑦 + 3𝑥
)and 𝑁 = (3𝑥𝑦2 − 3𝑦)
𝜕𝑀 −1 𝜕𝑁 −1
𝜕𝑦
= 3𝑥 2 𝑦2 and 𝜕𝑥
= 3𝑥2 𝑦2
𝜕𝑀 𝜕𝑁
Therefore 𝜕𝑦
= 𝜕𝑥
Hence equation is exact.
Solution is given by ∫ 𝑀𝑑𝑥 + ∫ 𝑡𝑒𝑟𝑚𝑠 𝑜𝑓 𝑁 𝑓𝑟𝑒𝑒 𝑓𝑟𝑜𝑚 𝑥 𝑑𝑦 = 𝑐
𝑦𝑖𝑒𝑙𝑑𝑠 1 2 1 𝑦𝑖𝑒𝑙𝑑𝑠 −1 2 1
→ ∫ ( 2 + ) 𝑑𝑥 + ∫ (− ) 𝑑𝑦 = 𝑐 → + 𝑙𝑜𝑔 𝑥 − 𝑙𝑜𝑔 𝑦= 𝑐
3𝑥 𝑦 3𝑥 3𝑦 3𝑥𝑦 3 3
Solve 𝑦 + 𝑝𝑥 = 𝑥 4 𝑝2
𝑦𝑖𝑒𝑙𝑑𝑠
Sol. Given equation is 𝑦 + 𝑝𝑥 = 𝑥 4 𝑝2 → 𝑦 = 𝑥 4 𝑝2 − 𝑝𝑥 ……………(i)
𝑑𝑦 𝑑𝑝 𝑑𝑝 𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑝 𝑑𝑝
= 4𝑥 3 𝑝2 + 2𝑥 4 𝑝 −𝑝−𝑥 → 𝑝 = 4𝑥 3 𝑝2 + 2𝑥 4 𝑝 −𝑝−𝑥
𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑑𝑥
𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑝 𝑑𝑝 𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑝 𝑑𝑝
→ 4𝑥 3 𝑝2 + 2𝑥 4 𝑝 − 2𝑝 − 𝑥 =0 → 2𝑥 3 𝑝 (2𝑝 + 𝑥 ) − (2𝑝 + 𝑥 ) = 0
𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑑𝑥
𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑝 𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑝
→ (2𝑥 3 𝑝 − 1) (2𝑝 + 𝑥 )=0 → (2𝑝 + 𝑥 ) = 0
𝑑𝑥 𝑑𝑥
𝑑𝑝
(By neglecting the factor(2𝑥 3 𝑝 − 1) = 0 𝑏𝑒𝑐𝑎𝑢𝑠𝑒 𝑖𝑡 𝑑𝑜𝑒𝑠𝑛𝑜𝑡 𝑐𝑜𝑛𝑡𝑎𝑖𝑛 𝑑𝑥
)
𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑥 𝑑𝑝 𝑦𝑖𝑒𝑙𝑑𝑠 𝑐
→ 2 + =0→ 𝑝= 2
𝑥 𝑝 𝑥
𝑐2 𝑐
Putting the value of p in (i), we get 𝑦 = 𝑥 4 𝑥 4 − 𝑥 2 𝑥
𝑦𝑖𝑒𝑙𝑑𝑠 𝑐
→ 𝑦 = 𝑐 2 − 𝑥is the required solution
𝑑𝑦
Bernoulli’s form: An equation of the form + 𝑃(𝑥)𝑦 = 𝑦 𝑛 𝑄(𝑥) is called Bernoulli’s form.
𝑑𝑥
𝑑𝑦
Dividing both sides by 𝑦 𝑛 , we get 𝑦 −𝑛 𝑑𝑥 + 𝑃(𝑥)𝑦1−𝑛 = 𝑄(𝑥) … … … … … … . . (𝑖)
𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑦 𝑑𝑧 𝑦𝑖𝑒𝑙𝑑𝑠 −𝑛 𝑑𝑦 1 𝑑𝑧
Put 𝑦1−𝑛 = 𝑧 → (1 − 𝑛)𝑦 −𝑛 = → 𝑦 = (1−𝑛) .
𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑑𝑥
1 𝑑𝑧 𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑧
∴ (𝑖)𝑏𝑒𝑐𝑜𝑚𝑒𝑠 (1−𝑛) 𝑑𝑥
+ 𝑃(𝑥)𝑧 = 𝑄(𝑥) → 𝑑𝑥
+ (1 − 𝑛)𝑃(𝑥)𝑧 = (1 − 𝑛)𝑄(𝑥).
It is Leibnitz form.
1. If the roots of Auxiliary equation are real and distinct say 𝑚1 , 𝑚2 , 𝑚3 , then
𝐶. 𝐹. = 𝑐1 𝑒 𝑚1 𝑥 + 𝑐2 𝑒 𝑚2 𝑥 + 𝑐3 𝑒 𝑚3 𝑥
2. If the roots of Auxiliary equation are real and distinct say 𝑚, 𝑚, 𝑚, 𝑚, 𝑚 , then
𝐶. 𝐹. = 𝑐1 𝑒 𝑚𝑥 + 𝑥𝑐2 𝑒 𝑚𝑥 + 𝑥 2 𝑐3 𝑒 𝑚𝑥 + 𝑥 3 𝑐4 𝑒 𝑚𝑥 + 𝑥 4 𝑐5 𝑒 𝑚𝑥
3. If the roots are a pair of complex conjugate numbers say 𝛼 ± 𝑖𝛽, then
𝐶. 𝐹. = 𝑒 𝛼𝑥 (𝑐1 cos 𝛽𝑥 + 𝑐2 sin 𝛽𝑥)
4. If the roots are repeated pair of complex conjugate numbers say 𝛼 ± 𝑖𝛽, 𝛼 ± 𝑖𝛽, 𝛼 ± 𝑖𝛽
then
𝐶. 𝐹. = 𝑒 𝛼𝑥 [(𝑐1 + 𝑥𝑐2 + 𝑥 2 𝑐3 ) cos 𝛽𝑥 + (𝑐1 + 𝑥𝑐2 + 𝑥 2 𝑐3 ) sin 𝛽𝑥]
Rules for Particular integral: Let the equation be of the form 𝐹(𝐷)𝑦 = 𝑓(𝑥).
1
𝑃. 𝐼. = 𝑓(𝑥)
𝐹(𝐷)
1. If 𝑓(𝑥, 𝑦) = 𝑒 𝑎𝑥
1
𝑒 𝑎𝑥
𝑃. 𝐼. =
𝐹(𝐷)
Put 𝐷 = 𝑎 provided denominator does not become zero. If denominator becomes zero,
then
1 𝑎𝑥
𝑃. 𝐼. = 𝑥 𝑑 𝑒
𝐹(𝐷)
𝑑𝐷
Put 𝐷 = 𝑎 provided denominator does not become zero. If denominator becomes zero,
then
1
𝑃. 𝐼. = 𝑥 2 𝑑2 𝑒 𝑎𝑥
𝑑𝐷 2
𝐹(𝐷)
Put 𝐷 = 𝑎 provided denominator does not become zero.
Continue this process till we get a non-zero denominator.
2. If 𝑓(𝑥) = cos(𝑎𝑥) 𝑜𝑟 sin(𝑎𝑥)
1
𝑃. 𝐼. = 𝑓(𝑥)
𝐹(𝐷)
Put 𝐷 2 = −(𝑎2 ) provided denominator does not become zero. If denominator
becomes zero, then
1
𝑃. 𝐼. = 𝑥 𝑑 𝑓(𝑥)
𝑑𝐷
𝐹(𝐷)
Put 𝐷 2 = −(𝑎2 ) provided denominator does not become zero. If denominator
becomes zero, then
1
𝑃. 𝐼. = 𝑥 2 𝑑2 𝑓(𝑥)
𝑑𝐷 2 𝐹(𝐷)
Put 𝐷 2 = −(𝑎2 ) provided denominator does not become zero.
3. If 𝑓(𝑥) = 𝑥 𝑚
1
𝑃. 𝐼. = 𝑥𝑚
𝐹(𝐷)
Use Binomial Theorem i.e.
𝑛(𝑛−1) 2 𝑛(𝑛−1)(𝑛−2) 3
If |𝑥| < 1 , (1 + 𝑥)𝑛 = 1 + 𝑛𝑥 + 𝑥 + 𝑥 + ⋯ … … … ….
2! 3!
−1 2 3
(1 + 𝑥) = 1 − 𝑥 + 𝑥 − 𝑥 + 𝑥 − ⋯ … … … ….4
i.
ii. (1 − 𝑥)−1 = 1 + 𝑥 + 𝑥 2 + 𝑥 3 + 𝑥 4 + ⋯ … … … ….
iii. (1 + 𝑥)−2 = 1 − 2𝑥 + 3𝑥 2 − 4𝑥 3 + 5𝑥 4 − ⋯ … … … ….
iv. (1 − 𝑥)−2 = 1 + 2𝑥 + 3𝑥 2 + 4𝑥 3 + 5𝑥 4 + ⋯ … … … ….
4. If 𝑓(𝑥) = 𝑒 𝑎𝑥 𝑉, where V is any function.
1 1
𝑃. 𝐼. = 𝑒 𝑎𝑥 𝑉 = 𝑒 𝑎𝑥 ( 𝑉)
𝐹(𝐷) 𝐹(𝐷)
5. General Rule: Factorize 𝐹(𝐷) into linear factors of the form (𝐷 − 𝐴). Apply partial
fractions.
1
𝑓(𝑥) = 𝑒 𝐴𝑥 ∫ 𝑒 −𝐴𝑥 𝑓(𝑥)𝑑𝑥
(𝐷 − 𝐴)
Method of Variation of Parameters
𝑑𝑛 𝑦 𝑑𝑛−1 𝑦 𝑑1 𝑦
𝑎𝑛 𝑥 𝑛 + 𝑎𝑛−1 𝑥 𝑛−1
+ ⋯ … … . . +𝑎1 𝑥 1
+ 𝑎0 𝑦 = 𝑓(𝑥)
𝑑𝑥 𝑛 𝑑𝑥 𝑛−1 𝑑𝑥 1
𝑑2 𝑦
𝑆𝑖𝑚𝑖𝑙𝑎𝑟𝑙𝑦 𝑤𝑒 𝑔𝑒𝑡 𝑥 2 = 𝐷(𝐷 − 1)𝑦
𝑑𝑥 2
Legendre’s equation:
𝑑𝑛 𝑦 𝑑𝑛−1 𝑦 𝑑1 𝑦
𝑎𝑛 (𝑎𝑥 + 𝑏)𝑛 + 𝑎𝑛−1 (𝑎𝑥 + 𝑏)𝑛−1
+ ⋯ … … . . +𝑎1 (𝑎𝑥 + 𝑏)1
+ 𝑎0 𝑦 = 𝑓(𝑥)
𝑑𝑥 𝑛 𝑑𝑥 𝑛−1 𝑑𝑥 1
Is called Cauchy’s homogenous equation.
𝑦𝑖𝑒𝑙𝑑𝑠 𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑧 𝑎
Put (𝑎𝑥 + 𝑏) = 𝑒 𝑧 → 𝑧 = log(𝑎𝑥 + 𝑏) → 𝑑𝑥
= (𝑎𝑥+𝑏)
𝑑2 𝑦
𝑆𝑖𝑚𝑖𝑙𝑎𝑟𝑙𝑦 𝑤𝑒 𝑔𝑒𝑡 (𝑎𝑥 + 𝑏)2 = 𝑎2 𝐷(𝐷 − 1)𝑦
𝑑𝑥 2
𝑑𝑦
Solve (𝑥 − 𝑎) 𝑑𝑥 + 3𝑦 = 12 (𝑥 − 𝑎)3
It is Leibnitz form.
3
Here 𝑃(𝑥) = (𝑥−𝑎) and 𝑄(𝑥) = 12 (𝑥 − 𝑎)2
3
∫(𝑥−𝑎)𝑑𝑥
Integrating factor is 𝑒 ∫ 𝑃(𝑥)𝑑𝑥 = 𝑒 = 𝑒 3𝑙𝑛(𝑥−𝑎) = (𝑥 − 𝑎)3
𝑦𝑖𝑒𝑙𝑑𝑠 𝑦𝑖𝑒𝑙𝑑𝑠
→ (𝑥 − 𝑎)3 𝑦 = ∫(𝑥 − 𝑎)3 12 (𝑥 − 𝑎)2 dx + c → (𝑥 − 𝑎)3 𝑦 = ∫ 12 (𝑥 − 𝑎)5 dx + c
𝑑𝑦
Solve𝑑𝑥 + 4𝑥𝑦 + 𝑥𝑦 3 = 0
𝑑𝑦 𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑦 𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑦
𝑑𝑥
+ 4𝑥𝑦 + 𝑥𝑦 3 = 0 → 𝑑𝑥
+ 4𝑥𝑦 = −𝑥𝑦 3 → 𝑦 −3 𝑑𝑥 + 4𝑥𝑦 −2 = −𝑥
𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑦 𝑑𝑧 𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑦 −1 𝑑𝑧
Put 𝑦 −2 = 𝑧 → −2𝑦 −3 = → 𝑦 −3 =
𝑑𝑥 𝑑𝑥 𝑑𝑥 2 𝑑𝑥
−1 𝑑𝑧 𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑧
Equation becomes 2 𝑑𝑥
+ 4𝑥𝑧 = −𝑥 → 𝑑𝑥
− 8𝑥𝑧 = 2𝑥 . It is Leibnitz form.
𝑦𝑖𝑒𝑙𝑑𝑠 2 2 𝑦𝑖𝑒𝑙𝑑𝑠 2
→ 𝑧𝑒 −4𝑥 = ∫ 2𝑥 (𝑒 −4𝑥 )dx + c → 𝑦 −2 𝑒 −4𝑥 = ∫(𝑒 −4𝑡 )dt + c
2
𝑦𝑖𝑒𝑙𝑑𝑠 2 𝑒 −4𝑡 𝑦𝑖𝑒𝑙𝑑𝑠 2 𝑒 −4𝑥 𝑦𝑖𝑒𝑙𝑑𝑠 −1 2
→ 𝑦 −2 𝑒 −4𝑥 = −4
+𝑐→ 𝑦 −2 𝑒 −4𝑥 = −4
+𝑐→ 𝑦 −2 = 4
+ 𝑐𝑒 −4𝑥
𝑑𝑦
Solve𝑑𝑥 − 𝑦 = 𝑦 2 (sin 𝑥 + cos 𝑥)
𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑦 𝑑𝑧 𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑦 𝑑𝑧
Put 𝑦 −1 = 𝑧 → −𝑦 −2 𝑑𝑥 = 𝑑𝑥 → 𝑦 −2 𝑑𝑥 = − 𝑑𝑥
𝑑𝑧 𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑧
Equation becomes − 𝑑𝑥 − 𝑧 = (sin 𝑥 + cos 𝑥) → 𝑑𝑥
+ 𝑧 = −(sin 𝑥 + cos 𝑥) . It is Leibnitz
form.
𝑦𝑖𝑒𝑙𝑑𝑠 𝑦𝑖𝑒𝑙𝑑𝑠 𝑒𝑥 𝑒𝑥 𝑒𝑥
→ 𝑧𝑒 𝑥 = ∫ −𝑒 𝑥 (sin 𝑥 + cos 𝑥)dx + c → = − [ (sin 𝑥 − cos 𝑥) + (cos 𝑥 + sin 𝑥)] + c
𝑦 2 2
𝑦𝑖𝑒𝑙𝑑𝑠 1 𝑐
→ 𝑦
= − sin 𝑥 + 𝑒 𝑥
𝑑𝑦 𝑦
Solve 𝑑𝑥 = 𝑥+√𝑥𝑦
𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑥 𝑑𝑧
Put 𝑥 1/2 = 𝑧 → 𝑥 −1/2 𝑑𝑦 = 2 𝑑𝑦
𝑑𝑧 𝑧 1 𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑧 𝑧 1
Equation becomes 2 − = → − = . It is Leibnitz form.
𝑑𝑦 𝑦 √𝑦 𝑑𝑦 2𝑦 2√𝑦
1 1
Here 𝑃(𝑦) = − and 𝑄(𝑦) =
2𝑦 2√𝑦
1
∫(−2𝑦)𝑑𝑦 1
Integrating factor is 𝑒 ∫ 𝑃(𝑦)𝑑𝑦 = 𝑒 = 𝑒 𝑙𝑜𝑔(1/√𝑦)=
√𝑦
𝑑𝑦
Solvetan 𝑦 𝑑𝑥 + tan 𝑥 = cos 𝑦 𝑐𝑜𝑠 2 𝑥 .
𝑑𝑦
Sol. tan 𝑦 𝑑𝑥 + tan 𝑥 = cos 𝑦 𝑐𝑜𝑠 2 𝑥
𝑑𝑦
Dividing both sides by cos 𝑦 , we get sec 𝑦 tan 𝑦 𝑑𝑥 + sec 𝑦 tan 𝑥 = 𝑐𝑜𝑠 2 𝑥
𝑦𝑖𝑒𝑙𝑑𝑠 𝑑𝑦 𝑑𝑧
Put sec 𝑦 = 𝑧 → sec 𝑦 tan 𝑦 𝑑𝑥 = 𝑑𝑦
𝑑𝑧
Equation becomes 𝑑𝑦
− 𝑧 tan 𝑥 = 𝑐𝑜𝑠 2 𝑥 . It is Leibnitz form.
𝑦𝑖𝑒𝑙𝑑𝑠 𝑦𝑖𝑒𝑙𝑑𝑠
→ 𝑧 cos 𝑥 = ∫ cos 𝑥 𝑐𝑜𝑠 2 𝑥 dx + c → 𝑧 cos 𝑥 = ∫ 𝑐𝑜𝑠 3 𝑥 dx + c
𝑦𝑖𝑒𝑙𝑑𝑠 1 sin 3x
→ sec 𝑦 cos 𝑥 = [ + 3 sin x] + c
4 3
Solve𝑦 ′′ − 4𝑦 ′ − 5𝑦 = 0
𝑦𝑖𝑒𝑙𝑑𝑠
→ 𝐷 = −1, 5
Particular Integral = 0
Solve 𝑦 ′′ + 2𝑦 ′ + 2𝑦 = 0
Auxiliary equation is 𝐷 2 + 2𝐷 + 2 = 0
𝑦𝑖𝑒𝑙𝑑𝑠
→ 𝐷 = −1 ± 𝑖
Particular Integral = 0
Particular Integral = 0
𝑑2 𝑥 𝑑𝑥
Solve 𝑑𝑡 2 + 5 𝑑𝑡 + 6𝑥 = 0 , given x(0)=0 and x’(0)=15
𝑦𝑖𝑒𝑙𝑑𝑠
Sol. Auxiliary equation is 𝐷 2 + 5𝐷 + 6 = 0 → 𝐷 = −2, −3
−2𝑡
Therefore complementary function is 𝑐1 𝑒 + 𝑐2 𝑒 −3𝑡
1
Particular Integral= 2 (0) = 0
𝐷 +5𝐷+6
Complete Solution = Complementary function + Particular Integral
𝑦𝑖𝑒𝑙𝑑𝑠
→ 𝑥 = 𝑐1 𝑒 −2𝑡 + 𝑐2 𝑒 −3𝑡
1 1
Particular Integral = (2𝐷−1)2 𝑒 𝑥/2 (𝑃𝑢𝑡 𝐷 = 2)
1 𝑥 2 𝑒 𝑥/2
∴ Particular Integral = 𝑥 2 8 𝑒 𝑥/2 = 8
𝑥
Solve (D2-3D+2)y=2excos2
𝑦𝑖𝑒𝑙𝑑𝑠
Sol.Auxiliary equation is 𝐷 2 − 3𝐷 + 2 = 0 → 𝐷 = 2, 1
Therefore complementary function is 𝑐1 𝑒 2𝑥 + 𝑐2 𝑒 𝑥
1 𝑥 1 𝑥 1 𝑥
Particular Integral = 2 2𝑒 𝑥 cos = 2𝑒 𝑥 [((𝐷+1)2 cos ] = 2𝑒 𝑥 [(𝐷2 ] cos
𝐷 −3𝐷+2 2 −3(𝐷+1)+2) 2 −𝐷) 2
1 1 𝑥 1 𝑥
(𝑃𝑢𝑡 𝐷 2 = − 4) = 2𝑒 𝑥 [ 1 ] cos 2 = −8𝑒 𝑥 [1+4𝐷] cos 2
(− −𝐷)
4
1 − 4𝐷 𝑥 1
= −8𝑒 𝑥 [ 2 ] cos (𝑃𝑢𝑡 𝐷 2 = − )
1 − 16𝐷 2 4
8 𝑥 1 𝑥
= − 5 𝑒 𝑥 (cos 2 − 2 sin 2)
1 1 2(𝐷−1)
Particular Integral = (𝐷−1)2 𝑥 sin 𝑥 = 𝑥 [𝐷2 −2𝐷+1 sin 𝑥] − [𝐷2 −2𝐷+1 sin 𝑥] (𝑃𝑢𝑡 𝐷 2 = −1)
1 2(𝐷 − 1) 1 1
= 𝑥[ sin 𝑥] − [ sin 𝑥] = 𝑥 cos 𝑥 − (𝐷 − 1) cos 𝑥 = 𝑥 cos 𝑥 + sin 𝑥 + cos 𝑥
−2𝐷 −2𝐷 2 2
1 2𝐷+5
= 𝑒 𝑥 [𝑥 [𝐷2 +5𝐷+6 sin 𝑥] − [(𝐷2 +5𝐷+6)2 sin 𝑥]] (𝑃𝑢𝑡 𝐷 2 = −1)
1 2𝐷 + 5
= 𝑒 𝑥 [𝑥 [ sin 𝑥] − [ sin 𝑥]]
5𝐷 + 5 (5𝐷 + 5)2
𝑥 (𝐷 − 1) 1 (2𝐷 + 5)
= 𝑒𝑥 [ [ 2 sin 𝑥] − [ 2 sin 𝑥]] (𝑃𝑢𝑡 𝐷 2 = −1)
5 𝐷 −1 25 𝐷 + 2𝐷 + 1
−𝑥 1 (2𝐷+5) −𝑥 1
= 𝑒 𝑥 [ 10 [(𝐷 − 1) sin 𝑥] − 25 [ 2𝐷
sin 𝑥]] = 𝑒 𝑥 [ 10 [(𝐷 − 1) sin 𝑥] + 50 [(2𝐷 + 5) cos 𝑥]]
−𝑥 1
= 𝑒 𝑥 [ 10 (cos 𝑥 − sin 𝑥) + 50 [(−2 sin 𝑥 + 5 cos 𝑥)]]
1 (2𝐷 + 4)
= 𝑒 2𝑥 [𝑥 [ cos 𝑥] − [ cos 𝑥]] (𝑃𝑢𝑡 𝐷 2 = −1)
𝐷 2 + 4𝐷 + 13 (𝐷 2 + 4𝐷 + 13)2
1 (2𝐷 + 4)
= [𝑒 2𝑥 𝑥 [ cos 𝑥] − [ cos 𝑥]]
4𝐷 + 12 (4𝐷 + 12)2
𝑥 (𝐷 − 3) 1 (2𝐷 + 4)
= 𝑒 2𝑥 [ [ 2 cos 𝑥] − [ 2 cos 𝑥]] (𝑃𝑢𝑡 𝐷 2 = −1)
4 𝐷 −9 16 𝐷 + 6𝐷 + 9
𝑥 (𝐷−3) 1 (2𝐷+4)
= 𝑒 2𝑥 [4 [ −10
cos 𝑥] − 16 [ 6𝐷+8 cos 𝑥]]
𝑥 (− sin 𝑥 − 3 cos 𝑥) 1
= 𝑒 2𝑥 [ [ ]+ [(12𝐷 2 + 8𝐷 − 32) cos 𝑥]]
4 −10 1600
𝑥 1
= 𝑒 2𝑥 [ (sin 𝑥 + 3 cos 𝑥) + (−12 cos 𝑥 − 8 sin 𝑥 − 32 cos 𝑥)]
40 1600
𝑥 1
= 𝑒 2𝑥 [ (sin 𝑥 + 3 cos 𝑥) − (44 cos 𝑥 + 8 sin 𝑥)]
40 1600
Sol. The given equation is 𝑦 ′′ + 16𝑦 = 32 sec 2𝑥by method of variation of parameters.
𝑦𝑖𝑒𝑙𝑑𝑠
Auxiliary equation is 𝐷 2 + 16 = 0 → 𝐷 = ±4𝑖
Particular Integral = 𝑢𝑦1 + 𝑣𝑦2 = 8 cos 2𝑥 cos 4𝑥 + (8 sin 2𝑥 − 4𝑙𝑜𝑔|sec2 𝑥 + tan2 𝑥|) sin4 𝑥
Particular Integral = 𝑢𝑦1 + 𝑣𝑦2 = −(𝑙𝑜𝑔|sec 𝑥 + tan 𝑥| − sin 𝑥)𝑒 𝑥 cos 𝑥 − cos 𝑥 𝑒 𝑥 sin 𝑥
2
Sol. 𝑦 ′′ − 𝑦 = 1+𝑒 𝑥
𝑦𝑖𝑒𝑙𝑑𝑠
Auxiliary equation is 𝐷 2 − 1 = 0 → 𝐷 = ±1
Here 𝑦1 = 𝑒 −𝑥 𝑎𝑛𝑑 𝑦2 = 𝑒 𝑥
𝑦1 𝑦2 𝑒 −𝑥 𝑒𝑥
𝑊 = |𝑦 ′ 𝑦2 ′| = | |=2
1 −𝑒 −𝑥 𝑒𝑥
𝑦2 𝑋𝑑𝑥 𝑒𝑥 2 𝑒𝑥
𝑢 = −∫ 𝑊
= −∫ 2 (1+𝑒 𝑥 )
𝑑𝑥 = − ∫ (1+𝑒 𝑥 ) 𝑑𝑥 = −log(1 + 𝑒 𝑥 )
𝑦1 𝑋𝑑𝑥 𝑒 −𝑥 2 𝑒 −𝑥
𝑣=∫ 𝑊
=∫ 2 (1+𝑒 𝑥 )
𝑑𝑥 = ∫ (1+𝑒 𝑥 ) 𝑑𝑥 Put 𝑒 −𝑥 = 𝑡. ∴ 𝑒 −𝑥 𝑑𝑥 = −𝑑𝑡
𝑒 −𝑥 𝑡2 𝑡2 − 1 + 1 1 𝑡2
∫ 𝑑𝑥 = − ∫ 𝑑𝑡 = − ∫ 𝑑𝑡 = − ∫(𝑡 − 1 + )𝑑𝑡 = − +𝑡
(1 + 𝑒 𝑥 ) (𝑡 + 1) (𝑡 + 1) (𝑡 + 1) 2
− log((𝑡 + 1)
𝑒 −2𝑥
=− + 𝑒 −𝑥 − log((𝑒 −𝑥 + 1)
2
𝑒 −2𝑥
Particular Integral = 𝑢𝑦1 + 𝑣𝑦2 = −log(1 + 𝑒 𝑥 ) 𝑒 −𝑥 + 𝑒 𝑥 [− + 𝑒 −𝑥 − log((𝑒 −𝑥 + 1)]
2
Solve4𝑥 2 𝑦 ′′ + 𝑦 = 0
𝑑2 𝑦
∴ 𝑊𝑒 𝑔𝑒𝑡 𝑥 2 = 𝐷(𝐷 − 1)𝑦
𝑑𝑥 2
𝑦𝑖𝑒𝑙𝑑𝑠
Hence the given equation becomes 4𝐷(𝐷 − 1)𝑦 + 𝑦 = 0 → (4𝐷 2 − 4𝐷 + 1)𝑦 = 0
𝑦𝑖𝑒𝑙𝑑𝑠 1 1
Auxiliary equation is 4𝐷 2 − 4𝐷 + 1 = 0 → 𝐷 = 2,2
1 1
Therefore complementary function is 𝑐1 𝑒 2𝑧 + 𝑧𝑐2 𝑒 2𝑧 = √𝑥(𝑐1 + 𝑐2 log 𝑥)
Solve 𝑥 2 𝑦 ′′ + 2𝑥𝑦 ′ − 2𝑦 = 0
𝑑2 𝑦 2
𝑆𝑖𝑚𝑖𝑙𝑎𝑟𝑙𝑦 𝑤𝑒 𝑔𝑒𝑡 𝑥 = 𝐷(𝐷 − 1)𝑦
𝑑𝑥 2
𝑦𝑖𝑒𝑙𝑑𝑠
Hence the given equation becomes 𝐷(𝐷 − 1)𝑦 + 2𝐷𝑦 − 2𝑦 = 0 → (𝐷 2 + 𝐷 − 2)𝑦 = 0
𝑦𝑖𝑒𝑙𝑑𝑠
Auxiliary equation is 𝐷 2 + 𝐷 − 2 = 0 → 𝐷 = 1, −2
𝑐
Therefore complementary function is 𝑐1 𝑒 𝑧 + 𝑐2 𝑒 −2𝑧 = 𝑐1 𝑥 + 𝑥22
Particular Integral = 0
Solve 𝑥 2 𝑦 ′′ + 𝑥𝑦 ′ + 𝑦 = 𝑥
𝑑2 𝑦 2
𝑆𝑖𝑚𝑖𝑙𝑎𝑟𝑙𝑦 𝑤𝑒 𝑔𝑒𝑡 𝑥 = 𝐷(𝐷 − 1)𝑦
𝑑𝑥 2
𝑦𝑖𝑒𝑙𝑑𝑠
Hence the given equation becomes 𝐷(𝐷 − 1)𝑦 + 𝐷𝑦 + 𝑦 = 𝑒 𝑧 → (𝐷 2 + 1)𝑦 = 𝑒 𝑧
𝑦𝑖𝑒𝑙𝑑𝑠
Auxiliary equation is 𝐷 2 + 1 = 0 → 𝐷 = ±𝑖
1 𝑥
= 𝑒𝑧 =
2 2
Solve 𝑥 2 𝑦 ′′ + 5𝑥𝑦 ′ + 3 = ln 𝑥
𝑑2 𝑦
𝑆𝑖𝑚𝑖𝑙𝑎𝑟𝑙𝑦 𝑤𝑒 𝑔𝑒𝑡 𝑥 2 = 𝐷(𝐷 − 1)𝑦
𝑑𝑥 2
Hence the given equation becomes
𝑦𝑖𝑒𝑙𝑑𝑠
𝐷(𝐷 − 1)𝑦 + 5𝐷𝑦 + 3𝑦 = 𝑧 → (𝐷 2 + 4𝐷 + 3)𝑦 = 𝑧
𝑦𝑖𝑒𝑙𝑑𝑠
Auxiliary equation is 𝐷 2 + 4𝐷 + 3 = 0 → 𝐷 = −3, −1
𝑐1 𝑐
Therefore complementary function is 𝑐1 𝑒 −𝑧 + 𝑐2 𝑒 −3𝑧 = 𝑥
+ 𝑥23
−1
1 1 1 𝐷2 4𝐷
Particular Integral = 𝑧= 𝐷2 4𝐷
𝑧 = [1 + ( + )] 𝑧
𝐷 2 +4𝐷+3 3(1+ + ) 3 3 3
3 3
2
1 𝐷 2 4𝐷 𝐷 2 4𝐷 1 4
= [1 − − + ( + ) − ⋯ … … ] 𝑧 = [𝑧 − 0 − + 0 − 0 … … … ]
3 3 3 3 3 3 3
𝑧 4 ln 𝑥 4
= − = −
3 9 3 9
𝑑2 𝑦 𝑑𝑦
Solve Cauchy’s homogenous linear eq.𝑥 2 𝑑𝑥2 − 𝑥 𝑑𝑥 + 2𝑦 = 𝑥𝑙𝑜𝑔𝑥
𝑑2 𝑦
𝑆𝑖𝑚𝑖𝑙𝑎𝑟𝑙𝑦 𝑤𝑒 𝑔𝑒𝑡 𝑥 2 = 𝐷(𝐷 − 1)𝑦
𝑑𝑥 2
Hence the given equation becomes
𝑦𝑖𝑒𝑙𝑑𝑠
𝐷(𝐷 − 1)𝑦 − 𝐷𝑦 + 2𝑦 = 𝑒 𝑧 𝑧 → (𝐷 2 − 2𝐷 + 1)𝑦 = 𝑒 𝑧 𝑧
𝑦𝑖𝑒𝑙𝑑𝑠
Auxiliary equation is 𝐷 2 − 2𝐷 + 1 = 0 → 𝐷 = 1,1
1 𝑧3 (log 𝑥)3
= 𝑒𝑧 [ 𝑧] = 𝑒 𝑧
= 𝑥
𝐷2 6 6
𝑑2 𝑦 𝑑𝑦
Solve Cauchy’s homogenous linear eq.𝑥 2 𝑑𝑥2 + 5𝑥 𝑑𝑥 − 5𝑦 = 24𝑥𝑙𝑜𝑔𝑥
𝑑2 𝑦
𝑆𝑖𝑚𝑖𝑙𝑎𝑟𝑙𝑦 𝑤𝑒 𝑔𝑒𝑡 𝑥 2 = 𝐷(𝐷 − 1)𝑦
𝑑𝑥 2
Hence the given equation becomes
𝑦𝑖𝑒𝑙𝑑𝑠
𝐷(𝐷 − 1)𝑦 + 5𝐷𝑦 − 5𝑦 = 24𝑒 𝑧 𝑧 → (𝐷 2 + 4𝐷 − 5)𝑦 = 24𝑒 𝑧 𝑧
𝑦𝑖𝑒𝑙𝑑𝑠
Auxiliary equation is 𝐷 2 + 4𝐷 − 5 = 0 → 𝐷 = −5, 1
𝑐
Therefore complementary function is 𝑐1 𝑒 −5𝑧 + 𝑐2 𝑒 𝑧 = 𝑥15 + 𝑐2 𝑥
1 1
Particular Integral = 𝐷2 +4𝐷−5 24𝑧𝑒 𝑧 = 24𝑒 𝑧 [((𝐷+1)2 +4(𝐷+1)−5) 𝑧]
1 24𝑒 𝑧 𝐷 −1 4𝑒 𝑧 𝐷 𝐷2 4𝑒 𝑧 1
= 24𝑒 𝑧 [ 2
𝑧] = (1 + ) 𝑧 = (1 − + … … … ) 𝑧 = (𝑧 − )
𝐷 + 6𝐷 6𝐷 6 𝐷 6 36 𝐷 6
2 (log 𝑥)2
𝑧 𝑧 log 𝑥
= 4𝑒 𝑧 ( − ) = 4𝑥 ( − )
2 6 2 6
𝑦𝑖𝑒𝑙𝑑𝑠
Sol. Put (1 + 𝑥) = 𝑒 𝑧 → 𝑧 = log(1 + 𝑥)
𝑦𝑖𝑒𝑙𝑑𝑠
The given equation becomes [𝐷(𝐷 − 1) + 𝐷 + 1]𝑦 = 4 cos 𝑧 → (𝐷 2 + 1)𝑦 = 4 cos 𝑧
𝑦𝑖𝑒𝑙𝑑𝑠
Auxiliary equation is 𝐷 2 + 1 = 0 → 𝐷 = ±𝑖
1
But it is case of failure. ∴ Particular Integral = 𝑧 4 cos 𝑧
2𝐷
SIMULTANEOUS EQUATIONS
(𝐷 2 + 10𝐷 + 24)𝑥 = 6𝑒 𝑡 − 𝑒 5𝑡
𝑦𝑖𝑒𝑙𝑑𝑠
Auxiliary equation is 𝐷 2 + 10𝐷 + 24 = 0 → 𝐷 = −4, −6
1 6𝑒 𝑡 𝑒 5𝑡
Particular Integral = 𝐷2 +10𝐷+24 (6𝑒 𝑡 − 𝑒 5𝑡 ) = 35
− 99
6𝑒 𝑡 𝑒 5𝑡
Putting the value of x in (i), we get (𝐷 + 5)(𝑐1 𝑒 −4𝑡 + 𝑐2 𝑒 −6𝑡 + 35
− 99
)+ 𝑦 = 𝑒𝑡
𝑦𝑖𝑒𝑙𝑑𝑠 6𝑒 𝑡 5𝑒 5𝑡 30𝑒 𝑡 5𝑒 5𝑡
→ −4𝑐1 𝑒 −4𝑡 − 6𝑐2 𝑒 −6𝑡 + − + 5𝑐1 𝑒 −4𝑡 + 5𝑐2 𝑒 −6𝑡 + − + 𝑦 = 𝑒𝑡
35 99 35 99
𝑦𝑖𝑒𝑙𝑑𝑠 36𝑒 𝑡 10𝑒 5𝑡 𝑦𝑖𝑒𝑙𝑑𝑠 𝑒 𝑡 10𝑒 5𝑡
→ 𝑐1 𝑒 −4𝑡 − 𝑐2 𝑒 −6𝑡 +
− + 𝑦 = 𝑒𝑡 → 𝑦 = −𝑐1 𝑒 −4𝑡 + 𝑐2 𝑒 −6𝑡 − +
35 99 35 99
Solve D𝑥 + 2𝑦 = − sin 𝑡 and −2𝑥 + 𝐷𝑦 = cos 𝑡 simultaneously.
Sol. Given equations are D𝑥 + 2𝑦 = − sin 𝑡 …………..(i) and −2𝑥 + 𝐷𝑦 = cos 𝑡 …………..(ii)
Multiplying (i) by 2 and (ii) by D and adding (ii) and (i), we get
(𝐷 2 − 4)𝑦 = −3 sin 𝑡
𝑦𝑖𝑒𝑙𝑑𝑠
Auxiliary equation is 𝐷 2 − 4 = 0 → 𝐷 = −2, 2
3 sin 𝑡
Putting the value of y in (ii), we get −2𝑥 + 𝐷 (𝑐1 𝑒 −2𝑡 + 𝑐2 𝑒 2𝑡 + 5
) = cos 𝑡
𝑦𝑖𝑒𝑙𝑑𝑠 3 cos 𝑡
→ −2𝑥 + (−2𝑐1 𝑒 −2𝑡 + 2𝑐2 𝑒 2𝑡 + ) = cos 𝑡
5
𝑦𝑖𝑒𝑙𝑑𝑠 cos 𝑡
→ 𝑥= + 𝑐1 𝑒 −2𝑡 − 𝑐2 𝑒 2𝑡
5
SPECIAL FUNCTIONS
2 x
t 2dt
Ans. Error function = erfx e .
0
Properties:
1. erf x erfx
2. erf x erfc x 1
Q.2 What are ordinary, regular & irregular singular points of an ordinary differential function?
d2y
Px Q x y 0 ; if both functions Px and Qx are
dy
2
dx dx
analytic at x x0 .
Irregular Singular Point: A singular point which is not regular, is called Irregular Singular Point.
y a0 a1 x a2 x 2 a3 x3
dy d2y
Step – IIFind and
dx dx2
dy d2y
Step – III Substitute the value of and in the given Differential equation.
dx dx2
1 x ddxy 2 x dy
2
Ans.
2
2
nn 1 y 0
dx
d2y
x nn 1 y 0
dy
x2 2
dx dx
m 1
Ans. Step – I Assume y a0 x a1 x
m
a2 x m 2
dy d2y
Step – IIFind the value of and
dx dx2
dy d2y
Step – III Substitute the value of y, and in the given differential equation.
dx dx2
Step – IV Equate the coefficient of lowest power of x equal to zero. This gives a
quadratic equation in m, which is known as Indicial equation.
Ans. Pn x
1 dn 2
n
2 n! dx n
x 1
n
Q.8 Show that x4
1
8P4 x 20 P2 x 7 P0 x .
35
Ans. P4 x
1
8
35 x 4 30 x 2 3
P2 x
1 2
2
3x 1
P0 x 1
Hence
1 1
35 8
1
8 35 x 4 30 x 2 3 20 3x 2 1 7 x 4
2
x 2 P2 x ;
2 1
x 3 P3 x x ; x P1 x ; P0 x 1
2 3
Ans.
5 5 3 3
2 3 2 1
x3 5 x 2 x 2 P3 x x 5 P2 x P1 x 2 P0 x
5 5 3 3
P3 x P2 x P1 x P0 x
2 10 8 1
5 3 5 3
Ans.
d n
dx
x J n x x n J n 1 x
Put n = 0,
d
J 0 x J1 x
dx
J 0 x J1 x
0 when α ß
0 xJ n x J n ßxdx 1 J
1
n 1 when α ß
2
2 x
0 e
t 2 dt
Ans. As we know erfx
Put x0
2 0 t 2 dt
erf 0
0 e 0
Ans. J n x represents the Bessel function which is the solution of the Bessel differential equation.
J 1 x
2
sin x .
2
x
Q.14 Prove that Pn 1 1 .
1
hn Pn x 1 2 xh h2
Ans. As we know, 2
n0
Put x = 1,
1
hn Pn 1 1 2h h2
2
n0
1 h
1
1 h h2 hn
Prove that Pn x 1 Pn x .
n
Q.15
1
As we know h Pn x 1 2 xh h
n 2
Ans. 2 ………………………………… (1)
n0
Replace x by –x in (1),
h Pn x 1 2 xh h2 2 ……………………………………………. (2)
1
n
n 0
Replace h by –h in (1),
h Pn x 1 2 xh h2 2 …………………………………………… (3)
1
n
n 0
n0 n0
Pn x 1 Pn x
n
0 if m n
1 Pm x Pn x dx 2
1
Ans.
if m n
2n 1
1 2 2 , ß 2 r 2 4 , 2r 4 , 2 n2 r 2 4
3 5
, ß 1, r 2 , n
2 4
3
n is not an integer. y x 2 c1 J 5 x 2 c2 J 5 x 2
4 4
x 1
z
2 z
Ans. The function e is called generating function.
x 1
z
e 2 z
= z n J n x
n
x 1
z
J n x is the coefficient of z in the expansion of e n 2 z
.
erf x
2 x
0 e
t 2
dt
Ans.
erf x
2 x
0 et
2
dt
2 x
0 e
t 2
dt
erfx
Prove that Pn 1 1 .
n
Q.20
We know that Pn x 1 Pn x
n
Ans.
Put x = 1,
Pn 1 1 Pn 1
n
1
n
As Pn 1 1
nn 1
Q.21 Prove that Pn1
2
Put x = 1,
Ans. We know J 0 J1
J 0 J1 J 3
1 1
2 4
3 1
J 0 J 3
4 4
4 J 0x 3J 0 x J 3 x 0
Q.23 Prove that J n3 J n5
2
n 4J n 4 .
x
Replace n by n + 4, we get:
2
n 4J n 4 J n 3 J n 5
x
As erf 0 0
erfc 0 1
xy 3 y xy 0 .
Ans. 1 2 3 , ß 2 r 2 1 , 2r 2 , 2 n2r 2 0
2, ß r 1, n2
𝑑
1. Prove that 𝑑𝑥 ⌈𝑥 𝑛 𝐽𝑛 (𝑥)⌉ = 𝑥 𝑛 𝐽𝑛−1 (𝑥)
Proof: Since
(−1)𝑘 𝑥 𝑛+2𝑘
𝐽𝑛 (𝑥) = ∑∞
𝑘=0 𝑘!Γ(𝑛+𝑘+1) (2 )
∞
𝑛
(−1)𝑘 𝑥 2𝑛+2𝑘
𝑥 𝐽𝑛 (𝑥) = ∑
2𝑛+2𝑘 𝑘! Γ(𝑛 + 𝑘 + 1)
𝑘=0
∞ ∞
𝑑 𝑛 (−1)𝑘 (2𝑛 + 2𝑘)𝑥 2𝑛+2𝑘−1 𝑛
(−1)𝑘 (𝑛 + 𝑘)𝑥 𝑛+2𝑘−1
⌈𝑥 𝐽𝑛 (𝑥)⌉ = ∑ = 𝑥 ∑ 𝑛+2𝑘−1
𝑑𝑥 2𝑛+2𝑘 𝑘! Γ(𝑛 + 𝑘 + 1) 2 𝑘! (𝑛 + 𝑘)Γ(𝑛 + 𝑘)
𝑘=0 𝑘=0
∞ ∞
𝑛
(−1)𝑘 𝑥 𝑛+2𝑘−1 (−1)𝑘 𝑥 𝑛+2𝑘−1
=𝑥 ∑ ( ) = 𝑥𝑛 ∑ ( )
𝑘! Γ(𝑛 + 𝑘) 2 𝑘! Γ(𝑛 − 1 + 𝑘 + 1) 2
𝑘=0 𝑘=0
= 𝑥 𝑛 𝐽𝑛−1 (𝑥)
𝑑
II) Prove that ⌈𝑥 −𝑛 𝐽𝑛 (𝑥)⌉ = −𝑥 −𝑛 𝐽𝑛+1 (𝑥)
𝑑𝑥
(−1)𝑘 𝑥 𝑛+2𝑘
Since 𝐽𝑛 (𝑥) = ∑∞
𝑘=0 𝑘!Γ(𝑛+𝑘+1) (2 )
∞
−𝑛
(−1)𝑘 𝑥 2𝑘
𝑥 𝐽𝑛 (𝑥) = ∑
2𝑛+2𝑘 𝑘! Γ(𝑛 + 𝑘 + 1)
𝑘=0
∞
𝑑 −𝑛 (−1)𝑘 (2𝑘)𝑥 2𝑘−1
⌈𝑥 𝐽𝑛 (𝑥)⌉ = ∑ 𝑛+2𝑘
𝑑𝑥 2 𝑘(𝑘 − 1)! Γ(𝑛 + 𝑘 + 1)
𝑘=0
∞
−𝑛
(−1)𝑘−1 𝑥 𝑛+2𝑘−1
=𝑥 ∑ 𝑛+3𝑘−1
2 (𝑘 − 1)! (𝑛 + 𝑘)Γ(𝑛 + 𝑘 + 1)
𝑘=0
∞
(−1)𝑟 𝑥 𝑛+2𝑘+1
= −𝑥 −𝑛 ∑ 𝑤ℎ𝑒𝑟𝑒 𝑟 = 𝑘 − 1
2𝑛+2𝑟+1 𝑟! Γ(𝑛 + 𝑟 + 2)
𝑘=0
∞
−𝑛
(−1)𝑟 𝑥 𝑛+1+2𝑟
= −𝑥 ∑ ( ) = −𝑥 −𝑛 𝐽𝑛+1 (𝑥)
𝑟! Γ(𝑛 + 1 + 𝑟 + 1) 2
𝑘=0
𝑛
III) Prove that 𝐽𝑛 , (𝑥) + 𝐽 (𝑥) = 𝐽𝑛−1 (𝑥)
𝑥 𝑛
Proof:
𝑑 𝑛
⌈𝑥 𝐽𝑛 (𝑥)⌉ = 𝑥 𝑛 𝐽𝑛−1 (𝑥)
𝑑𝑥
𝑛
IV) Prove that 𝐽𝑛 , (𝑥) − 𝐽 (𝑥)
𝑥 𝑛
= −𝐽𝑛+1 (𝑥)
𝑑
Proof: ⌈𝑥 𝑛 𝐽𝑛 (𝑥)⌉ = 𝑥 𝑛 𝐽𝑛−1 (𝑥)
𝑑𝑥
Multiplying by 𝑥 𝑛 , we have
𝑛
𝐽𝑛 , (𝑥) − 𝐽 (𝑥) = −𝐽𝑛+1 (𝑥)_ ______________________________________(2)
𝑥 𝑛
2𝑛
VI) Prove that 𝐽 (𝑥) = 𝐽𝑛−1 (𝑥) + 𝐽𝑛+1 (𝑥)
𝑥 𝑛
2𝑛
We get 𝐽 (𝑥)
𝑥 𝑛
= 𝐽𝑛−1 (𝑥) + 𝐽𝑛+1 (𝑥)
𝑛=0
∞ ∞
𝑛=0 𝑛=0
𝑛−1
Equating coefficients of ℎ on both sides,
x𝑃𝑛−1 (𝑥) − 𝑃𝑛−2 (𝑥) = 𝑛𝑃𝑛 (𝑥) − 2𝑥(𝑛 − 1)𝑃𝑛−1 (𝑥) + (𝑛 − 2)𝑃𝑛−2 (𝑥)
𝑛𝑃𝑛 (𝑥) = (2𝑛 − 1)𝑥𝑃𝑛−1 (𝑥) − (𝑛 − 1) 𝑃𝑛−2 (𝑥)
Replacing n by n+1, we get
(𝑛 + 1)𝑃𝑛+1 (𝑥) = (2𝑛 + 1)𝑥𝑃𝑛 (𝑥) − 𝑛𝑃𝑛−1 (𝑥)
−1
ii) we know that : (1 − 2𝑥ℎ + ℎ2 ) 2 = ∑∞ 𝑛
𝑛=0 ℎ 𝑃𝑛 (𝑥)
𝑛=0
∞ ∞
𝑛=0 𝑛=0
∞
−1 −3
(1 − 2𝑥ℎ + ℎ2 ) 2 (−2ℎ) = ∑ ℎ𝑛 𝑃′ 𝑛 (𝑥)
2
𝑛=0
∞
−3
(𝑥 − ℎ)(1 − 2𝑥ℎ + ℎ2 ) 2 = (𝑥 − ℎ) ∑ ℎ𝑛−1 ) 𝑃′ 𝑛 (𝑥)
𝑛=0
iii) From recurrence relation (i) (2𝑛 + 1)𝑥𝑃𝑛 (𝑥) = (𝑛 + 1)𝑃𝑛+1 (𝑥) + 𝑛𝑃𝑛−1 (𝑥)
(2𝑛 + 1)[𝑥𝑃𝑛 , (𝑥) + 𝑃𝑛 (𝑥)] = (𝑛 + 1)𝑃, 𝑛+1 (𝑥) + 𝑛𝑃, 𝑛−1 (𝑥) ………………….(1)
From (1) (2𝑛 + 1)[𝑛𝑃𝑛 (𝑥) + 𝑃, 𝑛−1 (𝑥) + 𝑃𝑛 (𝑥)] = (𝑛 + 1)𝑃, 𝑛+1 (𝑥) + 𝑛𝑃, 𝑛−1 (𝑥)
Ans. P4 x
1
8
35 x 4 30 x 2 3
P2 x
1 2
2
3x 1
P0 x 1
Hence
1 1
35 8
1
8 35 x 4 30 x 2 3 20 3x 2 1 7 x 4
2
2 3 2 1
x3 5 x 2 x 2 P3 x x 5 P2 x P1 x 2 P0 x
5 5 3 3
P3 x P2 x P1 x P0 x
2 10 8 1
5 3 5 3