function C=cmb(n,m) 
% cmb(n,m)  m combinations of n 
if min(n)<0 
   error('undefined for n < 0') 
end 
if min(m)<0 
   error('undefined for m < 0') 
end 
if max(m>n)==1 
   error('undefined for m > n') 
end 
C=fact(n)./(fact(n-m).*fact(m));