\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}\frac{\frac{\sqrt{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{\left(\alpha + \beta\right) + 2 \cdot 1}}}{\sqrt{\left(\alpha + \beta\right) + 2 \cdot 1}} \cdot \frac{\sqrt{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{\left(\alpha + \beta\right) + 2 \cdot 1}}}{\sqrt{\left(\alpha + \beta\right) + 2 \cdot 1}}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}double f(double alpha, double beta) {
double r99705 = alpha;
double r99706 = beta;
double r99707 = r99705 + r99706;
double r99708 = r99706 * r99705;
double r99709 = r99707 + r99708;
double r99710 = 1.0;
double r99711 = r99709 + r99710;
double r99712 = 2.0;
double r99713 = r99712 * r99710;
double r99714 = r99707 + r99713;
double r99715 = r99711 / r99714;
double r99716 = r99715 / r99714;
double r99717 = r99714 + r99710;
double r99718 = r99716 / r99717;
return r99718;
}
double f(double alpha, double beta) {
double r99719 = alpha;
double r99720 = beta;
double r99721 = r99719 + r99720;
double r99722 = r99720 * r99719;
double r99723 = r99721 + r99722;
double r99724 = 1.0;
double r99725 = r99723 + r99724;
double r99726 = 2.0;
double r99727 = r99726 * r99724;
double r99728 = r99721 + r99727;
double r99729 = r99725 / r99728;
double r99730 = sqrt(r99729);
double r99731 = sqrt(r99728);
double r99732 = r99730 / r99731;
double r99733 = r99732 * r99732;
double r99734 = r99728 + r99724;
double r99735 = r99733 / r99734;
return r99735;
}



Bits error versus alpha



Bits error versus beta
Results
Initial program 3.8
rmApplied add-sqr-sqrt4.2
Applied add-sqr-sqrt3.8
Applied times-frac3.9
Final simplification3.9
herbie shell --seed 2020033
(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)))