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