Affiliation:
1. Department of Mathematics, King Fahd University of Petroleum & Minerals, Dhahran 31261, Saudi Arabia
2. Department of Mathematical Sciences, Middle Tennessee State University, Murfreesboro, TN 37132, USA
Abstract
The two-parameter Mittag–Leffler function Eα,β is of fundamental importance in fractional calculus, and it appears frequently in the solutions of fractional differential and integral equations. However, the expense of calculating this function often prompts efforts to devise accurate approximations that are more cost-effective. When α>1, the monotonicity property is largely lost, resulting in the emergence of roots and oscillations. As a result, current rational approximants constructed mainly for α∈(0,1) often fail to capture this oscillatory behavior. In this paper, we develop computationally efficient rational approximants for Eα,β(−t), t≥0, with α∈(1,2). This process involves decomposing the Mittag–Leffler function with real roots into a weighted root-free Mittag–Leffler function and a polynomial. This provides approximants valid over extended intervals. These approximants are then extended to the matrix Mittag–Leffler function, and different implementation strategies are discussed, including using partial fraction decomposition. Numerical experiments are conducted to illustrate the performance of the proposed approximants.
Reference32 articles.
1. Kilbas, A., Srivastava, H., and Trujullo, J. (2006). Theory and Applications of Fractional Differential Equations, Elsevier.
2. Gorenflo, R., Kilbas, A.A., Mainardi, F., and Rogosin, S.V. (2014). Mittag-Leffler Functions. Related Topics and Applications, Springer.
3. Gorenflo, R., and Mainardi, F. (1996). Fractional Oscillations and Mittag-Leffler Functions, Citeseer.
4. Fractional relaxation-oscillation and fractional diffusion-wave phenomena;Mainardi;Chaos Solitons Fractals,1996
5. Dynamics of the fractional oscillator;Achar;Phys. A Stat. Mech. Its Appl.,2001