JUST THE MATHS SLIDES NUMBER 17.7 NUMERICAL MATHEMATICS 7 - - PDF document

“JUST THE MATHS” SLIDES NUMBER 17.7 NUMERICAL MATHEMATICS 7 (Numerical solution)

• f

(ordinary differential equations (B)) by A.J.Hobson

17.7.1 Picard’s method

UNIT 17.7 NUMERICAL MATHEMATICS 7 NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS (B) 17.7.1 PICARD’S METHOD This method of solution is one of successive approxima-

• tion. It is an iterative method in which the numerical

results become more and more accurate, the more times it is used. An approximation for y in terms of x is substituted into the right hand side of the differential equation, dy dx = f(x, y). The equation is then integrated with respect to x, giving y in terms of x as a second approximation. Given numerical values are substituted into the second approximation and the result rounded off to an assigned number of decimal places or significant figures. The iterative process is continued until two consecutive numerical solutions are the same when rounded off to the required number of decimal places.

1

A hint on notation Consider the differential equation dy dx = 3x2, given that y = y0 = 7 when x = x0 = 2. This can be solved exactly to give y = x3 + C, which requires that 7 = 23 + C. Hence, y − 7 = x3 − 23;

• r, in more general terms,

y − y0 = x3 − x3

0.

Thus,

y

y0 dy =

x

x0 3x2dx.

That is,

x

x0

dy dxdx =

x

x0 3x2dx.

In future, we shall integrate both sides of the given dif- ferential equation with respect to x, from x0 to x.

2

EXAMPLES

• 1. Given that

dy dx = x + y2, and that y = 0 when x = 0, determine the value of y when x = 0.3, correct to four places of decimals. Solution

x

x0

dy dxdx =

x

x0(x + y2)dx,

where x0 = 0 Hence, y − y0 =

x

x0(x + y2)dx,

where y0 = 0. That is, y =

x

0 (x + y2)dx.

(a) First Iteration Replace y by y0 in the function to be integrated. y1 =

x

0 xdx = x2

2 ≃ 0.0450 at x = 0.3 (b) Second Iteration Now we use dy dx = x + y2

1 = x + x4

4 .

3

Therefore,

x

dy dxdx =

x    x + x4

4

    dx,

which gives y − 0 = x2 2 + x5 20. y2 = x2 2 + x5 20 ≃ 0.0451 at x = 0.3 (c) Third Iteration Now we use dy dx = x + y2

2

= x + x4 4 + x7 20 + x10 400. Therefore,

x

dy dxdx =

x    x + x4

4 + x7 20 + x10 400

    dx,

which gives y − 0 = x2 2 + x5 20 + x8 160 + x11 4400. y3 = x2 2 + x5 20 + x8 160 + x11 4400 ≃ 0.0451 at x = 0.3

4

• 2. If

dy dx = 2 − y x, and y = 2 when x = 1, perform three iterations of Picard’s method to estimate a value for y when x = 1.2. Work to four places of decimals throughout and state how accurate is the result of the third iteration. (a) First Iteration

x

x0

dy dxdx =

x

x0

 2 − y

x

  dx,

where x0 = 1. That is, y − y0 =

x

x0

 2 − y

x

  dx,

where y0 = 2. Hence, y − 2 =

x

1

 2 − y

x

  dx.

Replacing y by y0 = 2 in the function being integrated, we have y − 2 =

x

1

  2 − 2

x

   dx.

Therefore, y = 2 + [2x − 2 ln x]x

1

= 2 + 2x − 2 ln x − 2 + 2 ln 1 = 2(x − ln x). y1 = 2(x − ln x) ≃ 2.0354 when x = 1.2

5

(b) Second Iteration dy dx = 2 − y1 x = 2 − 2(x − ln x) x = 2 ln x x . Hence,

x

1

dy dxdx =

x

1

2lnx x dx. That is, y − 2 =

• (ln x)2

x

1 = (ln x)2.

y2 = 2 + (ln x)2 ≃ 2.0332 when x = 1.2 (c) Third Iteration Finally, we use dy dx = 2 − y2 x = 2 − 2 x − (ln x)2 x . Hence,

x

1

dy dxdx =

x

1

   2 − 2

x − (lnx)2 x

    dx.

That is, y − 2 =

   2x − 2 ln x − (ln x)3

3

   

x 1

. y − 2 = 2x − 2 ln x − (ln x)3 3 − 2.

6

y3 = 2x − 2 ln x − (ln x)3 3 ≃ 2.0293 when x = 1.2 The results of the last two iterations are identical when rounded off to two places of decimals, namely 2.03.

7