Example 5.  Use the adaptive Simpson's rule to compute a numerical approximation to the integral  [Graphics:Images/AdaptiveQuadMod_gr_132.gif].  
Use the tolerances [Graphics:Images/AdaptiveQuadMod_gr_133.gif].  Compare with the analytic or "true value" of the integral.

Solution 5.

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

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

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


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

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


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

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


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

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


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

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



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

tol

0.001`

produces

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

tol

0.00001`

produces

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

tol

1.`*^-7

produces

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

true

value

is

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

 

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

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

[Graphics:../Images/AdaptiveQuadMod_gr_154.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_155.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_156.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004