Returns the logarithm of the factorial based on the gammaLn function.
double factorialLn(num n) => n < 0.0 ? double.nan : gammaLn(1.0 + n);