MITTAG_LEFFLER Generalized Mittag-Leffler function.
E = mittag_leffler(alpha, beta, z) evaluates the Mittag-Leffler function for each element of z with parameters alpha and beta, where alpha, beta are real scalars and alpha > 0.

E = mittag_leffler(alpha, beta, z, rho) evaluates the Mittag-Leffler function with accuracy rho. Other parameters are same as before. The default value of rho is 10^(-12).

The Generalized Mittag-Leffler function is defined as E_{alpha, beta}(z) = sum k from 0 to inf of Z^k/gamma(alpha*k + beta)

Reference:
(1) Rudolf Gorenflo, Yuri Luchko, and Joulia Loutchko Computation of the Mittag-Leffler function and its derivatives. Fractional Calculus & Applied Analysis. 2002, 5(4):12-15.
(2) MLF (alpha,beta,Z,P) written by Igor Podlubny, Martin Kacenak.