Function statrs::function::gamma::ln_gamma

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