On a novel iterative method to compute polynomial approximations to Bessel functions of the first kind and its connection to the solution of fractional diffusion/diffusion-wave problems
Abstract
We present an iterative method to obtain approximations to Bessel functions of the first kind Jp(x) (p > -1) via the repeated application of an integral operator to an initial seed function f0(x). The class of seed functions f0(x) leading to sets of increasingly accurate approximations fn(x) is considerably large and includes any polynomial. When the operator is applied once to a polynomial of degree s, it yields a polynomial of degree s + 2, and so the iteration of this operator generates sets of increasingly better polynomial approximations of increasing degree. We focus on the set of polynomial approximations generated from the seed function f0(x) = 1. This set of polynomials is useful not only for the computation of Jp(x) but also from a physical point of view, as it describes the long-time decay modes of certain fractional diffusion and diffusion-wave problems.
- Publication:
-
Journal of Physics A Mathematical General
- Pub Date:
- February 2011
- DOI:
- 10.1088/1751-8113/44/7/075203
- arXiv:
- arXiv:1101.2335
- Bibcode:
- 2011JPhA...44g5203B
- Keywords:
-
- Mathematical Physics;
- Condensed Matter - Statistical Mechanics;
- Mathematics - Classical Analysis and ODEs
- E-Print:
- 14 pages, 4 figures. To be published in J. Phys. A: Math. Theor