binom.m

来自「地震、测井方面matlab代码,解释的比较详细」· M 代码 · 共 18 行

function bin = binom(a,n,sign)
% Function computes n binomial coefficients
%      1 sign*a (sign^2)*a*(a-1)/2 (sign^3)*a*(a-1)*(a-2)/3! ...
% If a is a positive integer n can be dropped (it is then set to a+1)
%	bin = binom(a,n)
%
if exist('n') ~= 1; n=a+1; end
if exist('sign') ~= 1; sign=1; end
if n < 1 ; fprintf(' ERROR in binom; n < 1 \n'); return; end
bin=ones(n,1);
ai=a;
for i=1:n-1
  bin(i+1)=sign*bin(i)*ai/i;
  ai=ai-1;
end