Example 5.  Use Romberg integration to compute a numerical approximation to the integral .  
Use the tolerances .  Compare with the analytic or "true value" of the integral.

Solution 5.

5 (a). Plot the function over the interval  [0, 3].

5 (b). Construct the Romberg table using   tol = 0.001

5 (c). The last entry in the table is  .  Let's find it.

5 (d). Look at 10 digits in  .

5 (e). Are all 10 digits correct ?   Why ?  
NO.  The subroutine was called with the accuracy  .  

5 (f). Determine how to call the Romberg integration subroutine so that it will achieve 10 digits of accuracy.   

We suspect the following answer might have 10 digits of accuracy.

Since the last row was row 7, the sequential trapezoidal rule used the following number of function calls.

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

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

(c) John H. Mathews 2004