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.

134 CHAPTER 6. LINEAR MULTISTEP METHODS<br />

of y n that is on both sides of the equation should be solved for using a Newton iteration<br />

as has been discussed previously in connection with the backward Euler method.<br />

✬<br />

Algorithm 6.1. Prediction-Correction Method To solve the<br />

initial value problem<br />

✫<br />

y ′ = f(t, y), y(t 0 ) = y 0<br />

for the function y(t), at each time step, let EM n and IM n represent<br />

predictions at time-step n using an explicit method (EM) and an<br />

implicit method (IM), respectively.<br />

P Predict y (0)<br />

n<br />

= EM n<br />

Repeat n times (n = 1 for PEC and PECE methods)<br />

E Estimate: f (0)<br />

n = f(t n , y (0)<br />

n )<br />

C Correct: y (1)<br />

n = IM(t n , y (0)<br />

n ).<br />

PEC methods stop here after the first iteration. P(EC) n methods<br />

stop here after n iterations.<br />

E : Estimate: f n (1) = f(t n , y n<br />

(1) ) to be used in the next iteration.<br />

PECE and P(EC) n E methods stop here.<br />

✩<br />

✪<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!