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