Chapter 5 ODE Ashu

Chapter 5
Numerical Methods for Ordinary Differential Equations

INTRODUCTION
Numerical methods are becoming more and more important in engineering applications, simply
because of the difficulties encountered in finding exact analytical solutions but also, because of
the ease with which numerical techniques can be used in conjunction with modern high-speed
digital computers. Several numerical procedures for solving initial value problems involving
first-order ordinary differential equations are discussed in this chapter.
Definition: Solutions (particular) to differential equations may require satisfied certain defined
conditions and such conditions are called initial conditions if they are given at only one point of
the independent variable.
 The problem of solving an order ordinary differential equation together with initial
conditions is called an initial value problem (IVP).
 A first order initial value problem can be written as:
( ) ( ), ( ) …….. (1)
The single step methods for the solution of the IVP in (1):
1. Taylor series method
The Taylor series method is the fundamental numerical method for the solution of the IVP given
in (1).
Expanding ( ) in Taylor series about a point ,
( ) ( )( ( )
( ) ( ) ( )( ) ( ) ) …….. (2)
We denote the numerical solution and the exact solution at by and ( ), respectively.
Now, consider the interval [ ]. Substituting the length of the interval and
in (2), we obtain
( ) ( ) ( ) ( ) ( )( )
Neglecting the error term, we obtain the Taylor method as:

( )
……... (3)
is called an explicit single step Taylor series method of order .
1
Examples: Solve the following IVPs by using Taylor series method of 1st, 2nd and 3rd order.
a. ( ) in with .
b. ( ), ( ) in with .
c. , ( ) at and .
d. , [ ] with and when
Exercise:
i. Use the 2nd-order Taylor series method on ( ) for the IVP, ( )
Take h = 0.1.
ii. Use the fourth order Taylor series method with a single integration step to determine
( ) given that ( ) .
2. Euler’s Method
Euler’s explicit method (also called the forward Euler method) is a single-step, explicit method
for solving a first-order ordinary differential equation.
Let the differential equation is given by:
( ) with ( ) ( )
To solve the equation given in (4) for the values of at

integrate equation ( ) with one step size. ∫ ( ) (5)
Take ( ) ( ) in the range to , then we get
( ) ( ) ( )
Similarly, for the range to ( ) ( ), then
( ) ( ) ( )
In general form we have

( ) where
This is called Euler’s method.
Example: Using Euler’s method solve for of given that
2
and .
Exercise: Use Euler’s method to solve the following differential equation of
i. , ( ) and with step size .
ii. ( ) and with step size .
3. Modified Euler’s method

In case of Euler’s method we approximate ( ) by ( ) in the range to , but in this
case we approximate ( ) by means of Trapezoidal rule i.e.
( ) ( )
( ( ) ( ))
( )
Then, ( ( ) ( ))
Then, we get the iteration formula as:

( ) ( ) ))
( ( ) (
( )
and is the approximation of .
( )
The formula can be starting from by Euler’s formula.
( )
( ).
( ) ( )
* ( ) ( )+ is the ( ) a approximation of .
Example: Using Modified Euler’s method solves the following IVPs correct to 2 decimal places.
i. ( ) and at .
ii. given that and at .
iii. with ( ) and .
iv. ( ) and with step size .
4. Runge-Kutta Methods
Euler’s method is less efficient in practical problems since it requires to be small for obtaining
reasonable accuracy. The Runge-Kutta methods are designed to give greater accuracy and they
possess the advantage of requiring only the function values at some selected points in the
subinterval.
i. First-order Runge-Kutta method
The Euler’s method is called First-order Runge-Kutta method
3
( )
ii. Second- order Runge-Kutta method
The Euler’s modified method is called Second-order Runge-Kutta method
( )
where, ( ) and ( )
iii. Third- order Runge-Kutta method
( )
where, ( ), ( ) and ( ) such that
( )
iv. Fourth- order Runge-Kutta method: is the most commonly used.
( )
where, ( ), ( ), ( ) and
( )
Example:
i. Use the second order Runge-Kutta method with find and for
( )
ii. Use Runge-Kutta method of order two to integrate ( ) with ( ) from
to 0.5 in steps of . Keep four decimal places in the calculations.

iii. Use the Runge-Kutta method of order four with to obtain an approximation to
( ) for the solution of , ( ) . The exact solution is given by
. Determine the relative error and the percentage relative error.

Solution:
Exercise: solve , with correct to four decimal places and find ( )
and ( ) using Runge-Kutta method of
a. First order c. Third order
b. Second order d. Fourth order

Chapter 5 ODE Ashu

Uploaded by

Copyright:

Available Formats

Chapter 5 ODE Ashu

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Chapter 5 ODE Ashu

Uploaded by

Copyright:

Available Formats

Chapter 5

Numerical Methods for Ordinary Differential Equations

Neglecting the error term, we obtain the Taylor method as:

is called an explicit single step Taylor series method of order .

d. , [ ] with and when

i. Use the 2nd-order Taylor series method on ( ) for the IVP, ( )

To solve the equation given in (4) for the values of at

Take ( ) ( ) in the range to , then we get

Similarly, for the range to ( ) ( ), then

In general form we have

Example: Using Euler’s method solve for of given that

Exercise: Use Euler’s method to solve the following differential equation of

i. , ( ) and with step size .

ii. ( ) and with step size .

3. Modified Euler’s method

Then, we get the iteration formula as:

ii. given that and at .

iii. with ( ) and .

iv. ( ) and with step size .

where, ( ), ( ) and ( ) such that

ii. Use Runge-Kutta method of order two to integrate ( ) with ( ) from

to 0.5 in steps of . Keep four decimal places in the calculations.

( ) for the solution of , ( ) . The exact solution is given by

. Determine the relative error and the percentage relative error.

You might also like