DESC
DESC copied to clipboard
Make Dommaschk factorials more accurate
n! / m!
n! = gamma(n+1) gammaln(n+1) = log(n!) exp(gammaln(n+1) - gammaln(m+1)) = n!/m!