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.
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