Mittag-Leffer函数, Matlab内部函数

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

function u = Mitt_Leff( x,alph,beta,degree)
%MITT_LEFF Mittag-Leffer函数
%   E_alph_beta = \limits\sum^{n or \infty}_{k=0} \frac{x^n}{ \gamma( \alpha n + \beta ) }
%   degree:截取的项数（字符型），一般为‘20’即可
%   alph，beta：为Mittag-Leffer函数的参数
%   x：可以为向量或者标量
%   \alpha = \beta = 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