\frac{\frac{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}\begin{array}{l}
\mathbf{if}\;\alpha \le 2.4427574797641429 \cdot 10^{77}:\\
\;\;\;\;\frac{\frac{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + \left(2 \cdot 1 + 1\right)}\\
\mathbf{elif}\;\alpha \le 6.9212401795544938 \cdot 10^{125}:\\
\;\;\;\;\frac{\frac{\left(2 \cdot \frac{1}{{\alpha}^{2}} + 1\right) - 1 \cdot \frac{1}{\alpha}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + \left(2 \cdot 1 + 1\right)}\\
\mathbf{elif}\;\alpha \le 9.03381831342100199 \cdot 10^{202}:\\
\;\;\;\;\frac{\frac{\frac{\sqrt{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}}{\sqrt{\left(\alpha + \beta\right) + 2 \cdot 1}}}{\sqrt{\left(\alpha + \beta\right) + 2 \cdot 1}}}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}{\sqrt{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}} \cdot \left(\left(\alpha + \beta\right) + 2 \cdot 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\left(2 \cdot \frac{1}{{\alpha}^{2}} + 1\right) - 1 \cdot \frac{1}{\alpha}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + \left(2 \cdot 1 + 1\right)}\\
\end{array}double f(double alpha, double beta) {
double r119446 = alpha;
double r119447 = beta;
double r119448 = r119446 + r119447;
double r119449 = r119447 * r119446;
double r119450 = r119448 + r119449;
double r119451 = 1.0;
double r119452 = r119450 + r119451;
double r119453 = 2.0;
double r119454 = r119453 * r119451;
double r119455 = r119448 + r119454;
double r119456 = r119452 / r119455;
double r119457 = r119456 / r119455;
double r119458 = r119455 + r119451;
double r119459 = r119457 / r119458;
return r119459;
}
double f(double alpha, double beta) {
double r119460 = alpha;
double r119461 = 2.442757479764143e+77;
bool r119462 = r119460 <= r119461;
double r119463 = beta;
double r119464 = r119460 + r119463;
double r119465 = r119463 * r119460;
double r119466 = r119464 + r119465;
double r119467 = 1.0;
double r119468 = r119466 + r119467;
double r119469 = 2.0;
double r119470 = r119469 * r119467;
double r119471 = r119464 + r119470;
double r119472 = r119468 / r119471;
double r119473 = r119472 / r119471;
double r119474 = r119470 + r119467;
double r119475 = r119464 + r119474;
double r119476 = r119473 / r119475;
double r119477 = 6.921240179554494e+125;
bool r119478 = r119460 <= r119477;
double r119479 = 1.0;
double r119480 = 2.0;
double r119481 = pow(r119460, r119480);
double r119482 = r119479 / r119481;
double r119483 = r119469 * r119482;
double r119484 = r119483 + r119479;
double r119485 = r119479 / r119460;
double r119486 = r119467 * r119485;
double r119487 = r119484 - r119486;
double r119488 = r119487 / r119471;
double r119489 = r119488 / r119475;
double r119490 = 9.033818313421002e+202;
bool r119491 = r119460 <= r119490;
double r119492 = sqrt(r119468);
double r119493 = sqrt(r119471);
double r119494 = r119492 / r119493;
double r119495 = r119494 / r119493;
double r119496 = r119471 + r119467;
double r119497 = r119496 / r119492;
double r119498 = r119497 * r119471;
double r119499 = r119495 / r119498;
double r119500 = r119491 ? r119499 : r119489;
double r119501 = r119478 ? r119489 : r119500;
double r119502 = r119462 ? r119476 : r119501;
return r119502;
}



Bits error versus alpha



Bits error versus beta
Results
if alpha < 2.442757479764143e+77Initial program 0.4
rmApplied associate-+l+0.4
if 2.442757479764143e+77 < alpha < 6.921240179554494e+125 or 9.033818313421002e+202 < alpha Initial program 14.7
rmApplied associate-+l+14.7
Taylor expanded around inf 7.9
if 6.921240179554494e+125 < alpha < 9.033818313421002e+202Initial program 10.5
rmApplied add-sqr-sqrt10.6
Applied add-sqr-sqrt10.7
Applied add-sqr-sqrt10.6
Applied times-frac10.6
Applied times-frac10.7
Applied associate-/l*10.7
Simplified10.5
Final simplification2.5
herbie shell --seed 2020018
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/3"
:precision binary64
:pre (and (> alpha -1) (> beta -1))
(/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) (+ (+ (+ alpha beta) (* 2 1)) 1)))