\frac{\frac{\left(i \cdot \left(\left(\alpha + \beta\right) + i\right)\right) \cdot \left(\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1}\begin{array}{l}
\mathbf{if}\;\alpha \le 9.46141053092278866 \cdot 10^{209}:\\
\;\;\;\;\sqrt{\frac{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}}} \cdot \left(\sqrt{\frac{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}}} \cdot \frac{\frac{\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - \sqrt{1}}\right)\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}double f(double alpha, double beta, double i) {
double r188565 = i;
double r188566 = alpha;
double r188567 = beta;
double r188568 = r188566 + r188567;
double r188569 = r188568 + r188565;
double r188570 = r188565 * r188569;
double r188571 = r188567 * r188566;
double r188572 = r188571 + r188570;
double r188573 = r188570 * r188572;
double r188574 = 2.0;
double r188575 = r188574 * r188565;
double r188576 = r188568 + r188575;
double r188577 = r188576 * r188576;
double r188578 = r188573 / r188577;
double r188579 = 1.0;
double r188580 = r188577 - r188579;
double r188581 = r188578 / r188580;
return r188581;
}
double f(double alpha, double beta, double i) {
double r188582 = alpha;
double r188583 = 9.461410530922789e+209;
bool r188584 = r188582 <= r188583;
double r188585 = i;
double r188586 = beta;
double r188587 = r188582 + r188586;
double r188588 = r188587 + r188585;
double r188589 = r188585 * r188588;
double r188590 = 2.0;
double r188591 = r188590 * r188585;
double r188592 = r188587 + r188591;
double r188593 = r188589 / r188592;
double r188594 = 1.0;
double r188595 = sqrt(r188594);
double r188596 = r188592 + r188595;
double r188597 = r188593 / r188596;
double r188598 = sqrt(r188597);
double r188599 = r188586 * r188582;
double r188600 = r188599 + r188589;
double r188601 = r188600 / r188592;
double r188602 = r188592 - r188595;
double r188603 = r188601 / r188602;
double r188604 = r188598 * r188603;
double r188605 = r188598 * r188604;
double r188606 = 0.0;
double r188607 = r188584 ? r188605 : r188606;
return r188607;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if alpha < 9.461410530922789e+209Initial program 52.5
rmApplied add-sqr-sqrt52.5
Applied difference-of-squares52.5
Applied times-frac37.5
Applied times-frac35.1
rmApplied add-sqr-sqrt35.1
Applied associate-*l*35.1
if 9.461410530922789e+209 < alpha Initial program 64.0
rmApplied add-sqr-sqrt64.0
Applied difference-of-squares64.0
Applied times-frac56.3
Applied times-frac55.3
Taylor expanded around inf 42.2
Final simplification35.9
herbie shell --seed 2020003
(FPCore (alpha beta i)
:name "Octave 3.8, jcobi/4"
:precision binary64
:pre (and (> alpha -1) (> beta -1) (> i 1))
(/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1)))