\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}\begin{array}{l}
\mathbf{if}\;t \le -2.462699075631974 \cdot 10^{-37}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(a \cdot \left(c - b\right) + 0.83333333333333337 \cdot c\right)}}\\
\mathbf{elif}\;t \le 8.1421942784803702 \cdot 10^{-12}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \frac{\left(z \cdot \sqrt{t + a}\right) \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right) - t \cdot \left(\left(b - c\right) \cdot \left(\left(a \cdot a - \frac{5}{6} \cdot \frac{5}{6}\right) \cdot \left(t \cdot 3\right) - \left(a - \frac{5}{6}\right) \cdot 2\right)\right)}{t \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right)}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \log \left(e^{\frac{2}{t \cdot 3}}\right)\right)\right)}}\\
\end{array}double code(double x, double y, double z, double t, double a, double b, double c) {
return ((double) (x / ((double) (x + ((double) (y * ((double) exp(((double) (2.0 * ((double) (((double) (((double) (z * ((double) sqrt(((double) (t + a)))))) / t)) - ((double) (((double) (b - c)) * ((double) (((double) (a + ((double) (5.0 / 6.0)))) - ((double) (2.0 / ((double) (t * 3.0))))))))))))))))))));
}
double code(double x, double y, double z, double t, double a, double b, double c) {
double VAR;
if ((t <= -2.462699075631974e-37)) {
VAR = ((double) (x / ((double) (x + ((double) (y * ((double) exp(((double) (2.0 * ((double) (((double) (a * ((double) (c - b)))) + ((double) (0.8333333333333334 * c))))))))))))));
} else {
double VAR_1;
if ((t <= 8.14219427848037e-12)) {
VAR_1 = ((double) (x / ((double) (x + ((double) (y * ((double) exp(((double) (2.0 * ((double) (((double) (((double) (((double) (z * ((double) sqrt(((double) (t + a)))))) * ((double) (((double) (a - ((double) (5.0 / 6.0)))) * ((double) (t * 3.0)))))) - ((double) (t * ((double) (((double) (b - c)) * ((double) (((double) (((double) (((double) (a * a)) - ((double) (((double) (5.0 / 6.0)) * ((double) (5.0 / 6.0)))))) * ((double) (t * 3.0)))) - ((double) (((double) (a - ((double) (5.0 / 6.0)))) * 2.0)))))))))) / ((double) (t * ((double) (((double) (a - ((double) (5.0 / 6.0)))) * ((double) (t * 3.0))))))))))))))))));
} else {
VAR_1 = ((double) (x / ((double) (x + ((double) (y * ((double) exp(((double) (2.0 * ((double) (((double) (((double) (z * ((double) sqrt(((double) (t + a)))))) / t)) - ((double) (((double) (b - c)) * ((double) (((double) (a + ((double) (5.0 / 6.0)))) - ((double) log(((double) exp(((double) (2.0 / ((double) (t * 3.0))))))))))))))))))))))));
}
VAR = VAR_1;
}
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



Bits error versus c
Results
if t < -2.462699075631974e-37Initial program 3.1
Taylor expanded around inf 7.9
rmApplied add-exp-log7.9
Simplified4.2
if -2.462699075631974e-37 < t < 8.14219427848037e-12Initial program 5.5
rmApplied flip-+8.5
Applied frac-sub8.5
Applied associate-*r/8.5
Applied frac-sub6.6
if 8.14219427848037e-12 < t Initial program 2.7
rmApplied add-log-exp3.4
Final simplification5.0
herbie shell --seed 2020130
(FPCore (x y z t a b c)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2"
:precision binary64
(/ x (+ x (* y (exp (* 2.0 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))))))))