\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2} + 1}{2}\begin{array}{l}
\mathbf{if}\;\alpha \le 909606014637454464:\\
\;\;\;\;\frac{\frac{\sqrt[3]{\beta} \cdot \sqrt[3]{\beta}}{\sqrt[3]{\left(\alpha + \beta\right) + 2} \cdot \sqrt[3]{\left(\alpha + \beta\right) + 2}} \cdot \frac{\sqrt[3]{\beta}}{\sqrt[3]{\left(\alpha + \beta\right) + 2}} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2} - 1\right)}{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 code(double alpha, double beta) {
return ((((beta - alpha) / ((alpha + beta) + 2.0)) + 1.0) / 2.0);
}
double code(double alpha, double beta) {
double VAR;
if ((alpha <= 9.096060146374545e+17)) {
VAR = (((((cbrt(beta) * cbrt(beta)) / (cbrt(((alpha + beta) + 2.0)) * cbrt(((alpha + beta) + 2.0)))) * (cbrt(beta) / cbrt(((alpha + beta) + 2.0)))) - ((alpha / ((alpha + beta) + 2.0)) - 1.0)) / 2.0);
} else {
VAR = (((beta / ((alpha + beta) + 2.0)) - ((((4.0 / alpha) / alpha) + (-8.0 / pow(alpha, 3.0))) + (-2.0 / alpha))) / 2.0);
}
return VAR;
}



Bits error versus alpha



Bits error versus beta
Results
if alpha < 9.096060146374545e+17Initial program 0.6
rmApplied div-sub0.6
Applied associate-+l-0.6
rmApplied add-cube-cbrt0.8
Applied add-cube-cbrt0.6
Applied times-frac0.6
if 9.096060146374545e+17 < alpha Initial program 50.8
rmApplied div-sub50.8
Applied associate-+l-49.1
Taylor expanded around inf 18.5
Simplified18.5
Final simplification6.3
herbie shell --seed 2020100
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/1"
:precision binary64
:pre (and (> alpha -1) (> beta -1))
(/ (+ (/ (- beta alpha) (+ (+ alpha beta) 2)) 1) 2))