Global-Padé Approximation of the Three-Parameter Mittag-Leffler Function: Generalized Derivation and Numerical Implementation Issues

Abstract

In this paper, we discuss a framework for efficient construction of rational approximations for the three-parameter Mittag-Leffler function (MLF) based on the global-Padé approximation technique. Obtaining the coefficients of these approximants is usually the most tedious part in their construction as it requires delicate matching to procure a linear system for these coefficients. The main focus of this paper is the derivation of a novel generalized system that can be used to obtain the coefficients without having to perform the matching task. In particular, we illustrate the use of the generalized system in constructing rational approximants of various degrees subject to choices of feasible parameters. Inequalities providing bounds on permissible choices of these parameters are also obtained. Numerical experiments are conducted to illustrate the accuracy and efficiency of the approximants.

Keywords

Prabhakar function / Mittag-Leffler function (MLF) / Global rational approximation / Multi-term fractional differential equation / Fractional heat conduction equation / 33F05 / 33B10 / 41A20 / 41A21 / 41A80

References

[1] Abdelhedi W. Fractional differential equations with a ψ\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\psi $$\end{document}-Hilfer fractional derivative. Comput. Appl. Math.. 2021, 40: 53.

[2] Ali I, Malik NA. Hilfer fractional advection-diffusion equations with power-law initial condition; a numerical study using variational iteration method. Comput. Math. Appl.. 2014, 68(10): 1161-1179.

[3] Atkinson C, Osseiran A. Rational solutions for the time-fractional diffusion equation. SIAM J. Appl. Math.. 2011, 71192-106.

[4] Borhanifar A, Valizadeh S. Mittag-Leffler-Padé approximations for the numerical solution of space and time fractional diffusion equations. Int. J. Appl. Math. Res.. 2015, 4(4): 466.

[5] Furati KM, Mustapha K, Sarumi IO, Iyiola OS. Inverse source in two-parameter anomalous diffusion, numerical algorithms, and simulations over graded time meshes. Comput. Appl. Math.. 2021, 40125.

[6] Furati KM, Sarumi IO, Khaliq AQM. Fractional model for the spread of Covid-19 subject to government intervention and public perception. Appl. Math. Modell.. 2021, 95: 89-105.

[7] Garra, R., Garrappa, R.: The Prabhakar or three parameter Mittag-Leffler function: theory and application. Commun. Nonlinear Sci. Num. Simul. 56, 314–329 (2018)

[8] Garra, R., Gorenflo, R., Polito, F., Tomovski, Z.: Hilfer-Prabhakar derivatives and some applications. Appl. Math. Comput. 242, 576–589 (2014)

[9] Garrappa, R.: Numerical evaluation of two and three parameter Mittag-Leffler functions. SIAM J. Numer. Anal. 53(3), 1350–1369 (2015)

[10] Garrappa R, Mainnardi F. Models of dielectric relaxation based on completely monotone functions. Fract. Calc. Appl. Anal.. 2016, 1951105-1160.

[11] Garrappa R, Popolizio M. Generalized exponential time differencing methods for fractional order problems. Comput. Math. Appl.. 2011, 62: 876-890.

[12] Gorenflo R, Kilbas AA, Mainardi F, Rogosin SV. Mittag-Leffler Functions Related Topics and Applications. 2014, Berlin, Springer.

[13] Gorenflo R, Loutchko J, Luchko Y. Computation of the Mittag-Leffler function Eα,β(z)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$E_{\alpha,\beta }(z)$$\end{document} and its derivative. Fract. Calc. Appl. Anal.. 2002, 5(4): 491-518

[14] Ingo C, Barrick TR, Webb AG, Ronen I. Accurate Padé global approximations for the Mittag-Leffler function, its inverse, and its partial derivatives to efficiently compute convergent power series. Int. J. Appl. Comput. Math.. 2017, 3(2): 347-362.

[15] Iyiola OS, Asante-Asamani EO, Wade BA. A real distinct poles rational approximation of generalized Mittag-Leffler functions and their inverses: applications to fractional calculus. J. Comput. Appl. Math.. 2018, 330: 307-317.

[16] Kilbas AA, Srivastava HM, Trujullo JJ. Theory and Applications of Fractional Differential Equations. 2006, Amsterdam, Elsevier

[17] Luchko Y, Gorenflo R. An operational method for solving fractional differential equations with the Caputo derivatives. Acta Math. Vietnamica. 1999, 24: 207-233

[18] Pandey SC. The Lorenzo-Hartley’s function for fractional calculus and its applications pertaining to fractional order modelling of anomalous relaxation in dielectrics. Comput. Appl. Math.. 2018, 37: 2648-2666.

[19] Sarumi, I.O.: Reaction coefficient identification problem for a time-fractional diffusion equation. In: 2023 International Conference on Fractional Differentiation and Its Applications (ICFDA), Ajman, United Arab Emirates, pp. 1–6 (2023)

[20] Sarumi IO, Furati KM, Khaliq AQM. Highly accurate global Padé approximations of generalized Mittag-Leffler function and its inverse. J. Sci. Comput.. 2020, 82: 46.

[21] Sarumi, I.O., Furati, K.M., Khaliq, A.Q.M., Mustapha, K.: Generalized exponential time differencing schemes for stiff fractional systems with nonsmooth source term. J. Sci. Comput. 86, 23 (2021). https://doi.org/10.1007/s10915-020-01374-y

[22] Stanislavsky A, Weron K. A typical case of the dielectric relaxation responses and its fractional kinetic equation. Fract. Calc. Appl. Anal.. 2016, 19(1): 212-228.

[23] Starovoïtov AP, Starovoïtova NA. Padé approximants of the Mittag-Leffler functions. Sbornik Math.. 2007, 19871011-1023.

[24] Vieira N, Rodrigues MM, Ferreira M. Time-fractional diffusion equation with ψ\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\psi $$\end{document}-Hilfer derivative. Comput. Appl. Math.. 2022, 41(6): 230.

[25] Winitzki, S.: Uniform approximations for transcendental functions. In: Kumar, V., Gavrilova, M.L., Jeng Kenneth Tan, C., L’Ecuyer, P. (eds) Computational Science and Its Applications — ICCSA 2003, pp. 780–789. Springer, Berlin, Heidelberg (2003)

[26] Zakariya YF, Afolabi YO, Nuruddeen RI, Sarumi IO. Analytical solutions to fractional fluid flow and oscillatory process models. Fractal Fract.. 2018, 2181-12

[27] Zeng CB, Chen YQ. Global Padé approximations of the generalized Mittag-Leffler function and its inverse. Fract. Calc. Appl. Anal.. 2015, 18(6): 1492-1506.

[28] Zeng, F., Zhang, Z., Karniadakis, G.E.: Second-order numerical methods for multi-term fractional differential equations: smooth and non-smooth solutions. Comput. Methods Appl. Mech. Engrg. 327(1), 478–502 (2017)

RIGHTS & PERMISSIONS

Shanghai University