\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\begin{array}{l}
\mathbf{if}\;y \le -433455108479.84827 \lor \neg \left(y \le 1.1700160174606165\right):\\
\;\;\;\;\frac{\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)}}}{\frac{y}{x}}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{\frac{{\left(\frac{1}{a}\right)}^{1}}{e^{\log \left(\frac{1}{a}\right) \cdot t + b} \cdot {\left(\frac{1}{z}\right)}^{y}}}{y}\\
\end{array}double code(double x, double y, double z, double t, double a, double b) {
return ((double) (((double) (x * ((double) exp(((double) (((double) (((double) (y * ((double) log(z)))) + ((double) (((double) (t - 1.0)) * ((double) log(a)))))) - b)))))) / y));
}
double code(double x, double y, double z, double t, double a, double b) {
double VAR;
if (((y <= -433455108479.84827) || !(y <= 1.1700160174606165))) {
VAR = ((double) (((double) (((double) pow(((double) (1.0 / a)), 1.0)) / ((double) exp(((double) (((double) (y * ((double) log(((double) (1.0 / z)))))) + ((double) (((double) (((double) log(((double) (1.0 / a)))) * t)) + b)))))))) / ((double) (y / x))));
} else {
VAR = ((double) (x * ((double) (((double) (((double) pow(((double) (1.0 / a)), 1.0)) / ((double) (((double) exp(((double) (((double) (((double) log(((double) (1.0 / a)))) * t)) + b)))) * ((double) pow(((double) (1.0 / z)), y)))))) / y))));
}
return VAR;
}



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
if y < -433455108479.84827 or 1.1700160174606165 < y Initial program 0.0
Taylor expanded around inf 0.0
Simplified0.0
if -433455108479.84827 < y < 1.1700160174606165Initial program 3.7
Taylor expanded around inf 3.7
Simplified2.3
rmApplied *-un-lft-identity2.3
Applied times-frac2.6
Simplified2.6
Simplified3.0
Final simplification1.6
herbie shell --seed 2020122
(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))