x \cdot e^{y \cdot \left(\log z - t\right) + a \cdot \left(\log \left(1 - z\right) - b\right)}\left(x \cdot \sqrt{e^{y \cdot \left(\log z - t\right) + a \cdot \left(\log \left(1 - z\right) - b\right)}}\right) \cdot \sqrt{e^{y \cdot \left(\log z - t\right) + a \cdot \left(\log \left(1 - z\right) - b\right)}}double code(double x, double y, double z, double t, double a, double b) {
return (x * exp(((y * (log(z) - t)) + (a * (log((1.0 - z)) - b)))));
}
double code(double x, double y, double z, double t, double a, double b) {
return ((x * sqrt(exp(((y * (log(z) - t)) + (a * (log((1.0 - z)) - b)))))) * sqrt(exp(((y * (log(z) - t)) + (a * (log((1.0 - z)) - b))))));
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a



Bits error versus b
Results
Initial program 2.1
rmApplied add-sqr-sqrt2.1
Applied associate-*r*2.1
Final simplification2.1
herbie shell --seed 2020091
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaApprox from math-functions-0.1.5.2, B"
:precision binary64
(* x (exp (+ (* y (- (log z) t)) (* a (- (log (- 1 z)) b))))))