Evaluate the Mittag-Leffler (ML) function with two parameters ALPHA and BETA of the square matrix argument A
E = ML(A,ALPHA,BETA) evaluates the ML function with two parameters ALPHA and BETA at the square matrix A argument; ALPHA must be any real and positive scalar, BETA any real scalar and A any real or complex square matrix. If the second parameter BETA is missing, it is assumed BETA=1.
The ML function on the matrix argument A is evaluated by exploiting the Schur-Parlett algorithm and evaluating derivative of the scalar ML function by combining, by means of the derivatives balancing technique studied in , Taylor series, a summation formula based on the Prabhakar function and the numerical inversion of the Laplace transform obtained after generalizing the algorithm described in . For more details users are referred to . Please, cite this code as .
 R. Garrappa and M. Popolizio, Computing the matrix MittagLeffler function with applications to fractional calculus, submitted
 R. Garrappa, Numerical Evaluation of two and three parameter Mittag-Leffler functions, SIAM Journal of Numerical Analysis, 2015, 53(3), 1350-1369.
Inspired by: The Mittag-Leffler function