\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2} + 1}{2}\begin{array}{l}
\mathbf{if}\;\alpha \le 92108.895980526984:\\
\;\;\;\;\frac{\sqrt[3]{{\left(\frac{\beta}{\left(\alpha + \beta\right) + 2}\right)}^{3}} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2} - 1\right)}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{\frac{\left(\alpha + \beta\right) + 2}{\beta}} - \left(4 \cdot \frac{1}{{\alpha}^{2}} - \left(2 \cdot \frac{1}{\alpha} + 8 \cdot \frac{1}{{\alpha}^{3}}\right)\right)}{2}\\
\end{array}double f(double alpha, double beta) {
double r105068 = beta;
double r105069 = alpha;
double r105070 = r105068 - r105069;
double r105071 = r105069 + r105068;
double r105072 = 2.0;
double r105073 = r105071 + r105072;
double r105074 = r105070 / r105073;
double r105075 = 1.0;
double r105076 = r105074 + r105075;
double r105077 = r105076 / r105072;
return r105077;
}
double f(double alpha, double beta) {
double r105078 = alpha;
double r105079 = 92108.89598052698;
bool r105080 = r105078 <= r105079;
double r105081 = beta;
double r105082 = r105078 + r105081;
double r105083 = 2.0;
double r105084 = r105082 + r105083;
double r105085 = r105081 / r105084;
double r105086 = 3.0;
double r105087 = pow(r105085, r105086);
double r105088 = cbrt(r105087);
double r105089 = r105078 / r105084;
double r105090 = 1.0;
double r105091 = r105089 - r105090;
double r105092 = r105088 - r105091;
double r105093 = r105092 / r105083;
double r105094 = 1.0;
double r105095 = r105084 / r105081;
double r105096 = r105094 / r105095;
double r105097 = 4.0;
double r105098 = 2.0;
double r105099 = pow(r105078, r105098);
double r105100 = r105094 / r105099;
double r105101 = r105097 * r105100;
double r105102 = r105094 / r105078;
double r105103 = r105083 * r105102;
double r105104 = 8.0;
double r105105 = pow(r105078, r105086);
double r105106 = r105094 / r105105;
double r105107 = r105104 * r105106;
double r105108 = r105103 + r105107;
double r105109 = r105101 - r105108;
double r105110 = r105096 - r105109;
double r105111 = r105110 / r105083;
double r105112 = r105080 ? r105093 : r105111;
return r105112;
}



Bits error versus alpha



Bits error versus beta
Results
if alpha < 92108.89598052698Initial program 0.0
rmApplied div-sub0.0
Applied associate-+l-0.0
rmApplied add-cbrt-cube11.7
Applied add-cbrt-cube14.4
Applied cbrt-undiv14.4
Simplified0.0
if 92108.89598052698 < alpha Initial program 50.0
rmApplied div-sub50.0
Applied associate-+l-48.3
rmApplied clear-num48.3
Taylor expanded around inf 18.3
Final simplification5.9
herbie shell --seed 2020060
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/1"
:precision binary64
:pre (and (> alpha -1) (> beta -1))
(/ (+ (/ (- beta alpha) (+ (+ alpha beta) 2)) 1) 2))