function c = choose(n,k)

c = factorial(n)/(factorial(k) * factorial(n-k));
