\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2} + 1}{2}\begin{array}{l}
\mathbf{if}\;\alpha \le 45468567748.04238128662109375:\\
\;\;\;\;\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2} - \frac{\sqrt[3]{\frac{\alpha}{\left(\alpha + \beta\right) + 2} \cdot \frac{\alpha}{\left(\alpha + \beta\right) + 2} - 1 \cdot 1}}{\sqrt{\frac{\alpha}{\left(\alpha + \beta\right) + 2} + 1}} \cdot \frac{\sqrt[3]{\frac{\alpha}{\left(\alpha + \beta\right) + 2} \cdot \frac{\alpha}{\left(\alpha + \beta\right) + 2} - 1 \cdot 1} \cdot \sqrt[3]{\frac{\alpha}{\left(\alpha + \beta\right) + 2} \cdot \frac{\alpha}{\left(\alpha + \beta\right) + 2} - 1 \cdot 1}}{\sqrt{\frac{\alpha}{\left(\alpha + \beta\right) + 2} + 1}}}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2} - \left(\left(\frac{\frac{4}{\alpha}}{\alpha} + \frac{-8}{{\alpha}^{3}}\right) + \frac{-2}{\alpha}\right)}{2}\\
\end{array}double f(double alpha, double beta) {
double r91994 = beta;
double r91995 = alpha;
double r91996 = r91994 - r91995;
double r91997 = r91995 + r91994;
double r91998 = 2.0;
double r91999 = r91997 + r91998;
double r92000 = r91996 / r91999;
double r92001 = 1.0;
double r92002 = r92000 + r92001;
double r92003 = r92002 / r91998;
return r92003;
}
double f(double alpha, double beta) {
double r92004 = alpha;
double r92005 = 45468567748.04238;
bool r92006 = r92004 <= r92005;
double r92007 = beta;
double r92008 = r92004 + r92007;
double r92009 = 2.0;
double r92010 = r92008 + r92009;
double r92011 = r92007 / r92010;
double r92012 = r92004 / r92010;
double r92013 = r92012 * r92012;
double r92014 = 1.0;
double r92015 = r92014 * r92014;
double r92016 = r92013 - r92015;
double r92017 = cbrt(r92016);
double r92018 = r92012 + r92014;
double r92019 = sqrt(r92018);
double r92020 = r92017 / r92019;
double r92021 = r92017 * r92017;
double r92022 = r92021 / r92019;
double r92023 = r92020 * r92022;
double r92024 = r92011 - r92023;
double r92025 = r92024 / r92009;
double r92026 = 4.0;
double r92027 = r92026 / r92004;
double r92028 = r92027 / r92004;
double r92029 = 8.0;
double r92030 = -r92029;
double r92031 = 3.0;
double r92032 = pow(r92004, r92031);
double r92033 = r92030 / r92032;
double r92034 = r92028 + r92033;
double r92035 = -r92009;
double r92036 = r92035 / r92004;
double r92037 = r92034 + r92036;
double r92038 = r92011 - r92037;
double r92039 = r92038 / r92009;
double r92040 = r92006 ? r92025 : r92039;
return r92040;
}



Bits error versus alpha



Bits error versus beta
Results
if alpha < 45468567748.04238Initial program 0.2
rmApplied div-sub0.2
Applied associate-+l-0.2
rmApplied add-log-exp0.2
Applied add-log-exp0.2
Applied diff-log0.2
Simplified0.2
rmApplied flip--0.2
rmApplied add-sqr-sqrt0.2
Applied add-cube-cbrt0.2
Applied times-frac0.2
Applied exp-prod0.2
Applied log-pow0.2
Simplified0.2
if 45468567748.04238 < alpha Initial program 49.8
rmApplied div-sub49.8
Applied associate-+l-48.2
Taylor expanded around inf 18.3
Simplified18.3
Final simplification6.0
herbie shell --seed 2019322
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/1"
:precision binary64
:pre (and (> alpha -1) (> beta -1))
(/ (+ (/ (- beta alpha) (+ (+ alpha beta) 2)) 1) 2))