\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 -4.33536139562560236 \cdot 10^{-281} \lor \neg \left(t \le 2.7731803138233924 \cdot 10^{-105}\right):\\
\;\;\;\;\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) - \sqrt[3]{{\left(\frac{2}{t \cdot 3}\right)}^{3}}\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\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)}}}\\
\end{array}double f(double x, double y, double z, double t, double a, double b, double c) {
double r110737 = x;
double r110738 = y;
double r110739 = 2.0;
double r110740 = z;
double r110741 = t;
double r110742 = a;
double r110743 = r110741 + r110742;
double r110744 = sqrt(r110743);
double r110745 = r110740 * r110744;
double r110746 = r110745 / r110741;
double r110747 = b;
double r110748 = c;
double r110749 = r110747 - r110748;
double r110750 = 5.0;
double r110751 = 6.0;
double r110752 = r110750 / r110751;
double r110753 = r110742 + r110752;
double r110754 = 3.0;
double r110755 = r110741 * r110754;
double r110756 = r110739 / r110755;
double r110757 = r110753 - r110756;
double r110758 = r110749 * r110757;
double r110759 = r110746 - r110758;
double r110760 = r110739 * r110759;
double r110761 = exp(r110760);
double r110762 = r110738 * r110761;
double r110763 = r110737 + r110762;
double r110764 = r110737 / r110763;
return r110764;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r110765 = t;
double r110766 = -4.3353613956256024e-281;
bool r110767 = r110765 <= r110766;
double r110768 = 2.7731803138233924e-105;
bool r110769 = r110765 <= r110768;
double r110770 = !r110769;
bool r110771 = r110767 || r110770;
double r110772 = x;
double r110773 = y;
double r110774 = 2.0;
double r110775 = z;
double r110776 = a;
double r110777 = r110765 + r110776;
double r110778 = sqrt(r110777);
double r110779 = r110775 * r110778;
double r110780 = r110779 / r110765;
double r110781 = b;
double r110782 = c;
double r110783 = r110781 - r110782;
double r110784 = 5.0;
double r110785 = 6.0;
double r110786 = r110784 / r110785;
double r110787 = r110776 + r110786;
double r110788 = 3.0;
double r110789 = r110765 * r110788;
double r110790 = r110774 / r110789;
double r110791 = 3.0;
double r110792 = pow(r110790, r110791);
double r110793 = cbrt(r110792);
double r110794 = r110787 - r110793;
double r110795 = r110783 * r110794;
double r110796 = r110780 - r110795;
double r110797 = r110774 * r110796;
double r110798 = exp(r110797);
double r110799 = r110773 * r110798;
double r110800 = r110772 + r110799;
double r110801 = r110772 / r110800;
double r110802 = r110776 - r110786;
double r110803 = r110802 * r110789;
double r110804 = r110779 * r110803;
double r110805 = r110776 * r110776;
double r110806 = r110786 * r110786;
double r110807 = r110805 - r110806;
double r110808 = r110807 * r110789;
double r110809 = r110802 * r110774;
double r110810 = r110808 - r110809;
double r110811 = r110783 * r110810;
double r110812 = r110765 * r110811;
double r110813 = r110804 - r110812;
double r110814 = r110765 * r110803;
double r110815 = r110813 / r110814;
double r110816 = r110774 * r110815;
double r110817 = exp(r110816);
double r110818 = r110773 * r110817;
double r110819 = r110772 + r110818;
double r110820 = r110772 / r110819;
double r110821 = r110771 ? r110801 : r110820;
return r110821;
}



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 < -4.3353613956256024e-281 or 2.7731803138233924e-105 < t Initial program 3.3
rmApplied add-cbrt-cube3.3
Applied add-cbrt-cube4.4
Applied cbrt-unprod4.4
Applied add-cbrt-cube4.4
Applied cbrt-undiv4.5
Simplified4.5
if -4.3353613956256024e-281 < t < 2.7731803138233924e-105Initial program 6.4
rmApplied flip-+10.2
Applied frac-sub10.2
Applied associate-*r/10.2
Applied frac-sub8.4
Final simplification5.3
herbie shell --seed 2020062
(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)))))))))))