On the Numerical Computation of the Mittag–Leffler Function

Abstract

Abstract The Mittag–Leffler function (MLF) plays an important role in many applications of fractional calculus, establishing a connection between exponential and power law behaviors that characterize integer and fractional order phenomena, respectively. Nevertheless, the numerical computation of the MLF poses problems both of accuracy and convergence. In this paper, we study the calculation of the 2-parameter MLF by using polynomial computation and integral formulas. For the particular cases having Laplace transform (LT) the method relies on the inversion of the LT using the fast Fourier transform. Experiments with two other available methods compare also the computational time and accuracy. The 3-parameter MLF and its calculation are also considered.