\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 r131742 = x;
double r131743 = y;
double r131744 = z;
double r131745 = log(r131744);
double r131746 = r131743 * r131745;
double r131747 = t;
double r131748 = 1.0;
double r131749 = r131747 - r131748;
double r131750 = a;
double r131751 = log(r131750);
double r131752 = r131749 * r131751;
double r131753 = r131746 + r131752;
double r131754 = b;
double r131755 = r131753 - r131754;
double r131756 = exp(r131755);
double r131757 = r131742 * r131756;
double r131758 = r131757 / r131743;
return r131758;
}
double f(double x, double y, double z, double t, double a, double b) {
double r131759 = x;
double r131760 = y;
double r131761 = z;
double r131762 = log(r131761);
double r131763 = r131760 * r131762;
double r131764 = t;
double r131765 = 1.0;
double r131766 = r131764 - r131765;
double r131767 = a;
double r131768 = log(r131767);
double r131769 = r131766 * r131768;
double r131770 = r131763 + r131769;
double r131771 = b;
double r131772 = r131770 - r131771;
double r131773 = exp(r131772);
double r131774 = r131759 * r131773;
double r131775 = r131774 / r131760;
return r131775;
}



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