A recursive algorithm for the G transformation and accurate computation of incomplete Bessel functions
- Autores: Richard M. Slevinsky, Hassan Safouhi
- Localización: Applied numerical mathematics, ISSN-e 0168-9274, Vol. 60, Nº. 12, 2010, págs. 1411-1417
- Idioma: inglés
- DOI: 10.1016/j.apnum.2010.04.005
- Texto completo no disponible (Saber más ...)
-
Resumen
In the present contribution, we develop an efficient algorithm for the recursive computation of the View the MathML source transformation for the approximation of infinite-range integrals. Previous to this contribution, the theoretically powerful View the MathML source transformation was handicapped by the lack of an algorithmic implementation. Our proposed algorithm removes this handicap by introducing a recursive computation of the successive View the MathML source transformations with respect to the order n. This recursion, however, introduces the View the MathML source operator applied to the integrand. Consequently, we employ the Slevinsky�Safouhi formula I for the analytical and numerical developments of these required successive derivatives.
Incomplete Bessel functions, which pose as a numerical challenge, are computed to high pre-determined accuracies using the developed algorithm. The numerical results obtained show the high efficiency of the new method, which does not resort to any numerical integration in the computation.
