pub fn lnbeta(a: f64, b: f64) -> f64
This routine computes the logarithm of the Beta Function, \log(B(a,b)) subject to a and b not being negative integers.