Numerical Evaluation of Two and Three Parameter Mittag-Leffler Functions

The Mittag-Leffler (ML) function plays a fundamental role in fractional calculus but very few methods are available for its numerical evaluation. In this work we present a method for the efficient computation of the ML function based on the numerical inversion of its Laplace transform (LT): an optimal parabolic contour is selected on the basis of the distance and the strength of the singularities of the LT, with the aim of minimizing the computational effort and reducing the propagation of errors. Numerical experiments are presented to show accuracy and efficiency of the proposed approach. The application to the three parameter ML (also known as Prabhakar) function is also presented.

  • 1.  M. . Al-Bassam, Some existence theorems on differential equations of generalized order , J. Reine Angew. Math. , 218 ( 1965 ), pp. 70 -- 78 . ISIGoogle Scholar

  • 2.  J. C. Butcher, On the numerical inversion of Laplace and Mellin transforms, in Proceedings of the Conference on Data Processing and Automatic Computing Machines, Salisbury, Australia, 1957.Google Scholar

  • 3.  E. Capelas de Oliveira F. Mainardi and  J. Vaz , Jr., Models based on Mittag--Leffler functions for anomalous relaxation in dielectrics , Eur. Phys. J. Special Topics , 193 ( 2011 ), pp. 161 -- 171 . CrossrefISIGoogle Scholar

  • 4.  M. Caputo and  F. Mainardi , Linear models of dissipation in anelastic solids, Riv. Nuovo Cimento (Ser . II) , 1 ( 1971 ), pp. 161 -- 198 . Google Scholar

  • 5.  M. Caputo and  F. Mainardi , A new dissipation model based on memory mechanism , Fract. Calc. Appl. Anal. , 10 ( 2007 ), pp. 309 -- 324 ; reprinted from Pure Appl. Geophys., 9 (1971), pp. 134--147. Google Scholar

  • 6.  P. J. Davis and P. Rabinowitz, Methods of Numerical Integration, 2nd ed., Comput. Sci. Appl. Math., Academic Press Inc., Orlando, FL, 1984.Google Scholar

  • 7.  B. Dingfelder and  J. A. . Weideman, An improved Talbot method for numerical Laplace transform inversion , Numer. Algorithms , 68 ( 2015 ), pp. 167 -- 183 . CrossrefISIGoogle Scholar

  • 8.  M. M. Džrbašjan and  A. B. Nersesjan , Fractional derivatives and the Cauchy problem for differential equations of fractional order , Izv. Akad. Nauk Armjan. SSR Ser. Mat. , 3 ( 1968 ), pp. 3 -- 29 . Google Scholar

  • 9.  C. Fox , The asymptotic expansion of integral functions defined by generalized hypergeometric functionss , Proc. London Math. Soc. , 27 ( 1928 ), pp. 389 -- 400 . CrossrefGoogle Scholar

  • 10.  R. Garrappa, The Mittag--Leffler function, MATLAB Central File Exchange, 2014, file ID: 48154.Google Scholar

  • 11.  R. Garrappa I. Moret and  M. Popolizio , Solving the time-fractional Schrödinger equation by Krylov projection methods , J. Comput. Phys. , 293 ( 2015 ), pp. 115 -- 134 . CrossrefISIGoogle Scholar

  • 12.  R. Garrappa and  M. Popolizio , Generalized exponential time differencing methods for fractional order problems , Comput. Math. Appl. , 62 ( 2011 ), pp. 876 -- 890 . CrossrefISIGoogle Scholar

  • 13.  R. Garrappa and  M. Popolizio , Evaluation of generalized Mittag--Leffler functions on the real line , Adv. Comput. Math. , 39 ( 2013 ), pp. 205 -- 225 . CrossrefISIGoogle Scholar

  • 14.  I. P. Gavrilyuk and  V. L. Makarov , Exponentially convergent algorithms for the operator exponential with applications to inhomogeneous problems in Banach spaces , SIAM J. Numer. Anal. , 43 ( 2005 ), pp. 2144 -- 2171 . LinkISIGoogle Scholar

  • 15.  R. Gorenflo, A. A. Kilbas, F. Mainardi, and S. V. Rogosin, Mittag-Leffler functions, Related Topics and Applications, Springer Monogr. Math., Springer, Heidelberg, 2014.Google Scholar

  • 16.  R. Gorenflo J. Loutchko and  Y. Luchko , Computation of the Mittag-Leffler function $E_{\alpha,\beta}(z)$ and its derivative , Fract. Calc. Appl. Anal. , 5 ( 2002 ), pp. 491 -- 518 . Google Scholar

  • 17.  H. J. Haubold A. M. Mathai and  R. K. Saxena , Mittag-Leffler functions and their applications , J. Appl. Math., ( 2011 ), 298628 . Google Scholar

  • 18.  R. Hilfer and  H. J. Seybold , Computation of the generalized Mittag-Leffler function and its inverse in the complex plane , Integral Transforms Spec. Funct. , 17 ( 2006 ), pp. 637 -- 652 . CrossrefISIGoogle Scholar

  • 19.  E. Hille and  J. D. Tamarkin , On the theory of linear integral equations , Ann. of Math. (2) , 31 ( 1930 ), pp. 479 -- 528 . CrossrefGoogle Scholar

  • 20.  A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, North-Holland Math. Stud. 204, Elsevier Science B.V., Amsterdam, 2006.Google Scholar

  • 21.  M. López-Fernández C. Palencia and  A. Schädle , A spectral order method for inverting sectorial Laplace transforms , SIAM J. Numer. Anal. , 44 ( 2006 ), pp. 1332 -- 1350 . LinkISIGoogle Scholar

  • 22.  F. Mainardi, Fractional Calculus and Waves in Linear Viscoelasticity, Imperial College Press, London, 2010.Google Scholar

  • 23.  F. Mainardi and  R. Gorenflo , On Mittag-Leffler-type functions in fractional evolution processes , J. Comput. Appl. Math. , 118 ( 2000 ), pp. 283 -- 299 . CrossrefISIGoogle Scholar

  • 24.  M. G. Mittag-Leffler Sur ́grale de Laplace-Abel and  C. R. Acad . Paris (Ser . II) , 136 ( 1902 ), pp. 937 -- 939 . Google Scholar

  • 25.  M. G. Mittag-Leffler , Sopra la funzione ${E}_{\alpha}(x)$ , Rend. Accad. Lincei , 13 ( 1904 ), pp. 3 -- 5 . Google Scholar

  • 26.  I. Moret , A note on Krylov methods for fractional evolution problems , Numer. Funct. Anal. Optim. , 34 ( 2013 ), pp. 539 -- 556 . CrossrefISIGoogle Scholar

  • 27.  I. Moret and  P. Novati , On the convergence of Krylov subspace methods for matrix Mittag--Leffler functions , SIAM J. Numer. Anal. , 49 ( 2011 ), pp. 2144 -- 2164 . LinkISIGoogle Scholar

  • 28.  A. Murli and  M. Rizzardi , Algorithm 682: Talbot's method of the Laplace inversion problems , ACM Trans. Math. Software , 16 ( 1990 ), pp. 158 -- 168 . CrossrefISIGoogle Scholar

  • 29.  I. Podlubny, Fractional Differential Equations, An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications, Math. Sci. Engrg. 198, Academic Press Inc., San Diego, CA, 1999.Google Scholar

  • 30.  I. Podlubny and M. Kacenak, The MATLAB MLF Code, MATLAB Central File Exchange, 2001--2012, file ID: 8738.Google Scholar

  • 31.  T. R. Prabhakar , A singular integral equation with a generalized Mittag-Leffler function in the kernel , Yokohama Math. J. , 19 ( 1971 ), pp. 7 -- 15 . Google Scholar

  • 32.  H. Seybold and  R. Hilfer , Numerical algorithm for calculating the generalized Mittag-Leffler function , SIAM J. Numer. Anal. , 47 ( 2008 ), pp. 69 -- 88 . LinkISIGoogle Scholar

  • 33.  D. Sheen I. H. Sloan and  V. Thomée , A parallel method for time discretization of parabolic equations based on Laplace transformation and quadrature , IMA J. Numer. Anal. , 23 ( 2003 ), pp. 269 -- 299 . CrossrefISIGoogle Scholar

  • 34.  A. Talbot , The accurate numerical inversion of Laplace transforms , J. Inst. Math. Appl. , 23 ( 1979 ), pp. 97 -- 120 . CrossrefISIGoogle Scholar

  • 35.  L. N. Trefethen and  J. A. . Weideman, The exponentially convergent trapezoidal rule , SIAM Rev. , 56 ( 2014 ), pp. 385 -- 458 . LinkISIGoogle Scholar

  • 36.  L. N. Trefethen J. A. C. Weideman and  T. Schmelzer , Talbot quadratures and rational approximations , BIT , 46 ( 2006 ), pp. 653 -- 670 . CrossrefISIGoogle Scholar

  • 37.  F. G. Tricomi, Funzioni ipergeometriche confluenti, Edizione Cremonese, Roma, 1954.Google Scholar

  • 38.  J. A. C. Weideman and  Optimizing Talbot's contours for the inversion of the Laplace transform , SIAM J. Numer. Anal. , 44 ( 2006 ), pp. 2342 -- 2362 . LinkISIGoogle Scholar

  • 39.  J. A. . Weideman, Improved contour integral methods for parabolic PDEs , IMA J. Numer. Anal. , 30 ( 2010 ), pp. 334 -- 350 . CrossrefISIGoogle Scholar

  • 40.  J. A. C. Weideman and  L. N. Trefethen , Parabolic and hyperbolic contours for computing the Bromwich integral , Math. Comp. , 76 ( 2007 ), pp. 1341 -- 1356 . CrossrefISIGoogle Scholar

  • 41.  A. Wiman , Über den fundamentalsatz in der teorie der funktionen ${E}_{\alpha}(x)$ , Acta Math. , 29 ( 1905 ), pp. 191 -- 201 . CrossrefGoogle Scholar

  • 42.  Wolfram Research Inc. , Tricomi Confluent Hypergeometric Function , 1998 -- 2010 , http://functions.wolfram.com/HypergeometricFunctions/HypergeometricU/. Google Scholar

  • 43.  E. M. Wright, The asymptotic expansion of the generalized hypergeometric function, J. London Math. Soc., s-1-10 (1935), pp. 286--293.Google Scholar