WIT Press


A Mathematica Implementation Of Double Adaptive Quadrature

Price

Free (open access)

Paper DOI

10.2495/IMS970201

Volume

15

Pages

7

Published

1997

Size

411 kb

Author(s)

P. Favati, G. Florentino, G. Lotti & F. Romani

Abstract

This paper presents a Mathematica implementation of a recent algorithm for efficient high precision automatic quadrature. Extensive numerical testing has been performed on a comprehensive set of functions showing that our algorithm is more efficient than the built-in Mathematica Gauss-Kronrod quadrature algorithm. 1 Introduction We consider the problem of computing an approximation of a definite integral by using an automatic integration routine. In a couple of recent papers [1,2] a new algorithm was introduced (Double Adaptive Quadrature, hereafter DAQ) which is suitable for high precision adaptive automatic quadrature. The basis of DAQ, in a general adaptive scheme, is the combination of the two main strategies to improve the approximation of an integral, namely interval subdivision and application of more accurate formulas. The active subintervals are stored in a q

Keywords