Radial quadrature for multiexponential integrands.

We introduce a Gaussian quadrature, based on the polynomials that are orthogonal with respect to the weight function ln(2)x on the interval [0, 1], which is suitable for the evaluation of radial integrals. The quadrature is exact if the non-Jacobian part of the integrand is a linear combination of a geometric sequence of exponential functions. We find that the new scheme is a useful alternative to existing approaches, particularly for integrands that exhibit multiexponential behavior. Copyright 2003 Wiley Periodicals, Inc. J Comput Chem 24: 732-740, 2003