\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}double f(double x, double y, double z, double t, double a, double b) {
double r90565 = x;
double r90566 = y;
double r90567 = z;
double r90568 = log(r90567);
double r90569 = r90566 * r90568;
double r90570 = t;
double r90571 = 1.0;
double r90572 = r90570 - r90571;
double r90573 = a;
double r90574 = log(r90573);
double r90575 = r90572 * r90574;
double r90576 = r90569 + r90575;
double r90577 = b;
double r90578 = r90576 - r90577;
double r90579 = exp(r90578);
double r90580 = r90565 * r90579;
double r90581 = r90580 / r90566;
return r90581;
}
double f(double x, double y, double z, double t, double a, double b) {
double r90582 = x;
double r90583 = y;
double r90584 = z;
double r90585 = log(r90584);
double r90586 = r90583 * r90585;
double r90587 = t;
double r90588 = 1.0;
double r90589 = r90587 - r90588;
double r90590 = a;
double r90591 = log(r90590);
double r90592 = r90589 * r90591;
double r90593 = r90586 + r90592;
double r90594 = b;
double r90595 = r90593 - r90594;
double r90596 = exp(r90595);
double r90597 = r90582 * r90596;
double r90598 = r90597 / r90583;
return r90598;
}



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 1.9
Final simplification1.9
herbie shell --seed 2020046
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2"
:precision binary64
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1) (log a))) b))) y))