Computation of the Mittag-Leffler function \(E_{\alpha,\beta}(z)\) and its derivative. (English) Zbl 1027.33016

Summary: In this paper algorithms for numerical evaluation of the Mittag-Leffler function \[ E_{\alpha,\beta}(z)= \sum^\infty_{k=0} {z^k\over \Gamma(\beta+\alpha k)},\quad \alpha> 0,\quad \beta\in\mathbb{R},\quad z\in\mathbb{C} \] and its derivative for all values of the parameters \(\alpha> 0\), \(\beta\in\mathbb{R}\) and all values of the argument \(z\in\mathbb{C}\) C are presented. For different parts of the complex plane different numerical techniques are used. In every case we provide estimates for accuracy of the computation; numerous pictures showing the behaviour of the Mittag-Leffler function for different values of the parameters and on different lines in the complex plane are included. The ideas and techniques employed in the paper can be used for numerical evaluation of other functions of the hypergeometric type. In particular, the same method with some small modifications can be applied for the Wright function which plays a very important role in the theory of partial differential equations of fractional order.


33E12 Mittag-Leffler functions and generalizations
30E20 Integration, integrals of Cauchy type, integral representations of analytic functions in the complex plane
33F05 Numerical approximation and evaluation of special functions
65D20 Computation of special functions and constants, construction of tables