\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 r54683 = x;
double r54684 = y;
double r54685 = z;
double r54686 = log(r54685);
double r54687 = r54684 * r54686;
double r54688 = t;
double r54689 = 1.0;
double r54690 = r54688 - r54689;
double r54691 = a;
double r54692 = log(r54691);
double r54693 = r54690 * r54692;
double r54694 = r54687 + r54693;
double r54695 = b;
double r54696 = r54694 - r54695;
double r54697 = exp(r54696);
double r54698 = r54683 * r54697;
double r54699 = r54698 / r54684;
return r54699;
}
double f(double x, double y, double z, double t, double a, double b) {
double r54700 = x;
double r54701 = y;
double r54702 = z;
double r54703 = log(r54702);
double r54704 = r54701 * r54703;
double r54705 = t;
double r54706 = 1.0;
double r54707 = r54705 - r54706;
double r54708 = a;
double r54709 = log(r54708);
double r54710 = r54707 * r54709;
double r54711 = r54704 + r54710;
double r54712 = b;
double r54713 = r54711 - r54712;
double r54714 = exp(r54713);
double r54715 = r54700 * r54714;
double r54716 = r54715 / r54701;
return r54716;
}



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 2019326 +o rules:numerics
(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))