Example 4.  Use cubic spline quadrature to compute a numerical approximation to the integral  [Graphics:Images/SplineQuadMod_gr_146.gif].  
Use the tolerances [Graphics:Images/SplineQuadMod_gr_147.gif].  Compare with the analytic or "true value" of the integral.

Solution 4.

4 (a). Plot the function over the interval  [0, 2].

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

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

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

4 (b). Construct the cubic spline for 11 nodes and use it for quadrature.

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

[Graphics:../Images/SplineQuadMod_gr_152.gif]
[Graphics:../Images/SplineQuadMod_gr_153.gif]

4 (c). Construct the cubic spline for 21 nodes and use it for quadrature.

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

[Graphics:../Images/SplineQuadMod_gr_155.gif]
[Graphics:../Images/SplineQuadMod_gr_156.gif]

4 (d). Construct the cubic spline for 41 nodes and use it for quadrature.

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

[Graphics:../Images/SplineQuadMod_gr_158.gif]
[Graphics:../Images/SplineQuadMod_gr_159.gif]

4 (e). Construct the cubic spline for 41 nodes and use it for quadrature.

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

[Graphics:../Images/SplineQuadMod_gr_161.gif]
[Graphics:../Images/SplineQuadMod_gr_162.gif]

4 (f). Compare the results from parts b-d.

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

m sample points

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

11

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

21

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

41

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

81

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

 

4 (g). Use Mathematica to find the analytic solution to the integral, i.e. the "true value" of the integral.

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


[Graphics:../Images/SplineQuadMod_gr_170.gif]
[Graphics:../Images/SplineQuadMod_gr_171.gif]
[Graphics:../Images/SplineQuadMod_gr_172.gif]
[Graphics:../Images/SplineQuadMod_gr_173.gif]
[Graphics:../Images/SplineQuadMod_gr_174.gif]
[Graphics:../Images/SplineQuadMod_gr_175.gif]
[Graphics:../Images/SplineQuadMod_gr_176.gif]
[Graphics:../Images/SplineQuadMod_gr_177.gif]

4 (h). How close did our last numerical approximation using Romberg integration come to the "true value" of the integral.

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


[Graphics:../Images/SplineQuadMod_gr_179.gif]
[Graphics:../Images/SplineQuadMod_gr_180.gif]


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

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

[Graphics:../Images/SplineQuadMod_gr_183.gif]
[Graphics:../Images/SplineQuadMod_gr_184.gif]
[Graphics:../Images/SplineQuadMod_gr_185.gif]
[Graphics:../Images/SplineQuadMod_gr_186.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004