CH#4 Numerica Differential Equations and Difference Equations
CH#4 Numerica Differential Equations and Difference Equations
CHAPTER # 4
CHAPTER # 4
PART-I: NUMERICAL DIFFERENTIAL EQUATIONS
1): Euler’s Method:-
Aim:-
……………………..(1)
Procedure:-
or
∫ ∫
……………… (2)
This is called Euler’s iteration formula for finding numerical solution of first order first
degree differential equation.
For , we have:-
For we have:-
For , we have:-
, and so on.
This Euler’s iteration formula (2) can be used to find , where is known. On
substituting the value of in (2) we get . Similarly putting the value of
in (2), we obtain and so on.
Note
… ………………………….. ……………………………………….
Example# 1:-
Using Euler’s method, find an approximate value of y corresponding to , given that:
Solution:-
We have
“Table of values”
2.0 9.55
Example# 2(HW):-
Evaluate by Euler’s method from differential equation with
.
Answer:-
……………………..(1)
Procedure:-
OR
∫ ∫
∫ ∫
………… (2)
This is called Modified Euler’s iteration formula for finding numerical solution of first
order first degree differential equation.
Since which occurs on the right hand of eq. (2) can not be calculated since
Example# 3 :-
Solution :-
we have
“Table of values”
10 1.0 21.081
Example# 4(HW) :-
is given by:
where,
( )
( )
This is known as Runge –Kutta iteration formula. The error in this formula is of order .
This method have greater accuracy. No derivatives are required to be tabulated. It requires
only functional values at some selected points on the sub-intervals.
Example #5 :-
Apply Runge – Kutta method to find an approximate value of when , given that;
Solution:-
we have,
we have
( )
( )
We have,
( )
( )
and
Example# 6(HW) :-
Apply Runge – Kutta method (fourth order), to find an approximate value of when
, given that:
Answer:-
Example# 7(HW) :-
Find
Answer:-
CHAPTER # 4
PART-II: DIFFERENCE EQUATIONS
Difference Equation:-
function.
The order of a difference equation is the difference between the largest and smallest
The degree of the highest differences function appearing in a difference equation free
from radicals and functions is called the degree of the difference equation.
Examples
1 1 1
2 2 1
3 3 1
4 2 2
5 3 1
6 4 3
7 2 3
8 5 4
Solution of a difference equation is an expression for which satisfies the given difference
equation.
The general solution of a difference equation is that in which the number of arbitrary
A particular solution or particular integral is that solution which is obtained from the
general solution by giving particular values to the constants (using IVP & BVP).
In linear difference equations, occur in the first degree only and are not
multiplied together.
OR
where,
Auxillary Equation :-
Replace by in and equate it to zero, we get
This is called auxiliary equation. Let be the roots of the Eq. The
following cases will arise:
given by
where
√ ( )
3±i4
√ ( )
2,3,3,4,4,4,5,6
2,3,3,5,6,3±i4
√ ( )
Example 1:-
Solve: – .
Solution:-
The given equation can be written as
Example 2:-
Solve: –
Solution:-
The given equation can be written as
Example 3:-
Solve:
Solution:-
√ √ √
⇒
⇒ √
⇒
where,
√
√ ( )
Example 4:-
Solve: –
Solution:-
⇒ – , –
⇒ –
⇒ –
⇒ – –
⇒ –
Example 5:-
Solve: – –
Solution:-
– –
– –
Practice Problems:-
– – –
where,
⇒ is a function of
where,
and
There is no general procedure to find the particular solution of a difference equation. The
method depends on the form of the function i.e. In some simple cases we give the rules
(Exponential Case):-
1): Case-I: :-
1):
2):
3):
4):
5):
and so on.
If, we have
– –
⇒ But
⇒ etc.
– and ( – )
Ponit-3: The Remaining factor w.r.t. root 2 is ( – ) and Remaining factor w.r.t. root 4
(actually it is -2) and the Remaining factor w.r.t. root 4 is ( – ), never be equal to
Note:-
This means that have the factor – and can be written as:
1):
2):
3):
4):
5):
6):
7):
Note:-
1):
2):
3):
Example 1:-
Solve: –
Solution:-
The given equation can be written as
–
For :-
The auxiliary equation is
–
⇒
Roots are distinct and real. Hence the homogeneous solution is
( ) ( )
For :-
We have
⇒ ( ) ( )
Example 2:-
Solve: –
Solution:-
The given equation can be written as
For :-
The auxiliary equation is
⇒
⇒
Roots are repeated and real. Hence the homogeneous solution is
For :-
We have
(Here )
(Quadratic Case)
Example 3:-
Solve:
Solution:-
The given equation can be written as
For :-
The auxiliary equation is
⇒
⇒
Roots are distinct and real. Hence the homogeneous solution is
For :-
We have
⇒ (Here )
We have
⇒ (Linear Case)
Note:-
If
⇒ ( ) and ( )
(Oscillating Case):-
Then
OR
can be written as
and solve as in Case-I with . At the end separate real and imaginary parts
We have
Hence in general
Example 1:-
–
Solution:-
The given equation can be written as
:-
⇒
The homogenous solution is
For :-
( )
⇒
⇒
( )
Example 2:-
Solution:-
:-
⇒ ⇒
⇒ √ ⇒ ( )
⇒ = ( ) 1
⇒ =
For :-
We have
( )
⇒
⇒ =
Example 3:-
1):
2):
( )
4):
Factorial Polynomial :-
Derivative Formula
Power Formula
Zero Formula
2):
⇒ ( ) [ ]
⇒ ( )
3):
⇒ [ ]
⇒
⇒
(Polynomial Case):-
We have
Expansion Formulae of :-
1):
2):
Example-1:-
Solution:-
The given equation can be written as
:-
⇒ ⇒
:-
⇒ * +
⇒ [ ( ) ]
⇒ * +( )
⇒ [ ]( )
⇒ [ ( ) ]
⇒ [ ]
⇒ [ ]
⇒ [ ]
Example-2:-
:
Solution:-
The given equation can be written as
:-
⇒ ⇒ ⇒
Roots are distinct and real. Hence the homogeneous solution is
:-
We have
⇒ ( ) ( )
⇒ ( )
⇒ * + ( )
⇒ * ( ) +
⇒ * +( )
⇒ * +( )
⇒ [ ( )
( )]
⇒ [ ( ) ]
⇒ * +
⇒ ( )
⇒ ( )
, :-
(Mixture/Combine Case)
We have
Example-1:-
Solution:-
The given equation can be written as
:-
⇒ ⇒
Roots are distinct and real. Hence the homogeneous solution is
We have
Replace by , we get
⇒ * +
⇒ [ ( ) ]
⇒ [ ]
⇒ [ ]
⇒ [ ( ) ]
⇒ [ ]
⇒ [ ]
⇒ * +
Example-2:-
–
Solution:-
:-
⇒ ⇒
We have
⇒ =
Replace by , we get
⇒ ( )
⇒ ( )
⇒ ( )
⇒ ( )
⇒ ( ) ( )
⇒ [ ( ) ]
Example-3:-
–
Solution:-
:-
We have
⇒ ( )
⇒ ( ) ( ) ( )
1):
2): –
– ( )
4): ( – )
⇒
⇒
⇒
⇒
EXERCISE-4
PART-I: NUMERICAL DIFFERENTIAL EQUATIONS
Question # 1:-
Answer:-
Question # 2:-
Using Euler’s method find an approximate value of corresponding to , given that
Answer:-
Question # 3:-
Apply modified Euler’s method to solve and hence find
an approximate value of
Answer:-
Question # 4:-
Using Euler’s modified formula, find an approximate value of
taking the interval 0.02.
Answer:-
Question # 5:-
Apply RK method to find an approximate value of 0.2, given that:
.
Answer:-
Question # 6:-
Apply RK method to find an approximate value of 0.2, given that:
.
Answer:-
Question # 7:-
Use Euler’s method to compute from the differential equation
when
Answer:-
Question # 8:-
Given that with initial condition that find by Euler’s
modified method , , correct up to three decimal Places.
Answer:-
Question # 9:-
Use Euler’s method to evaluate correct up to two decimal Places, from
.
Answer:-
Question # 10:-
Given that by
Euler’s method.
Answer:-
Question # 11:-
Evaluate by modified Euler’s method from differential equation with
.
Answer:-
Question # 12:-
Evaluate by modified Euler’s method from differential equation with
.
Answer:-
Question # 13:-
Question #15:-
Use RK-method to compute from .
Answer:-
Question # 16:-
Compute using RK-method
with step-size
Answer:-
Question # 17:-
Use RK-method to compute from .
Answer:-
Question # 18:-
Use RK-method to compute from .
Answer:-
Question # 19:-
Compute by RK-method correct to five decimal Places, from the equation
Answer:-
Question # 20:-
Find the values of using modified Euler’s method with , given
that:
Answer:-
Question # 21:-
Find the values of by using modified Euler’s method
given that:
Answer:-
Question # 22:-
Find the values of by using modified Euler’s method with
given that:
Answer:-
Question # 23:-
Find the values of by using Euler’s and modified Euler’s
methods respectively if:
Answer:-
Question # 24:-
Find
Answer:-
Question # 25:-
Find
Answer:-
Question # 26:-
Find
Answer:-
Question # 27:-
Find
Answer:-
Question # 28:-
Find
Answer:-
.
Question # 29:-
Find
Answer:-
Question # 30:-
Find
Answer:-
Question # 31:-
Find
Answer:-
Question # 32:-
Find
Answer:-
Question # 33:-
Find
.
Answer:-
Question # 34:-
The following ordinary differential equation(ODE) is used to describe a Population of Fish:
Where and are +ve constants; is known as the intrinsic growth constant. It measures
the differences between the birth and death rates per population unit in the absence of
overcrowding. The ratio is called saturation constant. Take and .
Use initial condition, Find , using Modified Euler’s Formula for step size
. Compare your results with the analytical solution of the above differential equation:
EXERCISE-4
PART-II: DIFFERENCE EQUATIONS
Solve the following difference equations
Question # 39:-
Answer:-
Question # 40:-
Answer:-
Question # 41:-
Answer:-
{ ( ) ( )}
Question # 42:-
,
Answer:-
Question # 43:-
, ,
Answer:-
Question # 44:-
Answer:-
√ √
( ) ( )
Question # 45:-
Answer:-
{ ( ) ( )}
Question # 46:-
Answer:-
Question # 47:-
Answer:-
Question # 48:-
Answer:-
Question # 49:-
Answer:-
Question # 50:-
Answer:-
(√ )
Question # 51:-
Answer:-
Question # 52:-
Answer:-
Question # 53:-
Answer:-
Question # 54:-
–
Answer:-
( ) ( )
Question # 55:-
Answer:-
Question # 56:-
Answer:-
Question # 57:-
Answer:-
Question # 58:-
Answer:-
( )
Question # 59:-
.
Answer:-
( )
Question # 60:-
–
Answer:-
( )
Question # 61:-
Answer:-
Question # 62:-
Answer:-
Question # 63:-
Answer:-
Question # 64:-
Answer:-
( √ ) ( √ )
Question # 65:-
Answer:-
Question # 66:-
Answer:-
Question # 67:-
Answer:-
Question # 68:-
.
Answer:-
Question # 69:-
Answer:-
( )
Question # 70:-
Answer:-
Question # 71:-
Answer:-
( )
Question # 72:-
.
Answer:-
Question # 73:-
Answer:-
( )
Question # 74:-
.
Answer:-
Question # 75:-
The Discrete Mathematical Modeling of Forced Vibrating System of particular oscillation
Phenomena leads to the following Nonhomogeneous Difference Equation of the form:
Question # 76:-
The Discrete Mathematical Modeling of Forced Vibrating System of particular oscillation
Phenomena leads to the following Nonhomogeneous Difference Equation of the form:
Answer:-
( )
Question # 77:-
The Discrete Mathematical Modeling of Forced Vibrating System of particular oscillation
Phenomena leads to the following Nonhomogeneous Difference Equation of the form:
Question # 78:-
The Discrete Mathematical Modeling of Some Vibrating System of particular oscillation
Phenomena leads to the following Nonhomogeneous Difference Equation of the form:
Answer:-
( )
SUMMARY OF CH#4
PART-I: NUMERICAL DIFFERENTIAL EQUATIONS
1): Euler’s Method:-
Aim:-
The Euler’s iteration formula for finding numerical solution of first order first degree
differential equation.
where,
The order of a difference equation is the difference between the largest and smallest
arguments appearing in it.
Degree of a Difference Equation:-
The degree of the highest differences function appearing in a difference equation free
from radicals and functions is called the degree of the difference equation.
Solution of Difference Equation:-
Solution of a difference equation is an expression for which satisfies the given difference
equation.
General Solution of Difference Equation:-
The general solution of a difference equation is that in which the number of arbitrary
constants is equal to the order of the difference equation.
Particular Integral/Solution of Difference Equation:-
A particular solution or particular integral is that solution which is obtained from the
general solution by giving particular values to the constants (using IVP & BVP).
Linear Difference Equation:-
In linear difference equations, occur in the first degree only and are not
multiplied together.
General Linear Difference Equation of Order :-
A general linear difference equation of order is of the form:
Step-2:-
Write auxiliary equation:
Auxillary Equation :-
Replace by in and equate it to zero, we get
Step-3:-
Write the solutions according to the nature of the roots of auxiliary equation.
Case-I (Distinct Real Roots)
Dr. Jamil Book Series 1 Page 51
Advanced Mathematical Techniques CH#4 Numerical Diff. Eqs. & Difference Equations
( )
where
√ ( )
where,
is the homogeneous solution o
and
= is the particular solution which satisfy,
:-
1): Case-I:
Note:-
If
⇒ ( ) and ( )
OR
can be written as
and solve as in Case-I with . At the end separate real and imaginary parts
according to the cosine and sine case.
Factorial Polynomial :-
Derivative Formula
Power Formula
Zero Formula
We have
Expansion Formulae of :-
1):
2):
,
We have