factbygama.m

来自「Matlab语言常用算法程序集.rar 里面有很多有用的Matlab数值运」· M 代码 · 共 11 行

function y = factbygama(n)
format long;
y = 1;
if n <= 32 
    for i=1:n
        y = y*i;
    end
else
    y = exp(lngama(n+1));
end