Example 1.  Use Adaptive Simpson's rule to compute a numerical approximation to the integral  [Graphics:Images/AdaptQuadMod_gr_3.gif].  
Use the tolerances [Graphics:Images/AdaptQuadMod_gr_4.gif].  Compare with the analytic or "true value" of the integral.

Solution 1.

[Graphics:../Images/AdaptQuadMod_gr_5.gif]
[Graphics:../Images/AdaptQuadMod_gr_6.gif]
[Graphics:../Images/AdaptQuadMod_gr_7.gif] [Graphics:../Images/AdaptQuadMod_gr_8.gif]
[Graphics:../Images/AdaptQuadMod_gr_9.gif] [Graphics:../Images/AdaptQuadMod_gr_10.gif]
[Graphics:../Images/AdaptQuadMod_gr_11.gif] [Graphics:../Images/AdaptQuadMod_gr_12.gif]
[Graphics:../Images/AdaptQuadMod_gr_13.gif] [Graphics:../Images/AdaptQuadMod_gr_14.gif]
[Graphics:../Images/AdaptQuadMod_gr_15.gif] [Graphics:../Images/AdaptQuadMod_gr_16.gif]
[Graphics:../Images/AdaptQuadMod_gr_17.gif]

tol

0.001`

produces

[Graphics:../Images/AdaptQuadMod_gr_18.gif]

tol

0.00001`

produces

[Graphics:../Images/AdaptQuadMod_gr_19.gif]

tol

1.`*^-7

produces

[Graphics:../Images/AdaptQuadMod_gr_20.gif]

true

value

is

[Graphics:../Images/AdaptQuadMod_gr_21.gif]


Did the adapt subroutine behave as expected ?

The following graph illustrates the example.

[Graphics:../Images/AdaptQuadMod_gr_22.gif]

[Graphics:../Images/AdaptQuadMod_gr_23.gif]

[Graphics:../Images/AdaptQuadMod_gr_24.gif]

[Graphics:../Images/AdaptQuadMod_gr_25.gif]

[Graphics:../Images/AdaptQuadMod_gr_26.gif]

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2003