\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 9.265329209545527167418045692298918612185 \cdot 10^{-280}:\\
\;\;\;\;\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{elif}\;t \le 2.266034278735730289749003584874828813841 \cdot 10^{-170}:\\
\;\;\;\;\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)}}}\\
\mathbf{else}:\\
\;\;\;\;\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)}}\\
\end{array}double f(double x, double y, double z, double t, double a, double b, double c) {
double r6462889 = x;
double r6462890 = y;
double r6462891 = 2.0;
double r6462892 = z;
double r6462893 = t;
double r6462894 = a;
double r6462895 = r6462893 + r6462894;
double r6462896 = sqrt(r6462895);
double r6462897 = r6462892 * r6462896;
double r6462898 = r6462897 / r6462893;
double r6462899 = b;
double r6462900 = c;
double r6462901 = r6462899 - r6462900;
double r6462902 = 5.0;
double r6462903 = 6.0;
double r6462904 = r6462902 / r6462903;
double r6462905 = r6462894 + r6462904;
double r6462906 = 3.0;
double r6462907 = r6462893 * r6462906;
double r6462908 = r6462891 / r6462907;
double r6462909 = r6462905 - r6462908;
double r6462910 = r6462901 * r6462909;
double r6462911 = r6462898 - r6462910;
double r6462912 = r6462891 * r6462911;
double r6462913 = exp(r6462912);
double r6462914 = r6462890 * r6462913;
double r6462915 = r6462889 + r6462914;
double r6462916 = r6462889 / r6462915;
return r6462916;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r6462917 = t;
double r6462918 = 9.265329209545527e-280;
bool r6462919 = r6462917 <= r6462918;
double r6462920 = x;
double r6462921 = y;
double r6462922 = 2.0;
double r6462923 = z;
double r6462924 = cbrt(r6462917);
double r6462925 = r6462924 * r6462924;
double r6462926 = r6462923 / r6462925;
double r6462927 = a;
double r6462928 = r6462917 + r6462927;
double r6462929 = sqrt(r6462928);
double r6462930 = r6462929 / r6462924;
double r6462931 = r6462926 * r6462930;
double r6462932 = b;
double r6462933 = c;
double r6462934 = r6462932 - r6462933;
double r6462935 = 5.0;
double r6462936 = 6.0;
double r6462937 = r6462935 / r6462936;
double r6462938 = r6462927 + r6462937;
double r6462939 = 3.0;
double r6462940 = r6462917 * r6462939;
double r6462941 = r6462922 / r6462940;
double r6462942 = r6462938 - r6462941;
double r6462943 = r6462934 * r6462942;
double r6462944 = r6462931 - r6462943;
double r6462945 = r6462922 * r6462944;
double r6462946 = exp(r6462945);
double r6462947 = r6462921 * r6462946;
double r6462948 = r6462920 + r6462947;
double r6462949 = r6462920 / r6462948;
double r6462950 = 2.2660342787357303e-170;
bool r6462951 = r6462917 <= r6462950;
double r6462952 = r6462923 * r6462930;
double r6462953 = r6462927 - r6462937;
double r6462954 = r6462953 * r6462940;
double r6462955 = r6462952 * r6462954;
double r6462956 = r6462927 * r6462927;
double r6462957 = r6462937 * r6462937;
double r6462958 = r6462956 - r6462957;
double r6462959 = r6462958 * r6462940;
double r6462960 = r6462953 * r6462922;
double r6462961 = r6462959 - r6462960;
double r6462962 = r6462934 * r6462961;
double r6462963 = r6462925 * r6462962;
double r6462964 = r6462955 - r6462963;
double r6462965 = r6462925 * r6462954;
double r6462966 = r6462964 / r6462965;
double r6462967 = r6462922 * r6462966;
double r6462968 = exp(r6462967);
double r6462969 = r6462921 * r6462968;
double r6462970 = r6462920 + r6462969;
double r6462971 = r6462920 / r6462970;
double r6462972 = r6462951 ? r6462971 : r6462949;
double r6462973 = r6462919 ? r6462949 : r6462972;
return r6462973;
}



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 < 9.265329209545527e-280 or 2.2660342787357303e-170 < t Initial program 3.6
rmApplied add-cube-cbrt3.6
Applied times-frac2.3
if 9.265329209545527e-280 < t < 2.2660342787357303e-170Initial program 7.0
rmApplied add-cube-cbrt7.0
Applied times-frac7.3
rmApplied flip-+11.8
Applied frac-sub11.8
Applied associate-*r/11.8
Applied associate-*l/11.5
Applied frac-sub7.3
Final simplification2.9
herbie shell --seed 2019173
(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)))))))))))