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

Solution 1.

[Graphics:../Images/AdaptiveQuadMod_gr_5.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_6.gif]

1 (a). Plot the function over the interval  [0, 1.25].

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

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

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


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

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


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

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


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

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



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


[Graphics:../Images/AdaptiveQuadMod_gr_17.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_18.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_19.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_20.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_21.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_22.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_23.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_24.gif]



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

tol

0.001`

produces

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

tol

0.00001`

produces

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

tol

1.`*^-7

produces

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

true

value

is

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

 

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

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

[Graphics:../Images/AdaptiveQuadMod_gr_34.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_35.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_36.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004