\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 r52078 = x;
double r52079 = y;
double r52080 = z;
double r52081 = log(r52080);
double r52082 = r52079 * r52081;
double r52083 = t;
double r52084 = 1.0;
double r52085 = r52083 - r52084;
double r52086 = a;
double r52087 = log(r52086);
double r52088 = r52085 * r52087;
double r52089 = r52082 + r52088;
double r52090 = b;
double r52091 = r52089 - r52090;
double r52092 = exp(r52091);
double r52093 = r52078 * r52092;
double r52094 = r52093 / r52079;
return r52094;
}
double f(double x, double y, double z, double t, double a, double b) {
double r52095 = x;
double r52096 = y;
double r52097 = z;
double r52098 = log(r52097);
double r52099 = r52096 * r52098;
double r52100 = t;
double r52101 = 1.0;
double r52102 = r52100 - r52101;
double r52103 = a;
double r52104 = log(r52103);
double r52105 = r52102 * r52104;
double r52106 = r52099 + r52105;
double r52107 = b;
double r52108 = r52106 - r52107;
double r52109 = exp(r52108);
double r52110 = r52095 * r52109;
double r52111 = r52110 / r52096;
return r52111;
}



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
(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))