Module

for

Adams-Bashforth-Moulton Method for O.D.E.'s

   

    The methods of Euler, Heun, Taylor and Runge-Kutta are called single-step methods because they use only the information from one previous point to compute the successive point, that is, only the initial point  [Graphics:Images/AdamsBashforthMod_gr_1.gif]  is used to compute  [Graphics:Images/AdamsBashforthMod_gr_2.gif]  and in general  [Graphics:Images/AdamsBashforthMod_gr_3.gif]  is needed to compute  [Graphics:Images/AdamsBashforthMod_gr_4.gif].  After several points have been found it is feasible to use several prior points in the calculation.  The Adams-Bashforth-Moulton method uses [Graphics:Images/AdamsBashforthMod_gr_5.gif]  in the calculation of [Graphics:Images/AdamsBashforthMod_gr_6.gif].  This method is not self-starting;  four initial points  [Graphics:Images/AdamsBashforthMod_gr_7.gif], [Graphics:Images/AdamsBashforthMod_gr_8.gif], [Graphics:Images/AdamsBashforthMod_gr_9.gif],  and [Graphics:Images/AdamsBashforthMod_gr_10.gif] must be given in advance in order to generate the points [Graphics:Images/AdamsBashforthMod_gr_11.gif].  

    A desirable feature of a multistep method is that the local truncation error (L. T. E.) can be determined and a correction term can be included, which improves the accuracy of the answer at each step.  Also, it is possible to determine if the step size is small enough to obtain an accurate value for  [Graphics:Images/AdamsBashforthMod_gr_12.gif], yet large enough so that unnecessary and time-consuming calculations are eliminated.  If the code for the subroutine is fine-tuned, then the combination of a  predictor and corrector requires only two function evaluations of  f(t,y)  per step.  

 

Theorem  (Adams-Bashforth-MoultonMethod)  Assume that  f(t,y)  is continuous and satisfies a Lipschits condition in the variable  y,  and consider the  I. V. P. (initial value problem)

            
[Graphics:Images/AdamsBashforthMod_gr_13.gif] with [Graphics:Images/AdamsBashforthMod_gr_14.gif],  over the interval  [Graphics:Images/AdamsBashforthMod_gr_15.gif].
        
The Adams-Bashforth-Moulton method uses the formulas [Graphics:Images/AdamsBashforthMod_gr_16.gif],  and  

the predictor          [Graphics:Images/AdamsBashforthMod_gr_17.gif],  and

the corrector           [Graphics:Images/AdamsBashforthMod_gr_18.gif]  for  [Graphics:Images/AdamsBashforthMod_gr_19.gif]  

as an approximate solution to the differential equation using the discrete set of points  [Graphics:Images/AdamsBashforthMod_gr_20.gif].  

Remark.  The Adams-Bashforth-Moulton method is not a self-starting method.  Three additional starting values [Graphics:Images/AdamsBashforthMod_gr_21.gif]  must be given.  They are usually computed using the Runge-Kutta method.

Proof  Adams-Bashforth-Moulton Method  Adams-Bashforth-Moulton Method    

 

Theorem (Precision of Adams-Bashforth-MoultonMethod)  Assume that  [Graphics:Images/AdamsBashforthMod_gr_22.gif]  is the solution to the I.V.P.  [Graphics:Images/AdamsBashforthMod_gr_23.gif]  with  [Graphics:Images/AdamsBashforthMod_gr_24.gif].  If  [Graphics:Images/AdamsBashforthMod_gr_25.gif]  and  [Graphics:Images/AdamsBashforthMod_gr_26.gif]  is the sequence of approximations generated by Adams-Bashforth-Moulton method, then at each step, the local truncation error is of the order  [Graphics:Images/AdamsBashforthMod_gr_27.gif],  and the overall global truncation error  [Graphics:Images/AdamsBashforthMod_gr_28.gif] is of the order

        
[Graphics:Images/AdamsBashforthMod_gr_29.gif],  for  [Graphics:Images/AdamsBashforthMod_gr_30.gif].  


