\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2.0} + 1.0}{2.0}\begin{array}{l}
\mathbf{if}\;\alpha \le 139510554198006.0:\\
\;\;\;\;{e}^{\left(\log \left(\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0\right)}{2.0}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\sqrt[3]{\beta} \cdot \sqrt[3]{\beta}}{\sqrt[3]{\left(\alpha + \beta\right) + 2.0} \cdot \sqrt[3]{\left(\alpha + \beta\right) + 2.0}} \cdot \frac{\sqrt[3]{\beta}}{\sqrt[3]{\left(\alpha + \beta\right) + 2.0}} - \left(\left(\frac{\frac{4.0}{\alpha}}{\alpha} - \frac{2.0}{\alpha}\right) - \frac{8.0}{\left(\alpha \cdot \alpha\right) \cdot \alpha}\right)}{2.0}\\
\end{array}double f(double alpha, double beta) {
double r1490466 = beta;
double r1490467 = alpha;
double r1490468 = r1490466 - r1490467;
double r1490469 = r1490467 + r1490466;
double r1490470 = 2.0;
double r1490471 = r1490469 + r1490470;
double r1490472 = r1490468 / r1490471;
double r1490473 = 1.0;
double r1490474 = r1490472 + r1490473;
double r1490475 = r1490474 / r1490470;
return r1490475;
}
double f(double alpha, double beta) {
double r1490476 = alpha;
double r1490477 = 139510554198006.0;
bool r1490478 = r1490476 <= r1490477;
double r1490479 = exp(1.0);
double r1490480 = beta;
double r1490481 = r1490476 + r1490480;
double r1490482 = 2.0;
double r1490483 = r1490481 + r1490482;
double r1490484 = r1490480 / r1490483;
double r1490485 = r1490476 / r1490483;
double r1490486 = 1.0;
double r1490487 = r1490485 - r1490486;
double r1490488 = r1490484 - r1490487;
double r1490489 = r1490488 / r1490482;
double r1490490 = log(r1490489);
double r1490491 = pow(r1490479, r1490490);
double r1490492 = cbrt(r1490480);
double r1490493 = r1490492 * r1490492;
double r1490494 = cbrt(r1490483);
double r1490495 = r1490494 * r1490494;
double r1490496 = r1490493 / r1490495;
double r1490497 = r1490492 / r1490494;
double r1490498 = r1490496 * r1490497;
double r1490499 = 4.0;
double r1490500 = r1490499 / r1490476;
double r1490501 = r1490500 / r1490476;
double r1490502 = r1490482 / r1490476;
double r1490503 = r1490501 - r1490502;
double r1490504 = 8.0;
double r1490505 = r1490476 * r1490476;
double r1490506 = r1490505 * r1490476;
double r1490507 = r1490504 / r1490506;
double r1490508 = r1490503 - r1490507;
double r1490509 = r1490498 - r1490508;
double r1490510 = r1490509 / r1490482;
double r1490511 = r1490478 ? r1490491 : r1490510;
return r1490511;
}



Bits error versus alpha



Bits error versus beta
Results
if alpha < 139510554198006.0Initial program 0.3
rmApplied div-sub0.3
Applied associate-+l-0.3
rmApplied add-exp-log0.3
rmApplied pow10.3
Applied log-pow0.3
Applied exp-prod0.3
Simplified0.3
if 139510554198006.0 < alpha Initial program 50.2
rmApplied div-sub50.2
Applied associate-+l-48.6
rmApplied add-cube-cbrt48.7
Applied add-cube-cbrt48.6
Applied times-frac48.6
Taylor expanded around inf 18.4
Simplified18.4
Final simplification6.0
herbie shell --seed 2019154
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/1"
:pre (and (> alpha -1) (> beta -1))
(/ (+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0) 2.0))