bincoef.m

来自「一个非常实用的统计工具箱」· M 代码 · 共 14 行

function k = bincoef(n,N)
%BINCOEF  Binomial coefficients
%
%         k = bincoef(n,N)

%       Anders Holtsberg, 13-05-94
%       Copyright (c) Anders Holtsberg

if all(all(round(n)==n)) & all(all(round(N)==N))
   k = round(exp(gammaln(N+1) - gammaln(n+1) - gammaln(N-n+1)));
else
   k = exp(gammaln(N+1) - gammaln(n+1) - gammaln(N-n+1));
end