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(\left(\log 1 - \left(\frac{1}{2} \cdot \frac{{z}^{2}}{{1}^{2}} + 1 \cdot z\right)\right) - b\right)}}\right) \cdot \sqrt{e^{y \cdot \left(\log z - t\right) + a \cdot \left(\left(\log 1 - \left(\frac{1}{2} \cdot \frac{{z}^{2}}{{1}^{2}} + 1 \cdot z\right)\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) - ((0.5 * (pow(z, 2.0) / pow(1.0, 2.0))) + (1.0 * z))) - b)))))) * sqrt(exp(((y * (log(z) - t)) + (a * ((log(1.0) - ((0.5 * (pow(z, 2.0) / pow(1.0, 2.0))) + (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.0
Taylor expanded around 0 0.5
rmApplied add-sqr-sqrt0.5
Applied associate-*r*0.5
Final simplification0.5
herbie shell --seed 2020100
(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))))))