Numerical Solution of Ordinary Differential Equations
Numerical Solution of Ordinary Differential Equations
Numerical Solution of Ordinary Differential Equations
OF ORDINARY
DIFFERENTIAL
EQUATIONS
TAYLOR SERIES EXPANSION
EULER / MODIFIED EULER EXPANSION
RUNGE-KUTTA METHODS
NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL
EQUATIONS
Ordinary Differential Equations is reputed to be one
of the most beautiful parts of mathematics.
It is also an essential tool for modeling many
physical solutions:
Electrical Networks
Mechanical Networks
Bending Beams
Chemical Reactions
Pendulums
The Motion of Rotating Mass Around Another
Body
NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL
EQUATIONS
Advantage:
It can solve practically any Differential
Equation; provided the appropriate
method/technique is used.
Disadvantage:
We must recompute the entire table if the
initial conditions are changed!
TAYLOR SERIES EXPANSION
PROBLEM:
dy
y' F ( x , y )
dx
PROCEDURE:
Take the problem & differentiate
repeatedly until it will equal to zero or
until you reach the specified degree
or repetition.
Apply results to the solution for the
nth-order Taylors Series Expansion.
TAYLOR SERIES EXPANSION
SOLUTION:
h2 h3
yi 1 yi hF ( xi , yi ) F' ( xi , yi ) F" ( xi , yi )
2! 3!
h n n 1
... F ( xi , yi )
n!
h2 h3 hn
yi hyi ' yi ' ' yi ' ' ' ...
n'
yi 1 yi
2! 3! n!
PROBLEM:
dy
y' F ( x , y )
dx
EULER / MODIFIED EULER
METHOD
Disadvantage:
It is much less accurate compared to the Taylor
Series Expansion; however, by using small increments
or step sizes, the result is acceptable enough.
SOLUTION:
PROBLEM:
dy
y' F ( x , y )
dx
SECOND-ORDER
RUNGE-KUTTA METHOD
yi 1 yi k1 k 2 yi k1 k 2
1 1 1
2 2 2
k1 hF ( xi , yi )
k 2 hF ( xi h , yi k1 )
THIRD-ORDER
RUNGE-KUTTA METHOD
yi 1 yi k1 4 k 2 k3
1
6
k1 hF ( xi , yi )
h k1
k 2 hF ( xi , yi )
2 2
k3 hF ( xi h , yi k1 2k 2 )
FOURTH-ORDER
RUNGE-KUTTA METHOD
yi 1 yi k1 2k 2 2k3 k4
1
6
k1 hF ( xi , yi )
h k1
k 2 hF ( xi , yi )
2 2
h k2
k3 hF ( xi , yi )
2 2
k4 hF ( xi h , yi k3 )