function nb = Binomial(n,k)
    nb = factorial(n)/(factorial(k)*factorial(n-k));
end