Modifications of Romberg's integration method are applicable to functions with endpoint singularities. Several authors have reached different conclusions on the usefulness of these schemes due to the potential difficulties in calculating certain exponents in the asymptotic error expansion. In this paper, we consider a method based on the e-algorithm which does not require the user to supply these parameters. Tests show this luxury produces a method substantially better than the unmodified Romberg's method, but not as good as the modified procedure which assumes all the exponents are known. It is possible to design a scheme which incorporates the best of both methods and allows the user to decide how much information he wishes to provide. The algorithm is stable numerically in quadrature applications and converges for functions with endpoint singularities of an algebraic or logarithmic nature.
[1]
Israel Navot,et al.
A Further Extension of the Euler-Maclaurin Summation Formula
,
1962
.
[2]
L. Fox,et al.
On the Definite Integration of Singular Integrands
,
1970
.
[3]
I. Navot.
The Euler-Maclaurin functional for functions with a quasi-step discontinuity
,
1963
.
[4]
P. Wynn,et al.
On a device for computing the _{}(_{}) tranformation
,
1956
.
[5]
P. Wynn,et al.
On the Convergence and Stability of the Epsilon Algorithm
,
1966
.
[6]
L. Fox,et al.
Romberg integration for a class of singular integrands
,
1967,
Comput. J..
[7]
J. N. Lyness,et al.
Numerical quadrature and asymptotic expansions
,
1967
.
[8]
Israel Navot,et al.
An Extension of the Euler‐Maclaurin Summation Formula to Functions with a Branch Singularity
,
1961
.