wolffd@0: function c = choose(n,k)
wolffd@0: % CHOOSE The number of ways of choosing k things from n 
wolffd@0: % c = choose(n,k)
wolffd@0: 
wolffd@0: c = factorial(n)/(factorial(k) * factorial(n-k));