pub fn beta_e(a: f64, b: f64) -> (Value, Result)
This routine computes the Beta Function, B(a,b) = \Gamma(a)\Gamma(b)/\Gamma(a+b) subject to a and b not being negative integers.
