| Home > Publications > Reports > Numerical Analysis and Applied Mathematics (TW) |
TW 416
Joris Van Deun, Ronald Cools
Computing infinite range integrals of an arbitrary product of Bessel functions
Abstract
We present an algorithm to compute integrals of the form
∫0∞ xm ∏i=1k Jνi(ai x) dx
with Jνi(x) the Bessel function of the first kind and order νi. The parameter m is a real number such that ∑i νi +m > -1 and the coefficients ai are strictly positive real numbers. The main ingredients in this algorithm are the well-known asymptotic expansion for Jνi(x) and the observation that the infinite part of the integral can be approximated using the incomplete Gamma function Γ(a,z). Accurate error estimates are included in the algorithm, which is implemented as a Matlab program. report.pdf (201K) / mailto: J. Van Deun