\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2} + 1}{2}\begin{array}{l}
\mathbf{if}\;\alpha \le 8.43007903703742446 \cdot 10^{35}:\\
\;\;\;\;\frac{\left(\sqrt[3]{\frac{\beta}{\left(\alpha + \beta\right) + 2}} \cdot \log \left(e^{\sqrt[3]{\frac{\beta}{\left(\alpha + \beta\right) + 2}}}\right)\right) \cdot \log \left(e^{\sqrt[3]{\frac{\beta}{\left(\alpha + \beta\right) + 2}}}\right) - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2} - 1\right)}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(\sqrt[3]{\frac{\beta}{\left(\alpha + \beta\right) + 2}} \cdot \sqrt[3]{\frac{\beta}{\left(\alpha + \beta\right) + 2}}\right) \cdot \sqrt[3]{\frac{\beta}{\left(\alpha + \beta\right) + 2}} - \left(\frac{1}{{\alpha}^{2}} \cdot \left(4 - \frac{8}{\alpha}\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 <= 8.430079037037424e+35)) {
VAR = ((((cbrt((beta / ((alpha + beta) + 2.0))) * log(exp(cbrt((beta / ((alpha + beta) + 2.0)))))) * log(exp(cbrt((beta / ((alpha + beta) + 2.0)))))) - ((alpha / ((alpha + beta) + 2.0)) - 1.0)) / 2.0);
} else {
VAR = ((((cbrt((beta / ((alpha + beta) + 2.0))) * cbrt((beta / ((alpha + beta) + 2.0)))) * cbrt((beta / ((alpha + beta) + 2.0)))) - (((1.0 / pow(alpha, 2.0)) * (4.0 - (8.0 / alpha))) + (-2.0 / alpha))) / 2.0);
}
return VAR;
}



Bits error versus alpha



Bits error versus beta
Results
if alpha < 8.430079037037424e+35Initial program 1.6
rmApplied div-sub1.6
Applied associate-+l-1.6
rmApplied add-cube-cbrt1.6
rmApplied add-log-exp1.6
rmApplied add-log-exp1.6
if 8.430079037037424e+35 < alpha Initial program 51.1
rmApplied div-sub51.1
Applied associate-+l-49.4
rmApplied add-cube-cbrt49.4
Taylor expanded around inf 18.2
Simplified18.2
Final simplification6.4
herbie shell --seed 2020102
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/1"
:precision binary64
:pre (and (> alpha -1) (> beta -1))
(/ (+ (/ (- beta alpha) (+ (+ alpha beta) 2)) 1) 2))