Numerical Evaluation of Two and Three Parameter Mittag-Leffler Functions
Abstract
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. . Al-Bassam, Some existence theorems on differential equations of generalized order , J. Reine Angew. Math. , 218 ( 1965 ), pp. 70 -- 78 .
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.
3. , Jr., Models based on Mittag--Leffler functions for anomalous relaxation in dielectrics , Eur. Phys. J. Special Topics , 193 ( 2011 ), pp. 161 -- 171 .
4. , Linear models of dissipation in anelastic solids, Riv. Nuovo Cimento (Ser . II) , 1 ( 1971 ), pp. 161 -- 198 .
5. , 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.
6. P. J. Davis and P. Rabinowitz, Methods of Numerical Integration, 2nd ed., Comput. Sci. Appl. Math., Academic Press Inc., Orlando, FL, 1984.
7. . Weideman, An improved Talbot method for numerical Laplace transform inversion , Numer. Algorithms , 68 ( 2015 ), pp. 167 -- 183 .
8. , Fractional derivatives and the Cauchy problem for differential equations of fractional order , Izv. Akad. Nauk Armjan. SSR Ser. Mat. , 3 ( 1968 ), pp. 3 -- 29 .
9. , The asymptotic expansion of integral functions defined by generalized hypergeometric functionss , Proc. London Math. Soc. , 27 ( 1928 ), pp. 389 -- 400 .
10. R. Garrappa, The Mittag--Leffler function, MATLAB Central File Exchange, 2014, file ID: 48154.
11. , Solving the time-fractional Schrödinger equation by Krylov projection methods , J. Comput. Phys. , 293 ( 2015 ), pp. 115 -- 134 .
12. , Generalized exponential time differencing methods for fractional order problems , Comput. Math. Appl. , 62 ( 2011 ), pp. 876 -- 890 .
13. , Evaluation of generalized Mittag--Leffler functions on the real line , Adv. Comput. Math. , 39 ( 2013 ), pp. 205 -- 225 .
14. , Exponentially convergent algorithms for the operator exponential with applications to inhomogeneous problems in Banach spaces , SIAM J. Numer. Anal. , 43 ( 2005 ), pp. 2144 -- 2171 .
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.
16. , Computation of the Mittag-Leffler function $E_{\alpha,\beta}(z)$ and its derivative , Fract. Calc. Appl. Anal. , 5 ( 2002 ), pp. 491 -- 518 .
17. , Mittag-Leffler functions and their applications , J. Appl. Math., ( 2011 ), 298628 .
18. , Computation of the generalized Mittag-Leffler function and its inverse in the complex plane , Integral Transforms Spec. Funct. , 17 ( 2006 ), pp. 637 -- 652 .
19. , On the theory of linear integral equations , Ann. of Math. (2) , 31 ( 1930 ), pp. 479 -- 528 .
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.
21. , A spectral order method for inverting sectorial Laplace transforms , SIAM J. Numer. Anal. , 44 ( 2006 ), pp. 1332 -- 1350 .
22. F. Mainardi, Fractional Calculus and Waves in Linear Viscoelasticity, Imperial College Press, London, 2010.
23. , On Mittag-Leffler-type functions in fractional evolution processes , J. Comput. Appl. Math. , 118 ( 2000 ), pp. 283 -- 299 .
24. . Paris (Ser . II) , 136 ( 1902 ), pp. 937 -- 939 .
25. , Sopra la funzione ${E}_{\alpha}(x)$ , Rend. Accad. Lincei , 13 ( 1904 ), pp. 3 -- 5 .
26. , A note on Krylov methods for fractional evolution problems , Numer. Funct. Anal. Optim. , 34 ( 2013 ), pp. 539 -- 556 .
27. , On the convergence of Krylov subspace methods for matrix Mittag--Leffler functions , SIAM J. Numer. Anal. , 49 ( 2011 ), pp. 2144 -- 2164 .
28. , Algorithm 682: Talbot's method of the Laplace inversion problems , ACM Trans. Math. Software , 16 ( 1990 ), pp. 158 -- 168 .
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.
30. I. Podlubny and M. Kacenak, The MATLAB MLF Code, MATLAB Central File Exchange, 2001--2012, file ID: 8738.
31. , A singular integral equation with a generalized Mittag-Leffler function in the kernel , Yokohama Math. J. , 19 ( 1971 ), pp. 7 -- 15 .
32. , Numerical algorithm for calculating the generalized Mittag-Leffler function , SIAM J. Numer. Anal. , 47 ( 2008 ), pp. 69 -- 88 .
33. , A parallel method for time discretization of parabolic equations based on Laplace transformation and quadrature , IMA J. Numer. Anal. , 23 ( 2003 ), pp. 269 -- 299 .
34. , The accurate numerical inversion of Laplace transforms , J. Inst. Math. Appl. , 23 ( 1979 ), pp. 97 -- 120 .
35. . Weideman, The exponentially convergent trapezoidal rule , SIAM Rev. , 56 ( 2014 ), pp. 385 -- 458 .
36. , Talbot quadratures and rational approximations , BIT , 46 ( 2006 ), pp. 653 -- 670 .
37. F. G. Tricomi, Funzioni ipergeometriche confluenti, Edizione Cremonese, Roma, 1954.
38. contours for the inversion of the Laplace transform , SIAM J. Numer. Anal. , 44 ( 2006 ), pp. 2342 -- 2362 .
39. . Weideman, Improved contour integral methods for parabolic PDEs , IMA J. Numer. Anal. , 30 ( 2010 ), pp. 334 -- 350 .
40. , Parabolic and hyperbolic contours for computing the Bromwich integral , Math. Comp. , 76 ( 2007 ), pp. 1341 -- 1356 .
41. , Über den fundamentalsatz in der teorie der funktionen ${E}_{\alpha}(x)$ , Acta Math. , 29 ( 1905 ), pp. 191 -- 201 .
42.
Wolfram Research Inc. , Tricomi Confluent Hypergeometric Function , 1998 -- 2010 , http://functions.wolfram.com/HypergeometricFunctions/HypergeometricU/.43. E. M. Wright, The asymptotic expansion of the generalized hypergeometric function, J. London Math. Soc., s-1-10 (1935), pp. 286--293.


