Mittag-Leffer函数, Matlab内部函数

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

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