pub fn ln_gamma_approx(z: f64) -> f64