Numerical Evaluation of Mittag-Leffler Functions

The Mittag-Leffler function is computed via a quadrature approximation of a contour integral representation. We compare results for parabolic and hyperbolic contours, and give special attention to evaluation on the real line. The main point of difference with respect to similar approaches from the literature is the way that poles in the integrand are handled. Rational approximation of the Mittag-Leffler function on the negative real axis is also discussed.