Example 6.  Compare the 2D Trapezoidal and 2D Simpson rule approximations to   [Graphics:Images/SimpsonsRule2DMod_gr_153.gif]  that were calculated in Examples 1 and 4?

Solution 6.

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



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

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

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

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

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

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

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

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

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

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


We can compare the error in these approximations.

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



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

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

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

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


The error for the Trapezoidal 2D rule has the form    [Graphics:../Images/SimpsonsRule2DMod_gr_170.gif]  where h and k are the step sizes for the variables x and y, respectively.
The error for Simpson's 2D rule has the form    [Graphics:../Images/SimpsonsRule2DMod_gr_171.gif]  where h and k are the step sizes for the variables x and y, respectively.
For the above examples, we have the following results.

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



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

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

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

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

Remark.  Since both the step sizes were reduced by a factor of [Graphics:../Images/SimpsonsRule2DMod_gr_177.gif] the remainder term  [Graphics:../Images/SimpsonsRule2DMod_gr_178.gif] should be reduced by approximately [Graphics:../Images/SimpsonsRule2DMod_gr_179.gif].  

Remark.  Since both the step sizes were reduced by a factor of [Graphics:../Images/SimpsonsRule2DMod_gr_180.gif] the remainder term  [Graphics:../Images/SimpsonsRule2DMod_gr_181.gif] should be reduced by approximately [Graphics:../Images/SimpsonsRule2DMod_gr_182.gif].  

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



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

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

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

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

Therefore, both the 2D Trapezoidal and 2D Simpson's rules is behaving as predicted.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004