Numerical quadrature by the -algorithm

Author:
David K. Kahaner

Journal:
Math. Comp. **26** (1972), 689-693

MSC:
Primary 65D30

DOI:
https://doi.org/10.1090/S0025-5718-1972-0329210-X

MathSciNet review:
0329210

Full-text PDF Free Access

Abstract | References | Similar Articles | Additional Information

Abstract: 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 -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]**Philip J. Davis and Philip Rabinowitz,*Numerical integration*, Blaisdell Publishing Co. Ginn and Co., Waltham, Mass.-Toronto, Ont.-London, 1967. MR**0211604****[2]**F. L. Bauer, H. Rutishauser, and E. Stiefel,*New aspects in numerical quadrature*, Proc. Sympos. Appl. Math., Vol. XV, Amer. Math. Soc., Providence, R.I., 1963, pp. 199–218. MR**0174177****[3]**Israel Navot,*An extension of the Euler-Maclaurin summation formula to functions with a branch singularity*, J. Math. and Phys.**40**(1961), 271–276. MR**0140876****[4]**Israel Navot,*The Euler-Maclaurin functional for functions with a quasi-step discontinuity*, Math. Comp.**17**(1963), 337–345. MR**0155429**, https://doi.org/10.1090/S0025-5718-1963-0155429-1**[5]**I. Navot, ``A further extension of the Euler-Maclaurin summation formula,''*J. Mathematical Phys.*, v. 39, 1962, pp. 155-163.**[6]**J. N. Lyness and B. W. Ninham,*Numerical quadrature and asymptotic expansions*, Math. Comp.**21**(1967), 162–178. MR**0225488**, https://doi.org/10.1090/S0025-5718-1967-0225488-X**[7]**B. W. Ninham and J. N. Lyness,*Further asymptotic expansions for the error functional*, Math. Comp.**23**(1969), 71–83. MR**0242351**, https://doi.org/10.1090/S0025-5718-1969-0242351-0**[8]**L. Fox,*Romberg integration for a class of singular integrands*, Comput. J.**10**(1967), 87–93. MR**0213013**, https://doi.org/10.1093/comjnl/10.1.87**[9]**L. Fox and Linda Hayes,*On the definite integration of singular integrands*, SIAM Rev.**12**(1970), 449–457. MR**0270557**, https://doi.org/10.1137/1012084**[10]**T. Havie, ``On the practical application of the modified Romberg algorithm,''*Nordisk Tidskr. Informationsbehandling*, v. 7, 1967, pp. 103-113.**[11]**C. deBoor,*On Writing an Automatic Integration Algorithm*, Proc. Math. Software Sympos., Purdue University, West Lafayette, Ind., 1970.**[12]**D. Shanks,*An Analogy Between Transients and Mathematical Sequence and Some Nonlinear Sequence to Sequence Transforms Suggested by It*, Naval Ordnance Laboratory Memorandum 9994, July 1949.**[13]**P. Wynn,*On a device for computing the 𝑒_{𝑚}(𝑆_{𝑛}) tranformation*, Math. Tables Aids Comput.**10**(1956), 91–96. MR**0084056**, https://doi.org/10.1090/S0025-5718-1956-0084056-6**[14]**P. Wynn,*On the convergence and stability of the epsilon algorithm*, SIAM J. Numer. Anal.**3**(1966), no. 1, 91–122. MR**0207180**, https://doi.org/10.1137/0703007**[15]**P. Wynn,*The numerical transformation of slowly convergent series by methods of comparison. I*, Chiffres**4**(1961), 177–210. MR**0162350****[16]**D. Kahaner,*Numerical Quadrature by the -Algorithm*, Los Alamos Scientific Laboratory Report LA-DC-12844, September 1971.**[17]**P. Wynn,*Upon a Conjecture Concerning a Method for Solving Linear Equations, and Certain Other Matters*, MRC Technical Summary Report No. 626, 1966, p. 10.

Retrieve articles in *Mathematics of Computation*
with MSC:
65D30

Retrieve articles in all journals with MSC: 65D30

Additional Information

DOI:
https://doi.org/10.1090/S0025-5718-1972-0329210-X

Keywords:
-algorithm,
acceleration of convergence,
extrapolation,
numerical quadrature

Article copyright:
© Copyright 1972
American Mathematical Society