Example 4.  Use the adaptive Simpson's rule to compute a numerical approximation to the integral  [Graphics:Images/AdaptiveQuadMod_gr_91.gif]  that we investigated in example 2.  
The long solution is obtained if you add a print statement to investigate the in between computations.  
This subroutine is pedagogical and is intended to help us understand what's happening in a recursive program.
You would probably not want to always print out the in between steps, so you might want to re-execute the first version for some of your work.

Solution 4.

Now investigate the "recursive path" the subroutine takes.

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



[Graphics:../Images/AdaptiveQuadMod_gr_93.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_94.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_95.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_96.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_97.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_98.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_99.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_100.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_101.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_102.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_103.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_104.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_105.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_106.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_107.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_108.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_109.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_110.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_111.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_112.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_113.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_114.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_115.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_116.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_117.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_118.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_119.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_120.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_121.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_122.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_123.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_124.gif]


Now plot the intervals  [Graphics:../Images/AdaptiveQuadMod_gr_125.gif]  where Simpson Rule was applied in the form  

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

 

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

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

[Graphics:../Images/AdaptiveQuadMod_gr_129.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_130.gif]
[Graphics:../Images/AdaptiveQuadMod_gr_131.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004