\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}\;i \le 1.337632796530924991507977374043559908369 \cdot 10^{154}:\\
\;\;\;\;\frac{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\alpha + \beta\right) + 2 \cdot i} \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}}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}double f(double alpha, double beta, double i) {
double r200645 = i;
double r200646 = alpha;
double r200647 = beta;
double r200648 = r200646 + r200647;
double r200649 = r200648 + r200645;
double r200650 = r200645 * r200649;
double r200651 = r200647 * r200646;
double r200652 = r200651 + r200650;
double r200653 = r200650 * r200652;
double r200654 = 2.0;
double r200655 = r200654 * r200645;
double r200656 = r200648 + r200655;
double r200657 = r200656 * r200656;
double r200658 = r200653 / r200657;
double r200659 = 1.0;
double r200660 = r200657 - r200659;
double r200661 = r200658 / r200660;
return r200661;
}
double f(double alpha, double beta, double i) {
double r200662 = i;
double r200663 = 1.337632796530925e+154;
bool r200664 = r200662 <= r200663;
double r200665 = alpha;
double r200666 = beta;
double r200667 = r200665 + r200666;
double r200668 = r200667 + r200662;
double r200669 = r200662 * r200668;
double r200670 = 2.0;
double r200671 = r200670 * r200662;
double r200672 = r200667 + r200671;
double r200673 = r200669 / r200672;
double r200674 = r200666 * r200665;
double r200675 = r200674 + r200669;
double r200676 = r200675 / r200672;
double r200677 = 1.0;
double r200678 = sqrt(r200677);
double r200679 = r200672 - r200678;
double r200680 = r200676 / r200679;
double r200681 = r200673 * r200680;
double r200682 = r200672 + r200678;
double r200683 = r200681 / r200682;
double r200684 = 0.0;
double r200685 = r200664 ? r200683 : r200684;
return r200685;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if i < 1.337632796530925e+154Initial program 44.8
rmApplied add-sqr-sqrt44.8
Applied difference-of-squares44.8
Applied times-frac15.9
Applied times-frac11.7
rmApplied associate-*l/11.6
if 1.337632796530925e+154 < i Initial program 64.0
Taylor expanded around inf 61.9
Final simplification36.7
herbie shell --seed 2019352
(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)))