21.04.2015 Views

Lecture Notes in Differential Equations - Bruce E. Shapiro

Lecture Notes in Differential Equations - Bruce E. Shapiro

Lecture Notes in Differential Equations - Bruce E. Shapiro

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

361<br />

Figure 33.1: Illustration of Euler’s Method. A tangent l<strong>in</strong>e with slope<br />

f(t 0 , y 0 ) is constructed from (t 0 , y 0 ) forward a distance h = t 1 −t 0 <strong>in</strong> the t−<br />

direction to determ<strong>in</strong>ed y 1 . Then a l<strong>in</strong>e with slope f(t 1 , y 1 ) is constructed<br />

forward from (t 1 , y 1 ) to determ<strong>in</strong>e y 2 , and so forth. Only the first l<strong>in</strong>e is<br />

tangent to the actual solution; the subsequent l<strong>in</strong>es are only approximately<br />

tangent.<br />

y 2<br />

y 1<br />

y 0<br />

t 0 t 1 t 2<br />

S<strong>in</strong>ce y ′ (t) = f(t, y), we can approximate the left hand side of (33.10) by<br />

y ′ n(t n ) ≈ f(t n , y n ) (33.11)<br />

and hence<br />

y n+1 = y n + h n f(t n , y n ) (33.12)<br />

It is often the case that we use a fixed step size h = t j+1 − t j , <strong>in</strong> which case<br />

we have<br />

t j = t 0 + jh (33.13)<br />

In this case the Forward Euler’s method becomes<br />

y n+1 = y n + hf(t n , y n ) (33.14)<br />

The Forward Euler’s method is sometimes just called Euler’s Method.<br />

An alternate derivation of equation (33.12) is to expand the solution y(t)<br />

<strong>in</strong> a Taylor Series about the po<strong>in</strong>t t = t n :<br />

y(t n+1 ) = y(t n + h n ) = y(t n ) + h n y ′ (t n ) + h2 n<br />

2 y′′ (t n ) + · · · (33.15)<br />

= y(t n ) + h n f(t n , y( n )) + · · · (33.16)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!