\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 r60767 = x;
double r60768 = y;
double r60769 = z;
double r60770 = log(r60769);
double r60771 = r60768 * r60770;
double r60772 = t;
double r60773 = 1.0;
double r60774 = r60772 - r60773;
double r60775 = a;
double r60776 = log(r60775);
double r60777 = r60774 * r60776;
double r60778 = r60771 + r60777;
double r60779 = b;
double r60780 = r60778 - r60779;
double r60781 = exp(r60780);
double r60782 = r60767 * r60781;
double r60783 = r60782 / r60768;
return r60783;
}
double f(double x, double y, double z, double t, double a, double b) {
double r60784 = x;
double r60785 = y;
double r60786 = z;
double r60787 = log(r60786);
double r60788 = r60785 * r60787;
double r60789 = t;
double r60790 = 1.0;
double r60791 = r60789 - r60790;
double r60792 = a;
double r60793 = log(r60792);
double r60794 = r60791 * r60793;
double r60795 = r60788 + r60794;
double r60796 = b;
double r60797 = r60795 - r60796;
double r60798 = exp(r60797);
double r60799 = r60784 * r60798;
double r60800 = r60799 / r60785;
return r60800;
}



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
Final simplification2.0
herbie shell --seed 2019303
(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))