\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 r2986477 = x;
double r2986478 = y;
double r2986479 = z;
double r2986480 = log(r2986479);
double r2986481 = r2986478 * r2986480;
double r2986482 = t;
double r2986483 = 1.0;
double r2986484 = r2986482 - r2986483;
double r2986485 = a;
double r2986486 = log(r2986485);
double r2986487 = r2986484 * r2986486;
double r2986488 = r2986481 + r2986487;
double r2986489 = b;
double r2986490 = r2986488 - r2986489;
double r2986491 = exp(r2986490);
double r2986492 = r2986477 * r2986491;
double r2986493 = r2986492 / r2986478;
return r2986493;
}
double f(double x, double y, double z, double t, double a, double b) {
double r2986494 = x;
double r2986495 = a;
double r2986496 = log(r2986495);
double r2986497 = t;
double r2986498 = 1.0;
double r2986499 = r2986497 - r2986498;
double r2986500 = r2986496 * r2986499;
double r2986501 = z;
double r2986502 = log(r2986501);
double r2986503 = y;
double r2986504 = r2986502 * r2986503;
double r2986505 = r2986500 + r2986504;
double r2986506 = b;
double r2986507 = r2986505 - r2986506;
double r2986508 = exp(r2986507);
double r2986509 = cbrt(r2986508);
double r2986510 = r2986509 * r2986509;
double r2986511 = exp(1.0);
double r2986512 = sqrt(r2986511);
double r2986513 = pow(r2986512, r2986507);
double r2986514 = r2986513 * r2986513;
double r2986515 = cbrt(r2986514);
double r2986516 = r2986510 * r2986515;
double r2986517 = cbrt(r2986516);
double r2986518 = r2986510 * r2986517;
double r2986519 = r2986494 * r2986518;
double r2986520 = r2986519 / r2986503;
return r2986520;
}



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 +o rules:numerics
(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))