Absolute and Relative Error Control in Composite Interpolatory Quadrature: the CIRQUE Algorithm

Justin Steven Prentice


We introduce the CIRQUE algorithm, for approximating definite integrals of continuous, univariate, real-valued functions, using positive-coefficient composite interpolatory quadrature. CIRQUE estimates and controls absolute and/or relative error, without the need for a prior estimate of the magnitude of the integral. The limiting effects of roundoff error are catered for, and CIRQUE is able to provide estimates of error bounds as output. Moreover, if these bounds are deemed too large, it is a simple matter to rerun CIRQUE once to obtain an acceptable bound. We have demonstrated the algorithm using the Trapezium rule, Simpson's rule and four-point Gauss-Legendre quadrature.

Full Text:


DOI: http://dx.doi.org/10.5539/jmr.v3n3p63

Journal of Mathematics Research   ISSN 1916-9795 (Print)   ISSN 1916-9809 (Online)

Copyright © Canadian Center of Science and Education

To make sure that you can receive messages from us, please add the 'ccsenet.org' domain to your e-mail 'safe list'. If you do not receive e-mail in your 'inbox', check your 'bulk mail' or 'junk mail' folders.