\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}}{y}\frac{x \cdot \left(\left(\sqrt[3]{e^{\left(\log a \cdot \left(t - 1.0\right) + \log z \cdot y\right) - b}} \cdot \sqrt[3]{e^{\left(\log a \cdot \left(t - 1.0\right) + \log z \cdot y\right) - b}}\right) \cdot \sqrt[3]{\left(\sqrt[3]{e^{\left(\log a \cdot \left(t - 1.0\right) + \log z \cdot y\right) - b}} \cdot \sqrt[3]{e^{\left(\log a \cdot \left(t - 1.0\right) + \log z \cdot y\right) - b}}\right) \cdot \sqrt[3]{{\left(\sqrt{e}\right)}^{\left(\left(\log a \cdot \left(t - 1.0\right) + \log z \cdot y\right) - b\right)} \cdot {\left(\sqrt{e}\right)}^{\left(\left(\log a \cdot \left(t - 1.0\right) + \log z \cdot y\right) - b\right)}}}\right)}{y}double f(double x, double y, double z, double t, double a, double b) {
double r3185984 = x;
double r3185985 = y;
double r3185986 = z;
double r3185987 = log(r3185986);
double r3185988 = r3185985 * r3185987;
double r3185989 = t;
double r3185990 = 1.0;
double r3185991 = r3185989 - r3185990;
double r3185992 = a;
double r3185993 = log(r3185992);
double r3185994 = r3185991 * r3185993;
double r3185995 = r3185988 + r3185994;
double r3185996 = b;
double r3185997 = r3185995 - r3185996;
double r3185998 = exp(r3185997);
double r3185999 = r3185984 * r3185998;
double r3186000 = r3185999 / r3185985;
return r3186000;
}
double f(double x, double y, double z, double t, double a, double b) {
double r3186001 = x;
double r3186002 = a;
double r3186003 = log(r3186002);
double r3186004 = t;
double r3186005 = 1.0;
double r3186006 = r3186004 - r3186005;
double r3186007 = r3186003 * r3186006;
double r3186008 = z;
double r3186009 = log(r3186008);
double r3186010 = y;
double r3186011 = r3186009 * r3186010;
double r3186012 = r3186007 + r3186011;
double r3186013 = b;
double r3186014 = r3186012 - r3186013;
double r3186015 = exp(r3186014);
double r3186016 = cbrt(r3186015);
double r3186017 = r3186016 * r3186016;
double r3186018 = exp(1.0);
double r3186019 = sqrt(r3186018);
double r3186020 = pow(r3186019, r3186014);
double r3186021 = r3186020 * r3186020;
double r3186022 = cbrt(r3186021);
double r3186023 = r3186017 * r3186022;
double r3186024 = cbrt(r3186023);
double r3186025 = r3186017 * r3186024;
double r3186026 = r3186001 * r3186025;
double r3186027 = r3186026 / r3186010;
return r3186027;
}



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
rmApplied add-cube-cbrt2.1
rmApplied add-cube-cbrt2.1
rmApplied *-un-lft-identity2.1
Applied exp-prod2.1
Simplified2.1
rmApplied add-sqr-sqrt2.0
Applied unpow-prod-down2.0
Final simplification2.0
herbie shell --seed 2019168
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2"
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))