Computing Narrow Inclusions for Definite Integrals.
Abstract
The ACRITH computer program is used to implement algorithms which give narrow inclusions for definite integrals. Inclusions are obtained from familiar numerical quadrature formulas such as Gaussian or Newton-Cotes with remainder terms, or from the term-by-term integration of Taylor polynomials with remainder terms. Inclusions for the remainder terms are computed using automatic differentiation. The inclusions for the remainder terms are computed using automatic differentiation. The inclusions are valid if the integrand or the endpoints of the interval of integration are real-or interval-valued. Interval inclusions which contain only a few machine numbers are achieved by using ACRITH's accurate scalar product, by using order and subinterval adaptation, and by using special devices such as intersection of several estimates. Numerical examples show that such narrow, guaranteed bounds require about four times as long to compute as the estimates computed by the routine QAGS from QUADPACK. (Author)
Document Details
- Document Type
- Technical Report
- Publication Date
- Feb 01, 1986
- Accession Number
- ADA167491
Entities
People
- George F. Corliss
Organizations
- University of Wisconsin–Madison