\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.4919779066772806 \cdot 10^{-159} \lor \neg \left(t \le 2.8294745754139501 \cdot 10^{-110}\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(\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right) \cdot \left(b - c\right)\right) \cdot \left(\left(a + \frac{5}{6}\right) \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right) - \left(a - \frac{5}{6}\right) \cdot 2\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 r189110 = x;
double r189111 = y;
double r189112 = 2.0;
double r189113 = z;
double r189114 = t;
double r189115 = a;
double r189116 = r189114 + r189115;
double r189117 = sqrt(r189116);
double r189118 = r189113 * r189117;
double r189119 = r189118 / r189114;
double r189120 = b;
double r189121 = c;
double r189122 = r189120 - r189121;
double r189123 = 5.0;
double r189124 = 6.0;
double r189125 = r189123 / r189124;
double r189126 = r189115 + r189125;
double r189127 = 3.0;
double r189128 = r189114 * r189127;
double r189129 = r189112 / r189128;
double r189130 = r189126 - r189129;
double r189131 = r189122 * r189130;
double r189132 = r189119 - r189131;
double r189133 = r189112 * r189132;
double r189134 = exp(r189133);
double r189135 = r189111 * r189134;
double r189136 = r189110 + r189135;
double r189137 = r189110 / r189136;
return r189137;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r189138 = t;
double r189139 = -4.4919779066772806e-159;
bool r189140 = r189138 <= r189139;
double r189141 = 2.82947457541395e-110;
bool r189142 = r189138 <= r189141;
double r189143 = !r189142;
bool r189144 = r189140 || r189143;
double r189145 = x;
double r189146 = y;
double r189147 = 2.0;
double r189148 = z;
double r189149 = cbrt(r189138);
double r189150 = r189149 * r189149;
double r189151 = r189148 / r189150;
double r189152 = a;
double r189153 = r189138 + r189152;
double r189154 = sqrt(r189153);
double r189155 = r189154 / r189149;
double r189156 = r189151 * r189155;
double r189157 = b;
double r189158 = c;
double r189159 = r189157 - r189158;
double r189160 = 5.0;
double r189161 = 6.0;
double r189162 = r189160 / r189161;
double r189163 = r189152 + r189162;
double r189164 = 3.0;
double r189165 = r189138 * r189164;
double r189166 = r189147 / r189165;
double r189167 = r189163 - r189166;
double r189168 = r189159 * r189167;
double r189169 = r189156 - r189168;
double r189170 = r189147 * r189169;
double r189171 = exp(r189170);
double r189172 = r189146 * r189171;
double r189173 = r189145 + r189172;
double r189174 = r189145 / r189173;
double r189175 = r189148 * r189155;
double r189176 = r189152 - r189162;
double r189177 = r189176 * r189165;
double r189178 = r189175 * r189177;
double r189179 = r189150 * r189159;
double r189180 = r189163 * r189177;
double r189181 = r189176 * r189147;
double r189182 = r189180 - r189181;
double r189183 = r189179 * r189182;
double r189184 = r189178 - r189183;
double r189185 = r189150 * r189177;
double r189186 = r189184 / r189185;
double r189187 = r189147 * r189186;
double r189188 = exp(r189187);
double r189189 = r189146 * r189188;
double r189190 = r189145 + r189189;
double r189191 = r189145 / r189190;
double r189192 = r189144 ? r189174 : r189191;
return r189192;
}



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.4919779066772806e-159 or 2.82947457541395e-110 < t Initial program 2.6
rmApplied add-cube-cbrt2.6
Applied times-frac0.7
if -4.4919779066772806e-159 < t < 2.82947457541395e-110Initial program 6.8
rmApplied add-cube-cbrt6.8
Applied times-frac7.1
rmApplied flip-+11.0
Applied frac-sub11.0
Applied associate-*r/11.0
Applied associate-*l/10.6
Applied frac-sub7.0
rmApplied difference-of-squares7.0
Applied associate-*l*3.3
rmApplied associate-*r*2.6
Final simplification1.3
herbie shell --seed 2019195
(FPCore (x y z t a b c)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2"
(/ x (+ x (* y (exp (* 2.0 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))))))))