\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\frac{1}{\frac{y}{\left(x \cdot \frac{{\left(\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{\sqrt{a}}\right)}^{1}}{\sqrt[3]{e^{y \cdot \log \left(\frac{1}{z}\right) + \left(\log \left(\frac{1}{a}\right) \cdot t + b\right)}} \cdot \sqrt[3]{e^{y \cdot \log \left(\frac{1}{z}\right) + \left(\log \left(\frac{1}{a}\right) \cdot t + b\right)}}}\right) \cdot \frac{{\left(\frac{\sqrt[3]{1}}{\sqrt{a}}\right)}^{1}}{\sqrt[3]{e^{y \cdot \log \left(\frac{1}{z}\right) + \left(\log \left(\frac{1}{a}\right) \cdot t + b\right)}}}}}double f(double x, double y, double z, double t, double a, double b) {
double r89048 = x;
double r89049 = y;
double r89050 = z;
double r89051 = log(r89050);
double r89052 = r89049 * r89051;
double r89053 = t;
double r89054 = 1.0;
double r89055 = r89053 - r89054;
double r89056 = a;
double r89057 = log(r89056);
double r89058 = r89055 * r89057;
double r89059 = r89052 + r89058;
double r89060 = b;
double r89061 = r89059 - r89060;
double r89062 = exp(r89061);
double r89063 = r89048 * r89062;
double r89064 = r89063 / r89049;
return r89064;
}
double f(double x, double y, double z, double t, double a, double b) {
double r89065 = 1.0;
double r89066 = y;
double r89067 = x;
double r89068 = cbrt(r89065);
double r89069 = r89068 * r89068;
double r89070 = a;
double r89071 = sqrt(r89070);
double r89072 = r89069 / r89071;
double r89073 = 1.0;
double r89074 = pow(r89072, r89073);
double r89075 = z;
double r89076 = r89065 / r89075;
double r89077 = log(r89076);
double r89078 = r89066 * r89077;
double r89079 = r89065 / r89070;
double r89080 = log(r89079);
double r89081 = t;
double r89082 = r89080 * r89081;
double r89083 = b;
double r89084 = r89082 + r89083;
double r89085 = r89078 + r89084;
double r89086 = exp(r89085);
double r89087 = cbrt(r89086);
double r89088 = r89087 * r89087;
double r89089 = r89074 / r89088;
double r89090 = r89067 * r89089;
double r89091 = r89068 / r89071;
double r89092 = pow(r89091, r89073);
double r89093 = r89092 / r89087;
double r89094 = r89090 * r89093;
double r89095 = r89066 / r89094;
double r89096 = r89065 / r89095;
return r89096;
}



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
Taylor expanded around inf 2.0
Simplified1.3
rmApplied add-cube-cbrt1.4
Applied add-sqr-sqrt1.4
Applied add-cube-cbrt1.4
Applied times-frac1.4
Applied unpow-prod-down1.4
Applied times-frac1.4
Applied associate-*r*1.4
rmApplied clear-num1.4
Final simplification1.4
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))