Romberg Rule of Integration
Romberg Rule of Integration
Romberg Rule of Integration
What is integration?
Integration is the process of measuring the area under a function plotted on a graph. Why
would we want to integrate a function? Among the most common examples are finding the
velocity of a body from an acceleration function, and displacement of a body from a velocity
function. Throughout many engineering fields, there are (what sometimes seems like)
countless applications for integral calculus.
Sometimes, the evaluation of expressions involving these integrals can become daunting, if
not indeterminate. For this reason, a wide variety of numerical methods has been developed
to simplify the integral.
Here, we will discuss the Romberg rule of approximating integrals of the form
b
I f x dx (1)
a
where
f (x) is called the integrand
a lower limit of integration
b upper limit of integration
f x dx
a
07.04.1
07.04.2 Chapter 07.04
is given by
n
f i
b a
3
Et i 1
(2)
12n 2 n
where for each i , i is a point somewhere in the domain a i 1h, a ih , and
n
f i
the term i 1
can be viewed as an approximate average value of f x in a, b . This
n
leads us to say that the true error E t in Equation (2) is approximately proportional to
1
Et 2 (3)
n
b
for the estimate of f x dx using the n -segment trapezoidal rule.
a
Table 1 shows the results obtained for
30
140000
8 2000 ln 140000 2100t 9.8t dt
using the multiple-segment trapezoidal rule.
Approximate
n Et t % a %
Value
1 11868 807 7.296 ---
2 11266 205 1.854 5.343
3 11153 91.4 0.8265 1.019
4 11113 51.5 0.4655 0.3594
5 11094 33.0 0.2981 0.1669
6 11084 22.9 0.2070 0.09082
7 11078 16.8 0.1521 0.05482
8 11074 12.9 0.1165 0.03560
The true error for the 1-segment trapezoidal rule is 807 , while for the 2-segment rule, the
true error is 205 . The true error of 205 is approximately a quarter of 807 . The true
error gets approximately quartered as the number of segments is doubled from 1 to 2. The
same trend is observed when the number of segments is doubled from 2 to 4 (the true error
for 2-segments is 205 and for four segments is 51.5 ). This follows Equation (3).
This information, although interesting, can also be used to get a better approximation of the
integral. That is the basis of Richardson’s extrapolation formula for integration by the
trapezoidal rule.
Romberg rule of Integration 07.04.3
Example 1
The vertical distance in meters covered by a rocket from t 8 to t 30 seconds is given by
30
140000
x 2000 ln 9.8t dt
8 140000 2100t
a) Use Romberg’s rule to find the distance covered. Use the 2-segment and 4-segment
trapezoidal rule results given in Table 1.
b) Find the true error for part (a).
c) Find the absolute relative true error for part (a).
Solution
a) I 2 11266 m
I 4 11113 m
Using Richardson’s extrapolation formula for the trapezoidal rule, the true value is given by
I In
TV I 2 n 2 n
3
and choosing n 2 ,
I I2
TV I 4 4
3
07.04.4 Chapter 07.04
11113 11266
11113
3
11062 m
b) The exact value of the above integral is
30
140000
x 2000 ln 9.8t dt
8 140000 2100t
11061m
so the true error
Et True Value Approximate Value
11061 11062
1 m
c) The absolute relative true error, t , would then be
True Error
t 100
True Value
11061 11062
100
11061
0.00904%
Table 2 shows the Richardson’s extrapolation results using 1, 2, 4, and 8 segments. Results
are compared with those of the trapezoidal rule.
Table 2 Values obtained using Richardson’s extrapolation formula for the trapezoidal rule
for
30
140000
x 2000 ln 9.8t dt .
8 140000 2100t
Romberg Integration
Romberg integration is the same as Richardson’s extrapolation formula as given by
Equation (8) . However, Romberg used a recursive algorithm for the extrapolation as
follows.
The estimate of the true error in the trapezoidal rule is given by
Romberg rule of Integration 07.04.5
3 f i
Et
b a i 1
12n 2 n
Since the segment width, h , is given by
ba
h
n
Equation (2) can be written as
n
f
Et
h2 b a i 1
i
(9)
12 n
The estimate of true error is given by
Et Ch 2 (10)
It can be shown that the exact true error could be written as
Et A1h 2 A2 h 4 A3 h 6 ... (11)
and for small h ,
Et A1h 2 O h 4 (12)
Since we used Et Ch in the formula (Equation (12)), the result obtained from
2
Equation (10) has an error of O h 4 and can be written as
I 2 n R I 2 n I 2 n I n
3
I I
I 2 n 22n1 n (13)
4 1
where the variable TV is replaced by I 2 n R as the value obtained using Richardson’s
extrapolation formula. Note also that the sign is replaced by the sign =.
Hence the estimate of the true value now is
TV I 2n R Ch 4
Determine another integral value with further halving the step size (doubling the number of
segments),
I 4 n R I 4 n I 4 n I 2 n (14)
3
then
4
h
TV I 4 n R C
2
From Equation (13) and (14),
I 4 n R I 2 n R
TV I 4 n R
15
I 4n R I 2n R
I 4 n R (15)
431 1
07.04.6 Chapter 07.04
The above equation now has the error of O h6 . The above procedure can be further
improved by using the new values of the estimate of the true value that has the error of O h6
to give an estimate of O h8 .
Based on this procedure, a general expression for Romberg integration can be written as
I I
I k , j I k 1, j 1 k 1, j k11 k 1, j , k 2 (16)
4 1
The index k represents the order of extrapolation. For example, k 1 represents the values
obtained from the regular trapezoidal rule, k 2 represents the values obtained using the
true error estimate as O h 2 , etc. The index j represents the more and less accurate estimate
of the integral. The value of an integral with a j 1 index is more accurate than the value of
the integral with a j index.
For k 2 , j 1 ,
I 1, 2 I 1,1
I 2,1 I 1, 2
4 21 1
I1, 2 I1,1
I1, 2
3
For k 3 , j 1 ,
I 2, 2 I 2,1
I 3,1 I 2, 2
4 31 1
I 2, 2 I 2,1
I 2, 2 (17)
15
Example 2
The vertical distance in meters covered by a rocket from t 8 to t 30 seconds is given by
30
140000
x 2000 ln 9.8t dt
8 140000 2100t
Use Romberg’s rule to find the distance covered. Use the 1, 2, 4, and 8-segment trapezoidal
rule results as given in Table 1.
Solution
From Table 1, the needed values from the original the trapezoidal rule are
I1,1 11868
I1, 2 11266
I1,3 11113
I1, 4 11074
where the above four values correspond to using 1, 2, 4 and 8 segment trapezoidal rule,
respectively. To get the first order extrapolation values,
Romberg rule of Integration 07.04.7
I1, 2 I1,1
I 2,1 I1, 2
3
11266 11868
11266
3
11065
Similarly
I 1,3 I 1, 2
I 2, 2 I 1,3
3
11113 11266
11113
3
11062
I I
I 2,3 I1, 4 1, 4 1,3
3
11074 11113
11074
3
11061
For the second order extrapolation values,
I I 2,1
I 3,1 I 2, 2 2, 2
15
11062 11065
11062
15
11062
Similarly
I I
I 3, 2 I 2,3 2,3 2, 2
15
11061 11062
11061
15
11061
For the third order extrapolation values,
I I
I 4,1 I 3, 2 3, 2 3,1
63
11061 11062
11061
63
11061m
Table 3 shows these increasingly correct values in a tree graph.