pub fn ln_gamma(x: f64) -> f64
Computes the logarithm of the gamma function with an accuracy of 16 floating point digits. The implementation is derived from “An Analysis of the Lanczos Gamma Approximation”, Glendon Ralph Pugh, 2004 p. 116