\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 -9.473153273384903490502724565303245723582 \cdot 10^{-244} \lor \neg \left(t \le 9.548444342446725028979639354139199225475 \cdot 10^{-247}\right):\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z}{\sqrt[3]{t} \cdot \sqrt[3]{t}} \cdot \frac{\sqrt{t + a}}{\sqrt[3]{t}} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \frac{\left(z \cdot \frac{\sqrt{t + a}}{\sqrt[3]{t}}\right) \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right) - \left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right) \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)}{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right) \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right)}}}\\
\end{array}double f(double x, double y, double z, double t, double a, double b, double c) {
double r62718 = x;
double r62719 = y;
double r62720 = 2.0;
double r62721 = z;
double r62722 = t;
double r62723 = a;
double r62724 = r62722 + r62723;
double r62725 = sqrt(r62724);
double r62726 = r62721 * r62725;
double r62727 = r62726 / r62722;
double r62728 = b;
double r62729 = c;
double r62730 = r62728 - r62729;
double r62731 = 5.0;
double r62732 = 6.0;
double r62733 = r62731 / r62732;
double r62734 = r62723 + r62733;
double r62735 = 3.0;
double r62736 = r62722 * r62735;
double r62737 = r62720 / r62736;
double r62738 = r62734 - r62737;
double r62739 = r62730 * r62738;
double r62740 = r62727 - r62739;
double r62741 = r62720 * r62740;
double r62742 = exp(r62741);
double r62743 = r62719 * r62742;
double r62744 = r62718 + r62743;
double r62745 = r62718 / r62744;
return r62745;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r62746 = t;
double r62747 = -9.473153273384903e-244;
bool r62748 = r62746 <= r62747;
double r62749 = 9.548444342446725e-247;
bool r62750 = r62746 <= r62749;
double r62751 = !r62750;
bool r62752 = r62748 || r62751;
double r62753 = x;
double r62754 = y;
double r62755 = 2.0;
double r62756 = z;
double r62757 = cbrt(r62746);
double r62758 = r62757 * r62757;
double r62759 = r62756 / r62758;
double r62760 = a;
double r62761 = r62746 + r62760;
double r62762 = sqrt(r62761);
double r62763 = r62762 / r62757;
double r62764 = r62759 * r62763;
double r62765 = b;
double r62766 = c;
double r62767 = r62765 - r62766;
double r62768 = 5.0;
double r62769 = 6.0;
double r62770 = r62768 / r62769;
double r62771 = r62760 + r62770;
double r62772 = 3.0;
double r62773 = r62746 * r62772;
double r62774 = r62755 / r62773;
double r62775 = r62771 - r62774;
double r62776 = r62767 * r62775;
double r62777 = r62764 - r62776;
double r62778 = r62755 * r62777;
double r62779 = exp(r62778);
double r62780 = r62754 * r62779;
double r62781 = r62753 + r62780;
double r62782 = r62753 / r62781;
double r62783 = r62756 * r62763;
double r62784 = r62760 - r62770;
double r62785 = r62784 * r62773;
double r62786 = r62783 * r62785;
double r62787 = r62760 * r62760;
double r62788 = r62770 * r62770;
double r62789 = r62787 - r62788;
double r62790 = r62789 * r62773;
double r62791 = r62784 * r62755;
double r62792 = r62790 - r62791;
double r62793 = r62767 * r62792;
double r62794 = r62758 * r62793;
double r62795 = r62786 - r62794;
double r62796 = r62758 * r62785;
double r62797 = r62795 / r62796;
double r62798 = r62755 * r62797;
double r62799 = exp(r62798);
double r62800 = r62754 * r62799;
double r62801 = r62753 + r62800;
double r62802 = r62753 / r62801;
double r62803 = r62752 ? r62782 : r62802;
return r62803;
}



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 < -9.473153273384903e-244 or 9.548444342446725e-247 < t Initial program 3.0
rmApplied add-cube-cbrt3.0
Applied times-frac1.8
if -9.473153273384903e-244 < t < 9.548444342446725e-247Initial program 10.4
rmApplied add-cube-cbrt10.4
Applied times-frac10.5
rmApplied flip-+13.9
Applied frac-sub13.9
Applied associate-*r/13.9
Applied associate-*l/13.8
Applied frac-sub7.9
Final simplification2.4
herbie shell --seed 2019323
(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 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5 6)) (/ 2 (* t 3)))))))))))