\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}\;\beta \le 8.48955521481748054 \cdot 10^{157}:\\
\;\;\;\;\frac{\sqrt{\sqrt{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}}}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}{\sqrt{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{\left(\alpha + \beta\right) + 2 \cdot 1}}}} \cdot \frac{\sqrt{\frac{\frac{\sqrt{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}}{\sqrt{\left(\alpha + \beta\right) + 2 \cdot 1}}\\
\mathbf{else}:\\
\;\;\;\;\frac{0}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}\\
\end{array}double f(double alpha, double beta) {
double r403728 = alpha;
double r403729 = beta;
double r403730 = r403728 + r403729;
double r403731 = r403729 * r403728;
double r403732 = r403730 + r403731;
double r403733 = 1.0;
double r403734 = r403732 + r403733;
double r403735 = 2.0;
double r403736 = r403735 * r403733;
double r403737 = r403730 + r403736;
double r403738 = r403734 / r403737;
double r403739 = r403738 / r403737;
double r403740 = r403737 + r403733;
double r403741 = r403739 / r403740;
return r403741;
}
double f(double alpha, double beta) {
double r403742 = beta;
double r403743 = 8.48955521481748e+157;
bool r403744 = r403742 <= r403743;
double r403745 = alpha;
double r403746 = r403745 + r403742;
double r403747 = r403742 * r403745;
double r403748 = r403746 + r403747;
double r403749 = 1.0;
double r403750 = r403748 + r403749;
double r403751 = sqrt(r403750);
double r403752 = sqrt(r403751);
double r403753 = 2.0;
double r403754 = r403753 * r403749;
double r403755 = r403746 + r403754;
double r403756 = r403755 + r403749;
double r403757 = r403750 / r403755;
double r403758 = sqrt(r403757);
double r403759 = r403756 / r403758;
double r403760 = r403752 / r403759;
double r403761 = r403751 / r403755;
double r403762 = r403761 / r403755;
double r403763 = sqrt(r403762);
double r403764 = sqrt(r403755);
double r403765 = r403763 / r403764;
double r403766 = r403760 * r403765;
double r403767 = 0.0;
double r403768 = r403767 / r403756;
double r403769 = r403744 ? r403766 : r403768;
return r403769;
}



Bits error versus alpha



Bits error versus beta
Results
if beta < 8.48955521481748e+157Initial program 1.3
rmApplied add-sqr-sqrt1.4
Applied associate-/l*1.4
rmApplied sqrt-div1.4
Applied associate-/r/1.3
Applied *-un-lft-identity1.3
Applied *-un-lft-identity1.3
Applied add-sqr-sqrt1.3
Applied times-frac1.3
Applied times-frac1.3
Applied sqrt-prod1.4
Applied times-frac1.4
if 8.48955521481748e+157 < beta Initial program 15.2
Taylor expanded around inf 7.5
Final simplification2.4
herbie shell --seed 2020042
(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)))