\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2} + 1}{2}\begin{array}{l}
\mathbf{if}\;\alpha \le 58106.02365827491:\\
\;\;\;\;\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}}{e^{\log \left({\left(\left(\alpha + \beta\right) + 2\right)}^{\frac{1}{3}}\right)}} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2} - 1\right)}{2}\\
\mathbf{else}:\\
\;\;\;\;\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{1}{{\alpha}^{2}} \cdot \left(4 - \frac{8}{\alpha}\right) + \frac{-1}{\alpha} \cdot 2\right)}{2}\\
\end{array}double code(double alpha, double beta) {
return ((double) (((double) (((double) (((double) (beta - alpha)) / ((double) (((double) (alpha + beta)) + 2.0)))) + 1.0)) / 2.0));
}
double code(double alpha, double beta) {
double VAR;
if ((alpha <= 58106.023658274906)) {
VAR = ((double) (((double) (((double) (((double) (((double) (((double) cbrt(beta)) * ((double) cbrt(beta)))) / ((double) (((double) cbrt(((double) (((double) (alpha + beta)) + 2.0)))) * ((double) cbrt(((double) (((double) (alpha + beta)) + 2.0)))))))) * ((double) (((double) cbrt(beta)) / ((double) exp(((double) log(((double) pow(((double) (((double) (alpha + beta)) + 2.0)), 0.3333333333333333)))))))))) - ((double) (((double) (alpha / ((double) (((double) (alpha + beta)) + 2.0)))) - 1.0)))) / 2.0));
} else {
VAR = ((double) (((double) (((double) (((double) (((double) (((double) cbrt(beta)) * ((double) cbrt(beta)))) / ((double) (((double) cbrt(((double) (((double) (alpha + beta)) + 2.0)))) * ((double) cbrt(((double) (((double) (alpha + beta)) + 2.0)))))))) * ((double) (((double) cbrt(beta)) / ((double) cbrt(((double) (((double) (alpha + beta)) + 2.0)))))))) - ((double) (((double) (((double) (1.0 / ((double) pow(alpha, 2.0)))) * ((double) (4.0 - ((double) (8.0 / alpha)))))) + ((double) (((double) (-1.0 / alpha)) * 2.0)))))) / 2.0));
}
return VAR;
}



Bits error versus alpha



Bits error versus beta
Results
if alpha < 58106.02365827491Initial program 0.0
rmApplied div-sub0.0
Applied associate-+l-0.0
rmApplied add-cube-cbrt0.3
Applied add-cube-cbrt0.1
Applied times-frac0.1
rmApplied add-exp-log1.2
rmApplied pow1/31.2
if 58106.02365827491 < alpha Initial program 48.9
rmApplied div-sub48.8
Applied associate-+l-47.3
rmApplied add-cube-cbrt47.5
Applied add-cube-cbrt47.4
Applied times-frac47.4
Taylor expanded around inf 18.4
Simplified18.4
Final simplification6.9
herbie shell --seed 2020161
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/1"
:precision binary64
:pre (and (> alpha -1.0) (> beta -1.0))
(/ (+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0) 2.0))