Computation of infinite integrals involving Bessel functions of arbitrary order by the (D)over-bar-transformation

[摘要] The D-transformation due to the author is an effective extrapolation method for computing infinite oscillatory integrals of various kinds. In this work two new variants of this transformation are designed for computing integrals of the form integral(a)(infinity)y(t)C-v(t) dt, where g(x) is a nonoscillatory function and C-v(x) may be an arbitrary linear combination of the Bessel functions of the first and second kinds J(v)(x) and Y-v(x), of arbitrary real order v. When applied to such integrals, the D-transformation and its new variants are observed to produce very accurate results. It is also seen that their performance is very similar to that of the modified W-transformation due to the author, as extended in a recent work by Lucas and Stone with C-v(x) = J(v)(x). The present paper is concluded by stating the relevant convergence and stability results and by appending a numerical example.

[发布日期] 1997-02-03

[关键词] numerical integration;infinite oscillatory integrals;generalized Richardson extrapolation;Bessel functions;Hankel transforms