function r = binomcoeff2(n,k)
   prod((n-k+1:n)./(1:k))
end;