\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.97649787856798446 \cdot 10^{-84} \lor \neg \left(t \le -4.35171199149127813 \cdot 10^{-269}\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 r90595 = x;
double r90596 = y;
double r90597 = 2.0;
double r90598 = z;
double r90599 = t;
double r90600 = a;
double r90601 = r90599 + r90600;
double r90602 = sqrt(r90601);
double r90603 = r90598 * r90602;
double r90604 = r90603 / r90599;
double r90605 = b;
double r90606 = c;
double r90607 = r90605 - r90606;
double r90608 = 5.0;
double r90609 = 6.0;
double r90610 = r90608 / r90609;
double r90611 = r90600 + r90610;
double r90612 = 3.0;
double r90613 = r90599 * r90612;
double r90614 = r90597 / r90613;
double r90615 = r90611 - r90614;
double r90616 = r90607 * r90615;
double r90617 = r90604 - r90616;
double r90618 = r90597 * r90617;
double r90619 = exp(r90618);
double r90620 = r90596 * r90619;
double r90621 = r90595 + r90620;
double r90622 = r90595 / r90621;
return r90622;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r90623 = t;
double r90624 = -1.9764978785679845e-84;
bool r90625 = r90623 <= r90624;
double r90626 = -4.351711991491278e-269;
bool r90627 = r90623 <= r90626;
double r90628 = !r90627;
bool r90629 = r90625 || r90628;
double r90630 = x;
double r90631 = y;
double r90632 = 2.0;
double r90633 = z;
double r90634 = cbrt(r90623);
double r90635 = r90634 * r90634;
double r90636 = r90633 / r90635;
double r90637 = a;
double r90638 = r90623 + r90637;
double r90639 = sqrt(r90638);
double r90640 = r90639 / r90634;
double r90641 = r90636 * r90640;
double r90642 = b;
double r90643 = c;
double r90644 = r90642 - r90643;
double r90645 = 5.0;
double r90646 = 6.0;
double r90647 = r90645 / r90646;
double r90648 = r90637 + r90647;
double r90649 = 3.0;
double r90650 = r90623 * r90649;
double r90651 = r90632 / r90650;
double r90652 = r90648 - r90651;
double r90653 = r90644 * r90652;
double r90654 = r90641 - r90653;
double r90655 = r90632 * r90654;
double r90656 = exp(r90655);
double r90657 = r90631 * r90656;
double r90658 = r90630 + r90657;
double r90659 = r90630 / r90658;
double r90660 = r90633 * r90640;
double r90661 = r90637 - r90647;
double r90662 = r90661 * r90650;
double r90663 = r90660 * r90662;
double r90664 = r90637 * r90637;
double r90665 = r90647 * r90647;
double r90666 = r90664 - r90665;
double r90667 = r90666 * r90650;
double r90668 = r90661 * r90632;
double r90669 = r90667 - r90668;
double r90670 = r90644 * r90669;
double r90671 = r90635 * r90670;
double r90672 = r90663 - r90671;
double r90673 = r90635 * r90662;
double r90674 = r90672 / r90673;
double r90675 = r90632 * r90674;
double r90676 = exp(r90675);
double r90677 = r90631 * r90676;
double r90678 = r90630 + r90677;
double r90679 = r90630 / r90678;
double r90680 = r90629 ? r90659 : r90679;
return r90680;
}



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.9764978785679845e-84 or -4.351711991491278e-269 < t Initial program 3.7
rmApplied add-cube-cbrt3.7
Applied times-frac2.2
if -1.9764978785679845e-84 < t < -4.351711991491278e-269Initial program 5.9
rmApplied add-cube-cbrt5.9
Applied times-frac6.0
rmApplied flip-+8.8
Applied frac-sub8.8
Applied associate-*r/8.8
Applied associate-*l/8.7
Applied frac-sub6.0
Final simplification2.7
herbie shell --seed 2020047
(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)))))))))))