Mittag-Leffer函数, Matlab内部函数

Mittag-Leffer函数: $E_{alpha,eta}(x) = sumlimits^{infty}_{k=0} frac{x^k}{ Gamma( alpha k + eta ) }$
Matlab内部函数表示：

function u = Mitt_Leff( x,alph,beta,degree)
%MITT_LEFF Mittag-Leffer函数
%   E_alph_beta = limitssum^{n or infty}_{k=0} frac{x^n}{ gamma( alpha n + eta ) }
%   degree:截取的项数（字符型），一般为‘20’即可
%   alph，beta：为Mittag-Leffer函数的参数
%   x：可以为向量或者标量
%   alpha = eta = 1时，Mitt_Leff( x,1,1,degree ) = e^x.
%   Author: Hu Dongdong
n = str2double( degree );
len = length( x );
u = zeros( 1,len );
for i =1:len
    u( i ) = sum( x(i).^linspace( 0,n,n+1 ) ./ gamma( alph * linspace( 0,n,n+1 ) + beta  ) );
end

end