\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 -1.333804143538084395085776415371440122739 \cdot 10^{-42} \lor \neg \left(t \le 1.029580973042603122977889054947443375466 \cdot 10^{-146}\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 r95061 = x;
double r95062 = y;
double r95063 = 2.0;
double r95064 = z;
double r95065 = t;
double r95066 = a;
double r95067 = r95065 + r95066;
double r95068 = sqrt(r95067);
double r95069 = r95064 * r95068;
double r95070 = r95069 / r95065;
double r95071 = b;
double r95072 = c;
double r95073 = r95071 - r95072;
double r95074 = 5.0;
double r95075 = 6.0;
double r95076 = r95074 / r95075;
double r95077 = r95066 + r95076;
double r95078 = 3.0;
double r95079 = r95065 * r95078;
double r95080 = r95063 / r95079;
double r95081 = r95077 - r95080;
double r95082 = r95073 * r95081;
double r95083 = r95070 - r95082;
double r95084 = r95063 * r95083;
double r95085 = exp(r95084);
double r95086 = r95062 * r95085;
double r95087 = r95061 + r95086;
double r95088 = r95061 / r95087;
return r95088;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r95089 = t;
double r95090 = -1.3338041435380844e-42;
bool r95091 = r95089 <= r95090;
double r95092 = 1.0295809730426031e-146;
bool r95093 = r95089 <= r95092;
double r95094 = !r95093;
bool r95095 = r95091 || r95094;
double r95096 = x;
double r95097 = y;
double r95098 = 2.0;
double r95099 = z;
double r95100 = cbrt(r95089);
double r95101 = r95100 * r95100;
double r95102 = r95099 / r95101;
double r95103 = a;
double r95104 = r95089 + r95103;
double r95105 = sqrt(r95104);
double r95106 = r95105 / r95100;
double r95107 = r95102 * r95106;
double r95108 = b;
double r95109 = c;
double r95110 = r95108 - r95109;
double r95111 = 5.0;
double r95112 = 6.0;
double r95113 = r95111 / r95112;
double r95114 = r95103 + r95113;
double r95115 = 3.0;
double r95116 = r95089 * r95115;
double r95117 = r95098 / r95116;
double r95118 = r95114 - r95117;
double r95119 = r95110 * r95118;
double r95120 = r95107 - r95119;
double r95121 = r95098 * r95120;
double r95122 = exp(r95121);
double r95123 = r95097 * r95122;
double r95124 = r95096 + r95123;
double r95125 = r95096 / r95124;
double r95126 = r95099 * r95106;
double r95127 = r95103 - r95113;
double r95128 = r95127 * r95116;
double r95129 = r95126 * r95128;
double r95130 = r95103 * r95103;
double r95131 = r95113 * r95113;
double r95132 = r95130 - r95131;
double r95133 = r95132 * r95116;
double r95134 = r95127 * r95098;
double r95135 = r95133 - r95134;
double r95136 = r95110 * r95135;
double r95137 = r95101 * r95136;
double r95138 = r95129 - r95137;
double r95139 = r95101 * r95128;
double r95140 = r95138 / r95139;
double r95141 = r95098 * r95140;
double r95142 = exp(r95141);
double r95143 = r95097 * r95142;
double r95144 = r95096 + r95143;
double r95145 = r95096 / r95144;
double r95146 = r95095 ? r95125 : r95145;
return r95146;
}



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 < -1.3338041435380844e-42 or 1.0295809730426031e-146 < t Initial program 2.3
rmApplied add-cube-cbrt2.3
Applied times-frac0.6
if -1.3338041435380844e-42 < t < 1.0295809730426031e-146Initial program 7.0
rmApplied add-cube-cbrt7.0
Applied times-frac7.2
rmApplied flip-+11.2
Applied frac-sub11.2
Applied associate-*r/11.2
Applied associate-*l/10.9
Applied frac-sub8.0
Final simplification2.9
herbie shell --seed 2019235
(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)))))))))))