pub fn ln_gamma(x: f64) -> f64
Logarithm Gamma function
Use Lanczos approximation to implement Gamma function ($g=5, n=7$)