21.04.2015 Views

The Computable Differential Equation Lecture ... - Bruce E. Shapiro

The Computable Differential Equation Lecture ... - Bruce E. Shapiro

The Computable Differential Equation Lecture ... - 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.

128 CHAPTER 6. LINEAR MULTISTEP METHODS<br />

Figure 6.1: Outer border of regions of stability for Adams Methods. Left: Adams-<br />

Bashforth methods AB1 (forward Euler, solid black); AB2 (Thick solid red); AB3<br />

(Thick dashed purple); and AB4 (thick, short blue dashed). Right: Adams-Moulton<br />

Methods AM3 (thin, blue); AM4 (thin, long dashed purple); AM5 (short-dashed<br />

green); AM6 (thick, solid, orange).<br />

Solution. From equation 6.104<br />

hf n =<br />

2∑<br />

k=1<br />

1<br />

k ∇k y k (6.106)<br />

= ∇y n + 1 2 ∇2 y n (6.107)<br />

= y n − y n−1 + 1 2 (y n − 2y n−1 + y n−2 ) (6.108)<br />

= 3 2 y n − 2y n−1 + 1 2 y n−2 (6.109)<br />

2hf n = 3y n − 4y n−1 + y n−2 (6.110)<br />

6.6 Nyström and Milne Methods<br />

<strong>The</strong>se methods consider the integral<br />

y(t n ) = y(t n−2 ) +<br />

∫ tn<br />

f(t, y(t))dt<br />

t n−2<br />

(6.111)<br />

Substitution of Newton’s backward difference formula yields Nyström’s Method<br />

∑ν−1<br />

y n = y n−2 + h k j ∇ j f n−1 (6.112)<br />

j=0<br />

Math 582B, Spring 2007<br />

California State University Northridge<br />

c○2007, B.E.<strong>Shapiro</strong><br />

Last revised: May 23, 2007

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

Saved successfully!

Ooh no, something went wrong!