计算Mittag-Leffler函数的MATLAB程序

function E = ml（zalphabetagama）
%
% Evaluation of the Mittag-Leffler （ML） function with 1 2 or 3 parameters
% by means of the OPC algorithm [1]. The routine evaluates an approximation
% Et of the ML function E such that |E-Et|/（1+|E|） approx 1.0e-15   
%     
%
% E = ML（zalpha） evaluates the ML function with one parameter alpha for
% the corresponding elements of z; alpha must be a real and positive
% scalar. The one parameter ML function is defined as
%
% E = sum_{k=0}^{infty} z^k/Gamma（alpha*k+1）
%
% with Gamma the Euler‘s gamma function.
%
%
% E = ML（zalphabeta） evaluates the ML function with two parameters alpha
% and beta for the corresponding elements of z; alpha must be a real and
% positive scalar and beta a real scalar. The two parameters ML function is
%