\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 -281574713299796459520 \lor \neg \left(t \le 3.422950075620730720749696840768394601184 \cdot 10^{-8}\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) - \log \left(e^{\frac{2}{t \cdot 3}}\right)\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 r108925 = x;
double r108926 = y;
double r108927 = 2.0;
double r108928 = z;
double r108929 = t;
double r108930 = a;
double r108931 = r108929 + r108930;
double r108932 = sqrt(r108931);
double r108933 = r108928 * r108932;
double r108934 = r108933 / r108929;
double r108935 = b;
double r108936 = c;
double r108937 = r108935 - r108936;
double r108938 = 5.0;
double r108939 = 6.0;
double r108940 = r108938 / r108939;
double r108941 = r108930 + r108940;
double r108942 = 3.0;
double r108943 = r108929 * r108942;
double r108944 = r108927 / r108943;
double r108945 = r108941 - r108944;
double r108946 = r108937 * r108945;
double r108947 = r108934 - r108946;
double r108948 = r108927 * r108947;
double r108949 = exp(r108948);
double r108950 = r108926 * r108949;
double r108951 = r108925 + r108950;
double r108952 = r108925 / r108951;
return r108952;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r108953 = t;
double r108954 = -2.8157471329979646e+20;
bool r108955 = r108953 <= r108954;
double r108956 = 3.422950075620731e-08;
bool r108957 = r108953 <= r108956;
double r108958 = !r108957;
bool r108959 = r108955 || r108958;
double r108960 = x;
double r108961 = y;
double r108962 = 2.0;
double r108963 = z;
double r108964 = cbrt(r108953);
double r108965 = r108964 * r108964;
double r108966 = r108963 / r108965;
double r108967 = a;
double r108968 = r108953 + r108967;
double r108969 = sqrt(r108968);
double r108970 = r108969 / r108964;
double r108971 = r108966 * r108970;
double r108972 = b;
double r108973 = c;
double r108974 = r108972 - r108973;
double r108975 = 5.0;
double r108976 = 6.0;
double r108977 = r108975 / r108976;
double r108978 = r108967 + r108977;
double r108979 = 3.0;
double r108980 = r108953 * r108979;
double r108981 = r108962 / r108980;
double r108982 = exp(r108981);
double r108983 = log(r108982);
double r108984 = r108978 - r108983;
double r108985 = r108974 * r108984;
double r108986 = r108971 - r108985;
double r108987 = r108962 * r108986;
double r108988 = exp(r108987);
double r108989 = r108961 * r108988;
double r108990 = r108960 + r108989;
double r108991 = r108960 / r108990;
double r108992 = r108963 * r108970;
double r108993 = r108967 - r108977;
double r108994 = r108993 * r108980;
double r108995 = r108992 * r108994;
double r108996 = r108967 * r108967;
double r108997 = r108977 * r108977;
double r108998 = r108996 - r108997;
double r108999 = r108998 * r108980;
double r109000 = r108993 * r108962;
double r109001 = r108999 - r109000;
double r109002 = r108974 * r109001;
double r109003 = r108965 * r109002;
double r109004 = r108995 - r109003;
double r109005 = r108965 * r108994;
double r109006 = r109004 / r109005;
double r109007 = r108962 * r109006;
double r109008 = exp(r109007);
double r109009 = r108961 * r109008;
double r109010 = r108960 + r109009;
double r109011 = r108960 / r109010;
double r109012 = r108959 ? r108991 : r109011;
return r109012;
}



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.8157471329979646e+20 or 3.422950075620731e-08 < t Initial program 3.1
rmApplied add-cube-cbrt3.1
Applied times-frac0.1
rmApplied add-log-exp0.3
if -2.8157471329979646e+20 < t < 3.422950075620731e-08Initial program 5.5
rmApplied add-cube-cbrt5.5
Applied times-frac5.7
rmApplied flip-+8.3
Applied frac-sub8.4
Applied associate-*r/8.5
Applied associate-*l/8.3
Applied frac-sub5.8
Final simplification3.0
herbie shell --seed 2019322
(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)))))))))))