pub fn gamma(z: f64) -> f64
Gamma function Γ(z) via the Lanczos approximation, with the reflection formula for z < 0.5.