Scribd
Upload
235 views

Numerical Methods II

This document discusses numerical methods for solving ordinary differential equations. It introduces the Taylor series method and Euler's method as single-step methods for solving initial value problems. Euler's method approximates the solution by advancing in discrete steps, while Taylor series expresses the solution as a power series. The document also covers modified Euler's method and multi-step Runge-Kutta methods of order two and four, which provide more accurate approximations than Euler's method by using slopes at multiple points within each step. Examples are provided to demonstrate how to apply each method to solve basic initial value problems numerically.

Uploaded by

Nik Hen
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
235 views

Numerical Methods II

This document discusses numerical methods for solving ordinary differential equations. It introduces the Taylor series method and Euler's method as single-step methods for solving initial value problems. Euler's method approximates the solution by advancing in discrete steps, while Taylor series expresses the solution as a power series. The document also covers modified Euler's method and multi-step Runge-Kutta methods of order two and four, which provide more accurate approximations than Euler's method by using slopes at multiple points within each step. Examples are provided to demonstrate how to apply each method to solve basic initial value problems numerically.

Uploaded by

Nik Hen
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 9

NUMERICAL METHODS II

Initial Value Problems for Ordinary Differential Equations


Introduction:
The analytic methods of solving differential equations are applicable only to
limited class of equations. The differential equations appearing in physical
problems do not fall into the category of familiar types. Therefore one has to
study numerical method to solve such equations.

Basically a first order

differential equation of the form

y = f(x,y), y(x0 ) = y0

.. (1)

is to be solved numerically by different methods.

The methods for the solution

of the initial value problem (1) can be classified mainly into two types. They are
single step method and multi-step methods. In single step methods, the solution
at any point Xi+1 is obtained using the solution at only the previous point

xi

where as in multistep method the solution is obtained using the solution at a


number of previous points. These methods yield the solution in one of the two
forms:
(i)

A series for y in terms of powers of x, from which the value of y can be


obtained by direct substitution.

(ii)

A set of tabulated values of x and y.

Taylor Series Method :


Here, we describe the method to solve the initial value problems using Taylor
series method. This method is the fundamental numerical method for the solution
of initial value problems (1).
Consider the differential equation (1)
Let y = y(x) be a continuously differentiable function satisfying the equation
(1). Expanding y in terms of Taylor series around the point x = x0, we get

(x-x0 )2 '' (x-x0 )3 '''


'
y = y0 + (x-x0 )y0 +
y0 +
y0 +........... (2)
2!
3!
The value of the differential coefficients y'0 , y''0 ,

y'''0 , .......... at x = x0

can be computed from the equation y' = f(x,y)


Example 1. Use Taylor series method to find y(0.1) from the equation =
3 + 2 with y(0)=1. Find the value of y for x = 0.1 correct to 3 decimal places.
Solution. Given = 3 + 2 .(i)
Differentiating w.r.t x, we get

= 3 + 2 .(ii)

= 2 + 2( )2 (iii)
Given y(0)=1. So putting x=0 and y=1 we get from (i) (0) = 1
Putting x=0 and y=1 and (0) = 1 from (ii) we get (0) = 3 + 2(1)(1) = 5
Similarly, we get (0) = 12 from (iii).
Substituting these values in the above Taylor series expansion (2),we get
=1++

2
2!

(5) +

3
3!

(12) +.. .

Therefore y(0.1)= (0.1) = 1 + 0.1 +

(0.1)2
2!

(5) +

(0.1)3
3!

(12) = 1.127 correct to

3 Decimal places.
Example 2.

Use Taylor series method to find y(0.1) and y (0.2) from the

equation = with y(0)=1.


Solution. Given = . Therefore (0) = 0.
= , (0) = 1
= 2 , (0) = 0
= 3 , (0) = 3 Similarly, (0) = 0, (0) = 15
Substituting these values in the above Taylor series expansion (2),we get
=1+0+
=1

2
2

2
2!

(1) +

4
8

6
48

3
3!

(0) +

4
4!

(3) +

5
5!

(0) +

6
6!

(15) +.. .

Exercises: Use Taylor series method to solve the differential equations


numerically correct to three decimal places.
i)

= 2 + 3 , y(0)=1 at x=0.1

ii)

= (1 2 2 ), y(0)=1 at x=0.1 and x=0.2

Eulers Method:
In Taylor series method, we express a series for y in terms of powers of x,
from which the value of y can be obtained by direct substitution. As the
approximation is poor, we derive Euler method by using Taylor series method,
where the values of y are computed by short steps ahead for equal intervals
h of the independent variable.
Consider the differential equation
with the initial condition
If

