pub fn ln_gamma(x: f64) -> f64
Natural log of the Gamma function, ln Γ(x), via the Lanczos approximation (g = 7, 9 coefficients) with the reflection formula for x < 0.5. Accurate to ~15 significant figures for x > 0.
ln Γ(x)
x < 0.5
x > 0