Function beta 

pub fn beta<F>(a: F, b: F) -> F
where
    F: Float,

B(a, b) = Γ(a) Γ(b)/Γ(a + b), the beta function. A pole of Γ(a) or Γ(b) gives NaN, and a pole of Γ(a + b) alone gives 0.