function y=ncomb(a,b)

   y = fact(a)/fact(b)/fact(a-b);
