Returns the combinations based on the gamma function.
double combination(num n, num k) => factorial(n) / factorial(k) / factorial(n - k);