\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 -7.929621432689752 \cdot 10^{-232} \lor \neg \left(t \le 3.2534287537046085 \cdot 10^{-299}\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 r79953 = x;
double r79954 = y;
double r79955 = 2.0;
double r79956 = z;
double r79957 = t;
double r79958 = a;
double r79959 = r79957 + r79958;
double r79960 = sqrt(r79959);
double r79961 = r79956 * r79960;
double r79962 = r79961 / r79957;
double r79963 = b;
double r79964 = c;
double r79965 = r79963 - r79964;
double r79966 = 5.0;
double r79967 = 6.0;
double r79968 = r79966 / r79967;
double r79969 = r79958 + r79968;
double r79970 = 3.0;
double r79971 = r79957 * r79970;
double r79972 = r79955 / r79971;
double r79973 = r79969 - r79972;
double r79974 = r79965 * r79973;
double r79975 = r79962 - r79974;
double r79976 = r79955 * r79975;
double r79977 = exp(r79976);
double r79978 = r79954 * r79977;
double r79979 = r79953 + r79978;
double r79980 = r79953 / r79979;
return r79980;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r79981 = t;
double r79982 = -7.929621432689752e-232;
bool r79983 = r79981 <= r79982;
double r79984 = 3.2534287537046085e-299;
bool r79985 = r79981 <= r79984;
double r79986 = !r79985;
bool r79987 = r79983 || r79986;
double r79988 = x;
double r79989 = y;
double r79990 = 2.0;
double r79991 = z;
double r79992 = cbrt(r79981);
double r79993 = r79992 * r79992;
double r79994 = r79991 / r79993;
double r79995 = a;
double r79996 = r79981 + r79995;
double r79997 = sqrt(r79996);
double r79998 = r79997 / r79992;
double r79999 = r79994 * r79998;
double r80000 = b;
double r80001 = c;
double r80002 = r80000 - r80001;
double r80003 = 5.0;
double r80004 = 6.0;
double r80005 = r80003 / r80004;
double r80006 = r79995 + r80005;
double r80007 = 3.0;
double r80008 = r79981 * r80007;
double r80009 = r79990 / r80008;
double r80010 = r80006 - r80009;
double r80011 = r80002 * r80010;
double r80012 = r79999 - r80011;
double r80013 = r79990 * r80012;
double r80014 = exp(r80013);
double r80015 = r79989 * r80014;
double r80016 = r79988 + r80015;
double r80017 = r79988 / r80016;
double r80018 = r79991 * r79998;
double r80019 = r79995 - r80005;
double r80020 = r80019 * r80008;
double r80021 = r80018 * r80020;
double r80022 = r79995 * r79995;
double r80023 = r80005 * r80005;
double r80024 = r80022 - r80023;
double r80025 = r80024 * r80008;
double r80026 = r80019 * r79990;
double r80027 = r80025 - r80026;
double r80028 = r80002 * r80027;
double r80029 = r79993 * r80028;
double r80030 = r80021 - r80029;
double r80031 = r79993 * r80020;
double r80032 = r80030 / r80031;
double r80033 = r79990 * r80032;
double r80034 = exp(r80033);
double r80035 = r79989 * r80034;
double r80036 = r79988 + r80035;
double r80037 = r79988 / r80036;
double r80038 = r79987 ? r80017 : r80037;
return r80038;
}



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 < -7.929621432689752e-232 or 3.2534287537046085e-299 < t Initial program 3.3
rmApplied add-cube-cbrt3.3
Applied times-frac2.1
if -7.929621432689752e-232 < t < 3.2534287537046085e-299Initial program 12.0
rmApplied add-cube-cbrt12.0
Applied times-frac12.1
rmApplied flip-+15.1
Applied frac-sub15.1
Applied associate-*r/15.1
Applied associate-*l/15.0
Applied frac-sub8.9
Final simplification2.6
herbie shell --seed 2020035
(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)))))))))))