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 not only useful 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.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…