function x=choose(n,r)

x=prod((n-r+1):n)/factorial(r);