pub fn gamma<F: Float>(x: F) -> F
Γ(x). Γ(±0) is ±∞, a negative whole x and -∞ give NaN, and Γ(∞) is ∞.