\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2} + 1}{2}\begin{array}{l}
\mathbf{if}\;\alpha \le 246471310145705.31:\\
\;\;\;\;\frac{\frac{\beta \cdot \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2} \cdot \frac{\alpha}{\left(\alpha + \beta\right) + 2} + 1 \cdot \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2} + 1\right)\right) - \left(\left(\alpha + \beta\right) + 2\right) \cdot \left(\frac{1}{\left(\left(\alpha + \beta\right) + 2\right) \cdot \left(\left(\alpha + \beta\right) + 2\right)} \cdot \frac{{\alpha}^{3}}{\left(\alpha + \beta\right) + 2} - {1}^{3}\right)}{\left(\left(\alpha + \beta\right) + 2\right) \cdot \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2} \cdot \frac{\alpha}{\left(\alpha + \beta\right) + 2} + 1 \cdot \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2} + 1\right)\right)}}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2} - \left(\frac{4}{\alpha \cdot \alpha} - \left(\frac{8}{{\alpha}^{3}} + \frac{2}{\alpha}\right)\right)}{2}\\
\end{array}double f(double alpha, double beta) {
double r133986 = beta;
double r133987 = alpha;
double r133988 = r133986 - r133987;
double r133989 = r133987 + r133986;
double r133990 = 2.0;
double r133991 = r133989 + r133990;
double r133992 = r133988 / r133991;
double r133993 = 1.0;
double r133994 = r133992 + r133993;
double r133995 = r133994 / r133990;
return r133995;
}
double f(double alpha, double beta) {
double r133996 = alpha;
double r133997 = 246471310145705.3;
bool r133998 = r133996 <= r133997;
double r133999 = beta;
double r134000 = r133996 + r133999;
double r134001 = 2.0;
double r134002 = r134000 + r134001;
double r134003 = r133996 / r134002;
double r134004 = r134003 * r134003;
double r134005 = 1.0;
double r134006 = r134003 + r134005;
double r134007 = r134005 * r134006;
double r134008 = r134004 + r134007;
double r134009 = r133999 * r134008;
double r134010 = 1.0;
double r134011 = r134002 * r134002;
double r134012 = r134010 / r134011;
double r134013 = 3.0;
double r134014 = pow(r133996, r134013);
double r134015 = r134014 / r134002;
double r134016 = r134012 * r134015;
double r134017 = pow(r134005, r134013);
double r134018 = r134016 - r134017;
double r134019 = r134002 * r134018;
double r134020 = r134009 - r134019;
double r134021 = r134002 * r134008;
double r134022 = r134020 / r134021;
double r134023 = r134022 / r134001;
double r134024 = r133999 / r134002;
double r134025 = 4.0;
double r134026 = r133996 * r133996;
double r134027 = r134025 / r134026;
double r134028 = 8.0;
double r134029 = r134028 / r134014;
double r134030 = r134001 / r133996;
double r134031 = r134029 + r134030;
double r134032 = r134027 - r134031;
double r134033 = r134024 - r134032;
double r134034 = r134033 / r134001;
double r134035 = r133998 ? r134023 : r134034;
return r134035;
}



Bits error versus alpha



Bits error versus beta
Results
if alpha < 246471310145705.3Initial program 0.3
rmApplied div-sub0.3
Applied associate-+l-0.3
rmApplied flip3--0.3
Applied frac-sub0.3
Simplified0.3
Simplified0.3
rmApplied add-cube-cbrt0.4
Applied *-un-lft-identity0.4
Applied times-frac0.4
Applied unpow-prod-down0.4
Simplified0.4
Simplified0.3
if 246471310145705.3 < alpha Initial program 50.0
rmApplied div-sub50.0
Applied associate-+l-48.4
Taylor expanded around inf 18.6
Simplified18.6
Final simplification6.1
herbie shell --seed 2019199
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/1"
:pre (and (> alpha -1.0) (> beta -1.0))
(/ (+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0) 2.0))