Example 1  Investigate cubic spline quadrature for approximating the integral  [Graphics:Images/SplineQuadMod_gr_25.gif].  
Use  11, 21, 41 and 81 nodes.  Compare with the analytic or "true value" of the integral.

Solution 1.

[Graphics:../Images/SplineQuadMod_gr_26.gif]
[Graphics:../Images/SplineQuadMod_gr_27.gif]

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

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

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

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

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

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

[Graphics:../Images/SplineQuadMod_gr_32.gif]
[Graphics:../Images/SplineQuadMod_gr_33.gif]

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

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

[Graphics:../Images/SplineQuadMod_gr_35.gif]
[Graphics:../Images/SplineQuadMod_gr_36.gif]

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

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

[Graphics:../Images/SplineQuadMod_gr_38.gif]
[Graphics:../Images/SplineQuadMod_gr_39.gif]

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

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

[Graphics:../Images/SplineQuadMod_gr_41.gif]
[Graphics:../Images/SplineQuadMod_gr_42.gif]

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

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

m sample points

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

11

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

21

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

41

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

81

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

 

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

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


[Graphics:../Images/SplineQuadMod_gr_50.gif]
[Graphics:../Images/SplineQuadMod_gr_51.gif]
[Graphics:../Images/SplineQuadMod_gr_52.gif]
[Graphics:../Images/SplineQuadMod_gr_53.gif]
[Graphics:../Images/SplineQuadMod_gr_54.gif]
[Graphics:../Images/SplineQuadMod_gr_55.gif]
[Graphics:../Images/SplineQuadMod_gr_56.gif]
[Graphics:../Images/SplineQuadMod_gr_57.gif]

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

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


[Graphics:../Images/SplineQuadMod_gr_59.gif]
[Graphics:../Images/SplineQuadMod_gr_60.gif]


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

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

[Graphics:../Images/SplineQuadMod_gr_63.gif]
[Graphics:../Images/SplineQuadMod_gr_64.gif]
[Graphics:../Images/SplineQuadMod_gr_65.gif]
[Graphics:../Images/SplineQuadMod_gr_66.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004