/* nCr */
new; cls;
print nCr(5,5);
print nCr(5,4);
print nCr(5,3);
print nCr(5,2);
print nCr(5,1);
print nCr(5,0);



proc nCr(n,r);
   if r==0 or n==r;
      retp(1);
   endif;
   retp( round( exp(sumc(ln(seqa(n,-1,n)))-sumc(ln(seqa(r,-1,r)))-sumc(ln(seqa(n-r,-1,n-r)))) ) );
endp;