\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\frac{x \cdot \left(\left(\sqrt[3]{e^{\left(\log a \cdot \left(t - 1\right) + \log z \cdot y\right) - b}} \cdot \sqrt[3]{e^{\left(\log a \cdot \left(t - 1\right) + \log z \cdot y\right) - b}}\right) \cdot \sqrt[3]{\left(\sqrt[3]{e^{\left(\log a \cdot \left(t - 1\right) + \log z \cdot y\right) - b}} \cdot \sqrt[3]{e^{\left(\log a \cdot \left(t - 1\right) + \log z \cdot y\right) - b}}\right) \cdot \sqrt[3]{{\left(\sqrt{e}\right)}^{\left(\left(\log a \cdot \left(t - 1\right) + \log z \cdot y\right) - b\right)} \cdot {\left(\sqrt{e}\right)}^{\left(\left(\log a \cdot \left(t - 1\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 r2966722 = x;
double r2966723 = y;
double r2966724 = z;
double r2966725 = log(r2966724);
double r2966726 = r2966723 * r2966725;
double r2966727 = t;
double r2966728 = 1.0;
double r2966729 = r2966727 - r2966728;
double r2966730 = a;
double r2966731 = log(r2966730);
double r2966732 = r2966729 * r2966731;
double r2966733 = r2966726 + r2966732;
double r2966734 = b;
double r2966735 = r2966733 - r2966734;
double r2966736 = exp(r2966735);
double r2966737 = r2966722 * r2966736;
double r2966738 = r2966737 / r2966723;
return r2966738;
}
double f(double x, double y, double z, double t, double a, double b) {
double r2966739 = x;
double r2966740 = a;
double r2966741 = log(r2966740);
double r2966742 = t;
double r2966743 = 1.0;
double r2966744 = r2966742 - r2966743;
double r2966745 = r2966741 * r2966744;
double r2966746 = z;
double r2966747 = log(r2966746);
double r2966748 = y;
double r2966749 = r2966747 * r2966748;
double r2966750 = r2966745 + r2966749;
double r2966751 = b;
double r2966752 = r2966750 - r2966751;
double r2966753 = exp(r2966752);
double r2966754 = cbrt(r2966753);
double r2966755 = r2966754 * r2966754;
double r2966756 = exp(1.0);
double r2966757 = sqrt(r2966756);
double r2966758 = pow(r2966757, r2966752);
double r2966759 = r2966758 * r2966758;
double r2966760 = cbrt(r2966759);
double r2966761 = r2966755 * r2966760;
double r2966762 = cbrt(r2966761);
double r2966763 = r2966755 * r2966762;
double r2966764 = r2966739 * r2966763;
double r2966765 = r2966764 / r2966748;
return r2966765;
}



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