\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 r62978 = x;
double r62979 = y;
double r62980 = z;
double r62981 = log(r62980);
double r62982 = r62979 * r62981;
double r62983 = t;
double r62984 = 1.0;
double r62985 = r62983 - r62984;
double r62986 = a;
double r62987 = log(r62986);
double r62988 = r62985 * r62987;
double r62989 = r62982 + r62988;
double r62990 = b;
double r62991 = r62989 - r62990;
double r62992 = exp(r62991);
double r62993 = r62978 * r62992;
double r62994 = r62993 / r62979;
return r62994;
}
double f(double x, double y, double z, double t, double a, double b) {
double r62995 = x;
double r62996 = y;
double r62997 = z;
double r62998 = log(r62997);
double r62999 = r62996 * r62998;
double r63000 = t;
double r63001 = 1.0;
double r63002 = r63000 - r63001;
double r63003 = a;
double r63004 = log(r63003);
double r63005 = r63002 * r63004;
double r63006 = r62999 + r63005;
double r63007 = b;
double r63008 = r63006 - r63007;
double r63009 = exp(r63008);
double r63010 = r62995 * r63009;
double r63011 = r63010 / r62996;
return r63011;
}



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