Returns the logarithm of the combinations based on the gammaLn function.
double combinationLn(num n, num k) => factorialLn(n) - factorialLn(k) - factorialLn(n - k);