Function gamma 

pub fn gamma<F>(x: F) -> F
where
    F: Float,

Γ(x). Γ(±0) is ±∞, a negative whole x and -∞ give NaN, and Γ(∞) is ∞.