dy
f (x, y)
dx

y(x 0 ) y0 .

y(x) is the exact solution of the above equation, then the Taylors

series for y(x) around x x 0 is given by


y(x) y0 (x x 0 )y01

(x x 0 )2 11
y0 . . .
2!

Neglecting second and higher order terms, we get


y(x) y0 (x x 0 )y0' .

Denoting

x x 0 h , we get y(x 0 h) y0 h y(x 0 )' where

y(x 0 )' f (x 0 , y0 ) .

Taking x1 x 0 h , we get y1 y0 h f (x 0 , y0 ) .
Similarly

y2 y1 h f (x1 , y1 ) .

Proceeding in this way , for x n 1 x n h , we obtain the general formula


yn 1 yn h f (x n , yn ) .

Example 1. Using Eulers method, find an approximate value of y


corresponding to x=1, given that

= + 2 and y=1 when x=2.

Solution. (, ) = + 2 .
+1 = + ( , ) = + 0.1( + 2) Let us take n=10 and h=0.1.
The various calculations are arranged as follows.

+ 2 =

old y 0.1(dy dx) new y

1+0.1(3)=1.3

1.1 1.3

3.7

1.3+0.1(3.7)=1.67

1.2 1.67

4.54

1.67+0.1(4.54)=2.12

1.3 2.12

5.54

2.67

1.4 2.67

6.74

3.34

1.5 3.34

8.18

4.16

1.6 4.16

9.92

5.15

1.7 5.15

12

6.35

1.8 6.35

14.5

7.80

1.9 7.8

17.5

9.55

2.0 9.55

Thus the required approximate value of y(2)=9.55


Remark :
The process is very slow and to obtain reasonable accuracy with Eulers
method, we need to take a smaller value for h. Because of this restriction
on h, the method is unsuitable for practical use a modification of it, known
as modified Eulers method, which gives more accurate results.

Modified Eulers method :


Instead of approximating f (x, y) by f (x 0 , y0 ) , we approximate the integral
h
by means of Trapezoidal rule to obtain y1 y0 f (x 0 , y0 ) f (x1, y1(0) ) .
2

We thus obtain the iteration formula


y1(n 1) y0

h
f (x 0 , y 0 ) f (x1, y1(n) ) , n 0,1, 2,...
2

Where y1(n ) is the nth approximation to y1 , which can be found using


Eulers formula

y1(0) y0 hf (x 0 , y0 ) .

Example 1.Using modified Eulers method , Solve

= + ,y(0)=1

Choose h= 0.2, Compute y(0.4).


Solution : Here (, ) = + x0=0, y0=1. Let x1=x0+h = 0.2 and
x2=x1+h = 0.4. We have to find y1=y(x1)=y(0.2) and y2=y(x2)=y(0.4).
Now, f(x0,y0)=f(0.1)=0-1= -1.
Therefore by Eulers formula, y1(0)=y0+h f(x0,y0)= 1 + (0.2)(0+1)=1.2
(1)

By Modified Eulers formula, 1

(0)

= 0 + {(0 , 0 ) + (1 , 1 )}
2

=1+

0.2
2

{(0 + 1) + 0.2 + 1.2} = 1.2295

Next, again by Eulers formula,


y2(0)=y1+h f(x1,y1)= 1.2295 + (0.2)(0+1.2295)=1.4913
By Modified Eulers formula,

(1)

(0)

= 1 + {(1 , 1 ) + (2 , 2 )}

= 1.2295 +

0.2
2

{(0.2 + 1.2295) + 0.4 + 1.4913} = 1.5225

Note: In Euler method, the interval length h should be kept small and
hence these methods can be applied for tabulating y over a limited range
only.

Exercise
1. Using Eulers method, find the approximate value of y, when x=0.2 by
taking step length h=0.1. Given that
2. Solve the equation

= (1 + ),

= 1 ,

(0) = 1

(1) = 1 at x = 1.1 taking step

length h=0.05.
Runge-Kutta Methods
Introduction:
Eulers method is less efficient in practical problems since it requires h 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 value at some selected points on the subinterval. The basic idea of RK methods is to approximate the integral by a weighted average of slopes and
approximate slopes at a number of points in the interval [ xi, x i+1]
I. Runge Kutta Method of Second Order:

Consider the initial value problem

= (, ), ( ) = .

The second order R-K method formula is given by,

That is the value of y at x = xi is : = + ( + ).

Where,

= ( , )
= ( + , + ).

Example 1. Apply RK method of order two to find the value of y when x=0.2,
given that

= + y=1 when x=0.

