pub fn gamma(x: f32) -> f32