\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 r77987 = x;
double r77988 = y;
double r77989 = z;
double r77990 = log(r77989);
double r77991 = r77988 * r77990;
double r77992 = t;
double r77993 = 1.0;
double r77994 = r77992 - r77993;
double r77995 = a;
double r77996 = log(r77995);
double r77997 = r77994 * r77996;
double r77998 = r77991 + r77997;
double r77999 = b;
double r78000 = r77998 - r77999;
double r78001 = exp(r78000);
double r78002 = r77987 * r78001;
double r78003 = r78002 / r77988;
return r78003;
}
double f(double x, double y, double z, double t, double a, double b) {
double r78004 = x;
double r78005 = y;
double r78006 = z;
double r78007 = log(r78006);
double r78008 = r78005 * r78007;
double r78009 = t;
double r78010 = 1.0;
double r78011 = r78009 - r78010;
double r78012 = a;
double r78013 = log(r78012);
double r78014 = r78011 * r78013;
double r78015 = r78008 + r78014;
double r78016 = b;
double r78017 = r78015 - r78016;
double r78018 = exp(r78017);
double r78019 = r78004 * r78018;
double r78020 = r78019 / r78005;
return r78020;
}



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