function b=binomimao(n,k)
%BINLMIAL compute binomial coefficiet
%
%Usage:b=binomial(n,k)
%
%Parameters:         (n)
%               b=( )
%                 (k)
%
%Auther: Zhou Weida (Zhouwd@rsp.xidian.edu.cn)
if (nargin ~=2 | k>n|k<0)
%check correct number of arguments
    help binomial
else
    
    b = prod(k+1:n)/prod(1:(n-k);
    
end