1.1 Differential Equations and Their Classifications
1.1 Differential Equations and Their Classifications
1.1 Differential Equations and Their Classifications
Chapter 1
In this chapter, an overview of differential equations will be presented. A special attention is paid
to classifications and solution of linear differential equations.
𝑑2 𝑦 𝑑𝑦
1)
2
+ 𝑥𝑦( )2 = 0
𝑑𝑥 𝑑𝑥
𝑑 3 𝑦 𝑑 2 𝑦 𝑑𝑦
2) + + = 𝑦 2 tan 𝑡
𝑑𝑡 3 𝑑𝑡 2 𝑑𝑡
P a g e |2
𝜕𝑣 𝜕𝑣
1) + =𝑣
𝜕𝑠 𝜕𝑡
𝜕𝑧 𝜕2𝑧
2) +𝑧= 2
𝜕𝑥 𝜕𝑦
a) Every dependent variable and every derivative involved occur to the first degree only, and
b) No products of dependent variables and/or derivatives occur.
A differential equation which is not linear is called a non-linear differential equation. Some
examples for linearity and non-linearity of differential equations may be as follows. In all of the
ordinary differential equations in the below, y is the dependent variable, and x is the
independent variable.
𝑑2𝑦 dy
1) + 5 + 6𝑦 = 0 Linear
d𝑥 2 dx
𝑑4𝑦 2
𝑑3𝑦 𝑑𝑦
2)
4
+ 𝑥 3
+ 𝑥3 = 𝑥𝑒 𝑥 Linear
d𝑥 𝑑𝑥 𝑑𝑥
𝑑2𝑦 𝑑𝑦
3) 𝑥 2
+𝑦 = 𝑒𝑥 Non-linear
𝑑𝑥 𝑑𝑥
P a g e |3
𝑑2 𝑦 dy 3
4) + 5 (dx ) + 6𝑦= 𝒙𝟐 Non-linear
𝑑𝑥 2
𝑑2 𝑦 dy
5)
2
+ 5 + 6𝑦 2 = Sinx Non-linear
𝑑𝑥 dx
𝑑2 𝑦 dy 2
1) + 𝑥𝑦 ( ) = 0 O=2, O.D.E, Dep= y, Ind= x
𝑑𝑥 2 dx
𝑑4𝑥 𝑑2𝑥
2) + 5 + 3𝑥 = sin 𝑡 O=4, O.D.E, Dep= x, Ind= t
𝑑𝑡 4 𝑑𝑡 2
𝜕𝑣 𝜕𝑣
3) + =𝑣 O=1, P.D.E., Dep= v, Ind= s, t
𝜕𝑠 𝜕𝑡
𝑑2𝑦 2 𝑑𝑦
1) ( ) + (1 + )3 = 0 O=2, O.D.E., D=3, Dep=y, Ind=x
𝑑𝑡 2 𝑑𝑡
𝑑 3 𝑦 𝑑 2 𝑦 𝑑𝑦
2) + + = 𝑦 2 𝑡𝑎𝑛 𝑡 O=3, O.D.E., D=2, Dep=y, Ind=t
𝑑𝑡 3 𝑑𝑡 2 𝑑𝑡
P a g e |4
Exercises:
Classify each of the following differential equations as:PDE,ODE, Linear, Non-Linear. Determine
the order, and degree of them;
𝑑𝑦
1) + 𝑥𝑦 = 𝑥𝑒 𝑥 ODE, Linear, 1st order, Degree=1
𝑑𝑥
𝑑3𝑦
2) + 𝑥√𝑦 = sin x ODE, non-Linear, 3rd order, Degree =1
d𝑡 3
𝜕2𝑢 𝜕2𝑢
3) + =0 PDE, Linear, 2nd order, Degree =1
𝜕𝑥 2 𝜕𝑦 2
𝑑2𝑦
4) + 𝑥 sin 𝑦 = 0 ODE, non-Linear, 2nd order
d𝑥 2
𝑑2𝑦
5) + 𝑦 sin 𝑥 = 0 ODE, Linear, 2nd order, Degree =1
d𝑥 2
separated. Therefore, these equations are called Separable Differential Equations. In separable
differential equations, integrations of both sides yields an implicit solution straightforward, as
follows;
∫ 𝑔(𝑦)𝑑𝑦 = ∫ 𝑓(𝑥)𝑑𝑥 + 𝑐.
Example 1.1
∫ 9𝑦𝑑𝑦 = ∫ −4𝑥𝑑𝑥
9 2
𝑦 = −2𝑥 2 + 𝐶1
2
𝑥2 𝑦2
+ = C.
9 4
This solution is an equation of an ellipse family for C > 0.
𝜕𝑢 𝜕𝑢
= 𝑀 and = 𝑁.
𝜕𝑥 𝜕𝑦
P a g e |6
In such a case, the left hand side of the differential equation can be simplified to an exact
differential
𝜕𝑢 𝜕𝑢
𝑑𝑢 = 𝑑𝑥 + 𝑑𝑦
𝜕𝑥 𝜕𝑦
and the equation is reduced to 𝑑𝑢 = 0. That is why the equation is called exact. An implicit
solution is simply 𝑢(𝑥, 𝑦) = 𝐶.
A sufficient condition for the exactness of a first order ODE is obtained by the observation
𝜕𝑀 𝜕2𝑢 𝜕𝑁 𝜕2𝑢
= , =
𝜕𝑦 𝜕𝑦𝜕𝑥 𝜕𝑥 𝜕𝑥𝜕𝑦
For twice differentiable functions
𝜕2𝑢 𝜕2𝑢
=
𝜕𝑦𝜕𝑥 𝜕𝑥𝜕𝑦
Hence a sufficient condition for a first order differential equation M(𝑥, 𝑦)𝑑𝑥 + N(𝑥, 𝑦)𝑑𝑦 = 0
to be exact is
𝜕𝑀 𝜕𝑁
=
𝜕𝑦 𝜕𝑥
𝜕𝑢
i. = 𝑀 → 𝑢 = ∫ 𝑀𝑑𝑥 + 𝑘(𝑦)
𝜕𝑥
𝜕𝑢
ii. Use the equation 𝜕𝑦 = N to obtain𝑑𝑘/𝑑𝑦
𝑑𝑘
iii. integrate 𝑑𝑦 to get k
Example 1.2
Show that the following equation is exact and find the solution.
𝜕𝑀
𝑀 = 𝑥 3 + 3𝑥𝑦 2 ⟹ = 6𝑥𝑦
𝜕𝑦
P a g e |7
𝜕𝑁
𝑁 = 3𝑥 2 𝑦 + 𝑦 3 ⟹ = 6𝑥𝑦
𝜕𝑥
Hence
𝜕𝑀 𝜕𝑁
= = 6𝑥𝑦
𝜕𝑦 𝜕𝑥
𝑥 4 3𝑥 2 𝑦 2
𝑢 = ∫ 𝑀𝑑𝑥 + 𝑘(𝑦) = ∫(𝑥 3 + 3𝑥𝑦 2 )𝑑𝑥 + 𝑘(𝑦) = + + 𝑘 (𝑦 )
4 2
𝜕𝑢 𝑑𝑘
= 3𝑥 2 𝑦 +
𝜕𝑦 𝑑𝑦
𝜕𝑢
= 𝑁 = 3𝑥 2 𝑦 + 𝑦 3
𝜕𝑦
Comparison gives
𝑑𝑘 𝑦4
= 𝑦3 ⟹ 𝑘 = +𝑐
𝑑𝑦 4
Hence
𝑦 4 3𝑥 2 𝑦 2 𝑦 4
𝑢(𝑥, 𝑦) = + + +𝑐
4 2 4
𝑦 4 3𝑥 2 𝑦 2 𝑦 4
+ + =𝐶
4 2 4
First order equations that are linear in y and y are called linear first order equations. The
most general for of this kind of differential equations is
y p ( x) y q( x) .
e y p ( x )e y q ( x )e
p ( x ) dx p ( x ) dx p ( x ) dx
.
e p ( x ) dx y q( x)e p ( x ) dx .
ye p ( x ) dx dx .
p ( x ) dx
q ( x )e
Example 1.3
4
y y 8x 3 .
x
4
p( x)dx x dx 4 ln x ln x
4
.
y e ln x
4
8x e dx
3 ln x 4 1
4
1
8 x 7 dx 4 x 8 C .
x x
P a g e |9
In this section, solutions of second order differential equations of different kinds are discussed.
Examples
Initial and Boundary Value Problems for Second Order Linear Differential Equations
𝑦 ′′ + 𝑎𝑦 ′ + 𝑏𝑦 = 0
P a g e | 10
𝑟 2 𝑒 𝑟𝑥 + 𝑎. r𝑒 𝑟𝑥 + b. 𝑒 𝑟𝑥 = 0 , (𝑟 2 + 𝑟𝑎 + 𝑏)𝑒 𝑟𝑥 = 0
𝑒 𝑟𝑥 ≠ 0, hence
𝑦 = 𝑒 𝑟𝑥 ⟹ 𝑟 2 + ar + b = 0
Example 1.4
𝑦 ′′ + 𝑦 = 0, y(0) = 3, 𝑦(𝜋/2) = −3
𝑦 = 𝑒 𝑟𝑥 ⟹ (𝑟 2 + 1)𝑒 𝑟𝑥 = 0 ⟹ 𝑟1,2 = ∓𝑖
𝑦 = −3 𝑠𝑖𝑛𝑥 + 3𝑐𝑜𝑠𝑥
Example 1.5
Euler-Cauchy Equations
Euler-Cauchy Equations are second order linear equations with variable coefficients of the
special form
𝑥 2 𝑦 ′′ + 𝑎𝑥𝑦 ′ + 𝑏𝑦 = 𝑟(𝑥)
𝑥 = 𝑒 𝑧 → 𝑧 = ln 𝑥
𝑑 1 𝑑 𝑑𝑦 𝑑𝑦
= → 𝑥 =
𝑑𝑥 𝑥 𝑑𝑧 𝑑𝑥 𝑑𝑧
and
𝑑𝑦 𝑑2𝑦 1 𝑑2𝑦 𝑑𝑦 2
𝑑2𝑦 𝑑2𝑦 2
𝑑2𝑦 𝑑 2 𝑦 𝑑𝑦
+𝑥 2 = → 𝑥 + 𝑥 = → 𝑥 = −
𝑑𝑥 𝑑𝑥 𝑥 𝑑𝑧 2 𝑑𝑥 𝑑𝑥 2 𝑑𝑧 2 𝑑𝑥 2 𝑑𝑧 2 𝑑𝑧
𝑑 2 𝑦 𝑑𝑦 𝑑𝑦 𝑧)
𝑑2𝑦 𝑑𝑦
2
− + 𝑎 + 𝑏𝑦 = 𝑟 ( 𝑒 → 2
+ (𝑎 − 1) + 𝑏𝑦 = 𝑟(𝑒 𝑧 )
𝑑𝑧 𝑑𝑧 𝑑𝑧 𝑑𝑧 𝑑𝑧
The resulting equation is a second order linear equation with constant coefficients.
Example 1.6
Then
𝑑𝑦 𝑑𝑦 𝑑2𝑦
2
𝑑 2 𝑦 𝑑𝑦
𝑥 = , and 𝑥 = −
𝑑𝑥 𝑑𝑧 𝑑𝑥 2 𝑑𝑧 2 𝑑𝑧
𝑑 2 𝑦 𝑑𝑦 𝑑𝑦 𝑑
− + 7 + 13𝑦 = 0 → 𝑦" + 6𝑦′ + 13𝑦 = 0, ′=
𝑑𝑧 2 𝑑𝑧 𝑑𝑧 𝑑𝑧
P a g e | 12
𝑚2 + 6𝑚 + 13 = 0 ⟹ 𝑚1,2 = −3 ∓ 2𝑖
𝑦 = 𝑥 −3 [𝐴𝐶𝑜𝑠(2𝑙𝑛𝑥) + 𝐵𝑆𝑖𝑛(2𝑙𝑛𝑥)]
𝑦(𝑥 ) = 𝑦ℎ (𝑥 ) + 𝑦𝑝 (𝑥)
𝑦 ′′ + 𝑝(𝑥)𝑦 ′ + 𝑞(𝑥)𝑦 = 0
In previous section we have seen how to find the general solution of a homogeneous second order
linear differential equation. Now we will elaborate methods to find a particular solution of the
non-homogeneous differential equation.
If the non-homogeneity function 𝒓(𝒙) is from a special kind of functions that create only a finite
number of root functions upon successive differentiations, it is called a function of finite
derivatives.
P a g e | 13
Example 1.7
1) 𝒓(𝒙) = 𝐒𝐢𝐧𝟐𝐱 is a function of finite derivatives since upon successive differentiations create
only two root functions;
D = {𝐱 𝟓 , x 4 , x 3 , x 2 , x, 1}
1 1 1
𝐃 = {ln 𝑥 , , 2 , 3 , … }
x x x
Now let the set of the two linearly independent solutions of homogeneous second order
linear differential equation is
1) If the sets H and D do not have any common function, then we propose a particular
solution for the non-homogeneous differential equation as a linear combination of functions
in D with coefficients {a1 , a2 , … , an } to be determined
Example 1.8
Solve 𝑦 ′′ − 3𝑦 ′ + 2𝑦 = 𝑒 3𝑥
𝑟 2 − 3r + 2 = 0 ⟹ 𝑟1,2 = 1,2
𝑦ℎ = 𝐶1 𝑒 𝑥 + 𝐶2 𝑒 2𝑥 ; 𝐻 = {𝑒 𝑥 , 𝑒 2𝑥 }
where
D = {e3x }, and H ∩ D = ∅
𝑦𝑝 = 𝐶𝑒 3𝑥
𝑦𝑝 ′ = 3𝐶𝑒 3𝑥 ⟹ 𝑦𝑝 ′′ = 9𝐶𝑒 3𝑥
2𝐶𝑒 3𝑥 = 𝑒 3𝑥 ⟹ 𝑐 = 1/2
Therefore
1 3𝑥
𝑦𝑝 = 𝑒
2
1
𝑦 = 𝐶1 𝑒 𝑥 + 𝐶2 𝑒 2𝑥 + 𝑒 3𝑥
2
2) If the sets H and D are not disjoint, we multiply the set D till we get
H ∩ xm D = ∅
Example 1.9
Solve 𝑦 ′′ − 3𝑦 ′ + 2𝑦 = 𝑒 𝑥
P a g e | 15
𝑟 2 − 3r + 2 = 0 ⟹ 𝑟1,2 = 1,2
𝑦ℎ = 𝐶1 𝑒 𝑥 + 𝐶2 𝑒 2𝑥 ; 𝐻 = {𝑒 𝑥 , 𝑒 2𝑥 }
where
𝑦𝑝 = 𝐶𝑥𝑒 𝑥 ⟹ 𝑦𝑝 ′ = 𝐶 (𝑒 𝑥 + 𝑥𝑒 𝑥 ) ⟹ 𝑦𝑝 ′′ = 𝐶 (2𝑒 𝑥 + 𝑥𝑒 𝑥 )
𝐶 (2𝑒 𝑥 + 𝑥𝑒 𝑥 ) − 3𝐶 (𝑒 𝑥 + 𝑥𝑒 𝑥 ) + 2𝐶𝑥𝑒 𝑥 = 𝑒 𝑥
−𝐶𝑒 𝑥 = 𝑒 𝑥 ⟹ 𝐶 = −1
Therefore
𝑦𝑝 = −𝑥𝑒 𝑥
𝑦 = 𝐶1 𝑒 𝑥 + 𝐶2 𝑒 2𝑥 − 𝑥𝑒 𝑥
Example 1.10
The roots of the characteristic equation, and the homogenous solution are;
𝑟 2 + 2r + 5 = 0 ⟹ 𝑟1,2 = −1 ∓ 2𝑖
H ∩ D1 = ∅ , H ∩ D2 = ∅
P a g e | 16
Therefore
𝑦𝑝 = 𝐶𝑒 𝑥 + 𝐾 cos 2𝑥 + 𝑀 sin 2𝑥
𝑦𝑝 ′ = 𝐶𝑒 𝑥 − 2𝐾 sin 2𝑥 + 2𝑀 cos 2𝑥
𝑦𝑝 ′′ = 𝐶𝑒 𝑥 − 4𝐾 cos 2𝑥 − 4𝑀 sin 2𝑥
𝑦 ′′ + 2𝑦 ′ + 5𝑦 = 16𝑒 𝑥 + sin 2𝑥 →
8𝐶 = 16, → 𝐶 = 2
cos 2𝑥 : − 4𝐾 + 4𝑀 + 5𝐾 = 0 → 𝐾 + 4𝑀 = 0
sin 2𝑥 : − 4𝑀 − 4𝐾 + 5𝑀 = 1 → −4𝐾 + 𝑀 = 1
4 1
𝐾=− ;𝑀=
17 17
Therefore
4 1
𝑦𝑝 = − cos 2𝑥 + sin 2𝑥
17 17
4 1
𝑦 = 𝑒 −𝑥 (𝐴 cos 2𝑥 + 𝐵 sin 2𝑥 )+2𝑒 𝑥 − cos 2𝑥 + sin 2𝑥
17 17
with functions 𝑝 , 𝑞 and 𝑟 which are continuous in the given interval, the method of variation of
parameters is used to find a particular solution 𝑦𝑝 .
𝑦2 𝑟 𝑦1 𝑟
𝑦𝑝 (𝑥 ) = −𝑦1 ∫ 𝑑𝑥 + 𝑦2 ∫ 𝑑𝑥
|𝑊 | |𝑊 |
𝑦 ′′ + 𝑝(𝑥)𝑦 ′ + 𝑞(𝑥)𝑦 = 0
corresponding to original equation, and
𝑦1 𝑦2
𝑊 = 𝑦1 𝑦2′ − 𝑦1 ′𝑦2 → |𝑊 | = |𝑦 ′ 𝑦2′ |
1
Note: The Wronskian determinant is not zero when 𝑦1 , 𝑦2 are linearly independent and solution
is feasible.
Proof:
Assume
𝑢′ 𝑦1 + 𝑣 ′ 𝑦2 = 0. (1)
P a g e | 18
𝑦𝑝 ′′ = 𝑢′ 𝑦1 ′ + 𝑢𝑦1 ′′ + 𝑣 ′ 𝑦2 ′ + 𝑣𝑦2 ′′
since 𝑦1 and 𝑦2 are solutions of the differential equation, the above equation reduces to
𝑢′ 𝑦1 ′ + 𝑣 ′ 𝑦2 ′ = 𝑟(𝑥) (2)
and combining equation (2), with equation (1), a linear algebraic system of two linear equations
for 𝑢′ , 𝑣 ′ is obtained.
𝑦1 𝑦2 𝑢′ 0
[ ][ ]
𝑦1 ′ 𝑦2 ′ 𝑣′ = [𝑟(𝑥)].
0 𝑦2 𝑦 0
| | | 1′ | 𝑦 𝑟
′ 𝑟 𝑦2′ 𝑦2 𝑟 ′ 𝑦1 𝑟 1
𝑢 = =− , 𝑣 = = .
|𝑊 | 𝑊 |𝑊 | |𝑊 |
𝑦2 𝑟 𝑦1 𝑟
𝑢 = −∫ 𝑑𝑥, 𝑣=∫ 𝑑𝑥.
|𝑊 | |𝑊 |
𝑦2 𝑟 𝑦1 𝑟
𝑦𝑝 (𝑥) = 𝑢(𝑥 )𝑦1 (𝑥 ) + 𝑣 (𝑥 )𝑦2 (𝑥 ) = −𝑦1 ∫ 𝑑𝑥 + 𝑦2 ∫ 𝑑𝑥.
𝑊 𝑊
Example 1.11
Solve 𝑦 ′′ + 𝑦 = sec 𝑥
The roots of the characteristic equation, and the two linearly independent solutions are;
𝑦ℎ = 𝐴 cos x + 𝐵 sin x
cos x sin x
𝑊 (y1 , y2 ) = | | = cos 2 x + sin2 x = 1
− sin x cos x
Example 1.12
𝑒 −2𝑥
Solve 𝑦 ′′ + 4y′ + 4𝑦 = 𝑥2
The roots of the characteristic equation, and the two linearly independent solutions are;
yh = C1 e−2x + C2 xe−2x
y2 r y1 r
yp (x) = −y1 ∫ dx + y2 ∫ dx
W W
y1 y2 e−2x xe−2x
W = |y ′ | = e−4x − 2xe−4x + 2xe−4x = e−4x
1 y2 ′| = |−2e−2x e −2x
− 2xe−2x
x( e−2x )
yp (x) = −e−2x ln x − = −e−2x ln x − e−2x
x
Example 1.13
Solve
2𝑒 2𝑥
𝑦 ′′ − 4y′ + 5𝑦 =
sin(𝑥)
The roots of the characteristic equation, and the two linearly independent solutions are;
Example 1.14
The roots of the characteristic equation, and the two linearly independent solutions are;
𝑟 2 − 4𝑟 + 4 = 0 ⇒ 𝑟1,2 = 2 ⟹ 𝑦1 = 𝑒 2𝑥 𝑦2 = 𝑥𝑒 2𝑥
𝑦ℎ = 𝐶1 𝑒 2𝑥 + 𝐶2 𝑥𝑒 2𝑥
𝑥𝑒 2𝑥 6𝑥 −4 𝑒 2𝑥 𝑒 2𝑥 6𝑥 −4 𝑒 2𝑥
𝑦𝑝 = −𝑒 2𝑥 ∫ 𝑑𝑥 + 𝑥𝑒 2𝑥 ∫
𝑑𝑥
𝑒 4𝑥 𝑒 4𝑥
= −𝑒 2𝑥 ∫ 6𝑥 −3 𝑑𝑥 + 𝑥𝑒 2𝑥 ∫ 6𝑥 −4 𝑑𝑥 3𝑒 2𝑥 𝑥 −2 − 2𝑒 2𝑥 𝑥 −2 = 𝑒 2𝑥 𝑥 −2 .
𝑦𝑔𝑒𝑛 = 𝑦ℎ + 𝑦𝑝 = 𝐶1 𝑒 2𝑥 + 𝐶2 𝑥𝑒 2𝑥 + 𝑒 2𝑥 𝑥 −2 .
P a g e | 22
1.4 Exercises
Find the solutions of the following ODE's
1. (𝑥 2 − 4)𝑦 ′ = 2𝑥𝑦 + 6𝑥
2. 3𝑥𝑦 ′ − 𝑦 = ln 𝑥 + 1, 𝑥 > 0, 𝑦(1) = −4
𝑡
3. 1 + (𝑦 − 𝑐𝑜𝑠 𝑦) 𝑦 ′ = 0
4. 𝑦 ′ = 𝑒 𝑙𝑛 𝑥−𝑦+2
5. 𝑦 ′ (𝑒 −𝑦 + 2𝑥) − 𝑒 𝑦 = 0, 𝑦(0) = 0
Answers: (Check)
𝐶
1. 𝑦 = 𝑥 2−4 − 3
2. 𝑦 = −(ln 𝑥 + 4)
3. t = sin y + y −1 cos y + Cy −1
e2 x2
4. y = ln (C + )
2
−𝑦 ) 1 1
5. 𝑥 = −2𝑒 −2(1+𝑒 − 2 𝑒 −𝑦 + 4