\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.871833520333572462735300804758002311224 \cdot 10^{-43} \lor \neg \left(t \le 1.196097185809923326285406708328933897124 \cdot 10^{-103}\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 r98265 = x;
double r98266 = y;
double r98267 = 2.0;
double r98268 = z;
double r98269 = t;
double r98270 = a;
double r98271 = r98269 + r98270;
double r98272 = sqrt(r98271);
double r98273 = r98268 * r98272;
double r98274 = r98273 / r98269;
double r98275 = b;
double r98276 = c;
double r98277 = r98275 - r98276;
double r98278 = 5.0;
double r98279 = 6.0;
double r98280 = r98278 / r98279;
double r98281 = r98270 + r98280;
double r98282 = 3.0;
double r98283 = r98269 * r98282;
double r98284 = r98267 / r98283;
double r98285 = r98281 - r98284;
double r98286 = r98277 * r98285;
double r98287 = r98274 - r98286;
double r98288 = r98267 * r98287;
double r98289 = exp(r98288);
double r98290 = r98266 * r98289;
double r98291 = r98265 + r98290;
double r98292 = r98265 / r98291;
return r98292;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r98293 = t;
double r98294 = -1.8718335203335725e-43;
bool r98295 = r98293 <= r98294;
double r98296 = 1.1960971858099233e-103;
bool r98297 = r98293 <= r98296;
double r98298 = !r98297;
bool r98299 = r98295 || r98298;
double r98300 = x;
double r98301 = y;
double r98302 = 2.0;
double r98303 = z;
double r98304 = a;
double r98305 = r98293 + r98304;
double r98306 = sqrt(r98305);
double r98307 = r98303 * r98306;
double r98308 = r98307 / r98293;
double r98309 = b;
double r98310 = c;
double r98311 = r98309 - r98310;
double r98312 = 5.0;
double r98313 = 6.0;
double r98314 = r98312 / r98313;
double r98315 = r98304 + r98314;
double r98316 = 3.0;
double r98317 = r98293 * r98316;
double r98318 = r98302 / r98317;
double r98319 = 3.0;
double r98320 = pow(r98318, r98319);
double r98321 = cbrt(r98320);
double r98322 = r98315 - r98321;
double r98323 = r98311 * r98322;
double r98324 = r98308 - r98323;
double r98325 = r98302 * r98324;
double r98326 = exp(r98325);
double r98327 = r98301 * r98326;
double r98328 = r98300 + r98327;
double r98329 = r98300 / r98328;
double r98330 = r98304 - r98314;
double r98331 = r98330 * r98317;
double r98332 = r98307 * r98331;
double r98333 = r98304 * r98304;
double r98334 = r98314 * r98314;
double r98335 = r98333 - r98334;
double r98336 = r98335 * r98317;
double r98337 = r98330 * r98302;
double r98338 = r98336 - r98337;
double r98339 = r98311 * r98338;
double r98340 = r98293 * r98339;
double r98341 = r98332 - r98340;
double r98342 = r98293 * r98331;
double r98343 = r98341 / r98342;
double r98344 = r98302 * r98343;
double r98345 = exp(r98344);
double r98346 = r98301 * r98345;
double r98347 = r98300 + r98346;
double r98348 = r98300 / r98347;
double r98349 = r98299 ? r98329 : r98348;
return r98349;
}



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.8718335203335725e-43 or 1.1960971858099233e-103 < t Initial program 2.3
rmApplied add-cbrt-cube2.3
Applied add-cbrt-cube2.3
Applied cbrt-unprod2.3
Applied add-cbrt-cube2.3
Applied cbrt-undiv2.3
Simplified2.3
if -1.8718335203335725e-43 < t < 1.1960971858099233e-103Initial program 6.9
rmApplied flip-+10.3
Applied frac-sub10.3
Applied associate-*r/10.3
Applied frac-sub7.7
Final simplification4.3
herbie shell --seed 2019362
(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)))))))))))