Solution. Here (, ) = +

x0=0, y0=1 , h=0.2, f(x0,y0)=1

k1= hf(x0,y0)=0.2
k 2 = hf(x0 + h, y0 + k1 ) = 0.2(0.2,1.2) = 0.2(1.4) = 0.28
1

Then y(0.2) = y1 = y0 + (k1 + k 2 ) = 1 + (0.2 + 0.28) = 1.24


II.Runge-Kutta Method of fourth order:
The most commonly used RK method is a method which uses four slopes and is
called R-K method of fourth order. The method is given by :
Consider the Ordinary Differential Equation

= (, ), ( ) = .

= + . The fourth-order Runge-Kutta

We need to find y at
method formula is given by

= ( + ) = + ( + + + )

Where

= ( , )

= ( + , +
= ( + , +

)
)

= ( + , + ).
Example 1. Using 4th order Runge Kutta method , solve = + 2 with
y(0)=1 at x= 0.2 in steps of length h=0.1
Solution. Here (, ) = + 2

x0=0, y0=1 , h=0.1, x1= 0+h=0.1,

x2 = 0.1+0.1=0.2 To fond y2=y(x2)=y(0.2).


k1 = h f(x0, y0)=(0.1)(0+12)=0.1
k2 = h f(x0+h/2, y0+k1/2)=(0.1)f(0.05, 1.05)=0.11525
k3 = h f(x0+h/2, y0+k2/2)=(0.1)f(0.05, 1.057625)=0.116857
k4 = h f(x0+h, y0+k3)=(0.1)f(0.1, 1.116857)=0.13474. Then using RK formula,

= ( + ) = + ( + + + ) we get

y1 = y(0.1)= 1+1/6(0.1 +2 (0.11525 +2(0.116857)+0.13474) )=0.11649


k1 = h f(x1, y1)=(0.1)f(0.1,1.1165)=0.134657

k2 = h f(x0+h/2, y0+k1/2)=(0.1)f(0.15, 1.1838)=0.15514


k3 = h f(x0+h/2, y0+k2/2)=(0.1)f(0.15, 1.1941)=0.15759
k4 = h f(x0+h, y0+k3)=(0.1)f(0.2, 1.27409)=0.18233. Then using RK formula,
1

y2 = y(x1 + h) = y1 + (k1 + 2k 2 + 2k 3 + k 4 ) we get


6

y2 = y(0.2)= 1.1165+1/6(0.134657 +2 (0.15514 +2(0.15759)+0.118233) )


=1.27358.
Example 2. Using 4th order Runge Kutta method , solve = 1 + 2 with
y(0)=0 at x= 0.4 in steps of length h=0.2
Solution. Here (, ) = 1 + 2

x0=0, y0=0 , h=0.2, x1= 0+h=0.2,

x2 = 0.4 To fond y2=y(x2)=y(0.4).


k1 = h f(x0, y0)=(0.2)f(0,0)=0.2
k2 = h f(x0+h/2, y0+k1/2)=h{1+(y0+k1/2)}2 =(0.2){1+(0.1)2}=0.202
k3 = h f(x0+h/2, y0+k2/2)= (0.2){1+(0.101)2}=0.2020402
k4 = h f(x0+h, y0+k3)= (0.2){1+(0.2020402)2}=0.208164. Then using RK
1

formula, y1 = y(x0 + h) = y0 + (k1 + 2k 2 + 2k 3 + k 4 ) we get


6

y1 = y(0.2)= 1+1/6(0.2 +2 (0.202 +2(0.2020402)+0.208164) )=0.2027074


Again,
k1 = h f(x1, y1)=h(1+12 )=(0.2)(1+(0.2027074)2)=0.208218
k2 = h f(x0+h/2, y0+k1/2)=(0.2) (1+(0.3068164)2)=0.2188272
k3 = h f(x0+h/2, y0+k2/2)=(0.2) (1+(0.312121)2)=0.21948319
k4 = h f(x0+h, y0+k3)=(0.1) (1+(0.4221913)2)=0.235649.
1

Then using RK formula, y2 = y(x1 + h) = y1 + (k1 + 2k 2 + 2k 3 + k 4 )


6

we get
y2 = y(0.4)= 0.2027074+1/6(0.208218 +2 (0.2188272 +2(0.21948319) +
0.235649) ) =0.42279.

Exercise :

1) Using R.K Method of order 4, evaluate y(0.1).


Given that

= , (0) = 2

2) Using R.K Method of order 4, evaluate y(0.1) and y(0.2).


Given that

= 3 + 2, (0) = 0

You might also like