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

Solution 6.

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

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

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


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

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


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

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


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

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


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

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



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

tol

0.001`

produces

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

tol

0.00001`

produces

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

tol

1.`*^-7

produces

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

true

value

is

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

 

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

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

[Graphics:../Images/AdaptiveQuadMod_gr_179.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_180.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_181.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004