The error at the right end of the interval is called the final global error  

        
[Graphics:Images/AdamsBashforthMod_gr_31.gif].  

Proof  Adams-Bashforth-Moulton Method  Adams-Bashforth-Moulton Method    

 

Animations (Adams-Bashforth-Moulton Method  Adams-Bashforth-Moulton Method).  Internet hyperlinks to animations.

 

Algorithm (Adams-Bashforth-Moulton Method).  To approximate the solution of the initial value problem [Graphics:Images/AdamsBashforthMod_gr_32.gif] with [Graphics:Images/AdamsBashforthMod_gr_33.gif]  over  [Graphics:Images/AdamsBashforthMod_gr_34.gif]  at a discrete set of points using the formulas:

use the predictor      [Graphics:Images/AdamsBashforthMod_gr_35.gif]  

and  the corrector       [Graphics:Images/AdamsBashforthMod_gr_36.gif]  for  [Graphics:Images/AdamsBashforthMod_gr_37.gif].  

 

Computer Programs  Adams-Bashforth-Moulton Method  Adams-Bashforth-Moulton Method    

 

Mathematica Subroutine (Adams-Bashforth-Moulton Method).

[Graphics:Images/AdamsBashforthMod_gr_38.gif]

Example 1.  Solve the I.V.P.  [Graphics:Images/AdamsBashforthMod_gr_39.gif].  
Solution 1.

 

Example 2.  Use Mathematica to find the analytic solution and graph for the I.V.P.  [Graphics:Images/AdamsBashforthMod_gr_50.gif].  
Solution 2.

 

Example 3.  Plot the error for Adams-Bashforth-Moulton's method.
Solution 3.

 

Example 4.  Reduce the step size by  [Graphics:Images/AdamsBashforthMod_gr_70.gif] and see what happens to the error.
Recalculate points for Adams-Bashforth-Moulton's method, and the analytic solution using twice as many subintervals.
Then Plot the error for Adams-Bashforth-Moulton's method.
Solution 4.

 

Example 5.   Solve  [Graphics:Images/AdamsBashforthMod_gr_81.gif]  with  [Graphics:Images/AdamsBashforthMod_gr_82.gif]  over  [Graphics:Images/AdamsBashforthMod_gr_83.gif].
Solution 5.

 

Example 6.  Use Mathematica to find the analytic solution and graph for the I.V.P.  [Graphics:Images/AdamsBashforthMod_gr_91.gif].  
Solution 6.

 

Example 7.  Plot the absolute value of the error for Adams-Bashforth-Moulton's method.
Solution 7.

 

Example 8.  Reduce the step size by  [Graphics:Images/AdamsBashforthMod_gr_109.gif] and see what happens to the error.
Recalculate points for Adams-Bashforth-Moulton's method, and the analytic solution using twice as many subintervals.
Then Plot the error for Adams-Bashforth-Moulton's method.
Solution 8.

 

Example 9.  Solve the I.V.P.  [Graphics:Images/AdamsBashforthMod_gr_120.gif].  
Solution 9.

 

Example 10. Use Mathematica to find the analytic solution and graph for the I.V.P.  [Graphics:Images/AdamsBashforthMod_gr_145.gif].  
Solution 10.

 

Various Scenarios and Animations for the Adams-Bashforth-Moulton Method.

Example 11.  Solve the I.V.P.  [Graphics:Images/AdamsBashforthMod_gr_165.gif].    Compute the Adams-Bashforth-Moulton solution to the I.V.P.
Solution 11.

 

Animations (Adams-Bashforth-Moulton Method  Adams-Bashforth-Moulton Method).  Internet hyperlinks to animations.

 

Old Lab Project  (Adams-Bashforth-Moulton Method  Adams-Bashforth-Moulton Method).  Internet hyperlinks to an old lab project.  

 

Research Experience for Undergraduates

Adams-Bashforth-Moulton Method  Adams-Bashforth-Moulton Method  Internet hyperlinks to web sites and a bibliography of articles.

 

Download this Mathematica Notebook Adams-Bashforth-Moulton Method for O.D.E.'s

 

Return to Numerical Methods - Numerical Analysis

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004