\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\left(\sqrt[3]{\frac{x \cdot \frac{{\left(\frac{1}{a}\right)}^{1}}{e^{y \cdot \log \left(\frac{1}{z}\right) + \left(\log \left(\frac{1}{a}\right) \cdot t + b\right)}}}{y}} \cdot \sqrt[3]{\left(\sqrt[3]{\frac{x \cdot \frac{{\left(\frac{1}{a}\right)}^{1}}{e^{y \cdot \log \left(\frac{1}{z}\right) + \left(\log \left(\frac{1}{a}\right) \cdot t + b\right)}}}{y}} \cdot \sqrt[3]{\frac{x \cdot \frac{{\left(\frac{1}{a}\right)}^{1}}{e^{y \cdot \log \left(\frac{1}{z}\right) + \left(\log \left(\frac{1}{a}\right) \cdot t + b\right)}}}{y}}\right) \cdot \sqrt[3]{\frac{x \cdot \frac{{\left(\frac{1}{a}\right)}^{1}}{e^{y \cdot \log \left(\frac{1}{z}\right) + \left(\log \left(\frac{1}{a}\right) \cdot t + b\right)}}}{y}}}\right) \cdot \sqrt[3]{\frac{x \cdot \frac{{\left(\frac{1}{a}\right)}^{1}}{e^{y \cdot \log \left(\frac{1}{z}\right) + \left(\log \left(\frac{1}{a}\right) \cdot t + b\right)}}}{y}}double f(double x, double y, double z, double t, double a, double b) {
double r90043 = x;
double r90044 = y;
double r90045 = z;
double r90046 = log(r90045);
double r90047 = r90044 * r90046;
double r90048 = t;
double r90049 = 1.0;
double r90050 = r90048 - r90049;
double r90051 = a;
double r90052 = log(r90051);
double r90053 = r90050 * r90052;
double r90054 = r90047 + r90053;
double r90055 = b;
double r90056 = r90054 - r90055;
double r90057 = exp(r90056);
double r90058 = r90043 * r90057;
double r90059 = r90058 / r90044;
return r90059;
}
double f(double x, double y, double z, double t, double a, double b) {
double r90060 = x;
double r90061 = 1.0;
double r90062 = a;
double r90063 = r90061 / r90062;
double r90064 = 1.0;
double r90065 = pow(r90063, r90064);
double r90066 = y;
double r90067 = z;
double r90068 = r90061 / r90067;
double r90069 = log(r90068);
double r90070 = r90066 * r90069;
double r90071 = log(r90063);
double r90072 = t;
double r90073 = r90071 * r90072;
double r90074 = b;
double r90075 = r90073 + r90074;
double r90076 = r90070 + r90075;
double r90077 = exp(r90076);
double r90078 = r90065 / r90077;
double r90079 = r90060 * r90078;
double r90080 = r90079 / r90066;
double r90081 = cbrt(r90080);
double r90082 = r90081 * r90081;
double r90083 = r90082 * r90081;
double r90084 = cbrt(r90083);
double r90085 = r90081 * r90084;
double r90086 = r90085 * r90081;
return r90086;
}



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 1.8
Taylor expanded around inf 1.9
Simplified1.2
rmApplied add-cube-cbrt1.3
rmApplied add-cube-cbrt1.3
Final simplification1.3
herbie shell --seed 2020018
(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))