This is a MATLAB routine for evaluating the MittagLeffler function with two parameters (sometimes also called generalized exponential function).
The MittagLeffler function with two parameters plays an important role and appears frequently in solutions of fractional differential equations (i.e. differential equations containing fractional derivatives).
USAGE:
MLF(alpha,beta,Z,P) is the MittagLeffler function E_{alpha,beta}(Z) evaluated with accuracy 10^(P) for each element of Z.
Update 20090325:
(1) Now Z can also be a twodimensional array.
(2) Addressed the issue reported by Li Jackie. It was caused by rounding errors during computations.
Update 20120907:
(1) Addressed the potential issue reported by Brian Bigler.
(2) Code clean up.
