\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.0}\begin{array}{l}
\mathbf{if}\;\beta \le 2.8543550670306893 \cdot 10^{+154}:\\
\;\;\;\;\sqrt{\frac{1}{\frac{\left(\left(\beta + \alpha\right) + 2 \cdot i\right) - \sqrt{1.0}}{\frac{\alpha \cdot \beta + i \cdot \left(i + \left(\beta + \alpha\right)\right)}{\left(\beta + \alpha\right) + 2 \cdot i}}} \cdot \frac{\frac{i \cdot \left(i + \left(\beta + \alpha\right)\right)}{\left(\beta + \alpha\right) + 2 \cdot i}}{\sqrt{1.0} + \left(\left(\beta + \alpha\right) + 2 \cdot i\right)}} \cdot \sqrt{\frac{1}{\frac{\left(\left(\beta + \alpha\right) + 2 \cdot i\right) - \sqrt{1.0}}{\frac{\alpha \cdot \beta + i \cdot \left(i + \left(\beta + \alpha\right)\right)}{\left(\beta + \alpha\right) + 2 \cdot i}}} \cdot \frac{\frac{i \cdot \left(i + \left(\beta + \alpha\right)\right)}{\left(\beta + \alpha\right) + 2 \cdot i}}{\sqrt{1.0} + \left(\left(\beta + \alpha\right) + 2 \cdot i\right)}}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}double f(double alpha, double beta, double i) {
double r4841748 = i;
double r4841749 = alpha;
double r4841750 = beta;
double r4841751 = r4841749 + r4841750;
double r4841752 = r4841751 + r4841748;
double r4841753 = r4841748 * r4841752;
double r4841754 = r4841750 * r4841749;
double r4841755 = r4841754 + r4841753;
double r4841756 = r4841753 * r4841755;
double r4841757 = 2.0;
double r4841758 = r4841757 * r4841748;
double r4841759 = r4841751 + r4841758;
double r4841760 = r4841759 * r4841759;
double r4841761 = r4841756 / r4841760;
double r4841762 = 1.0;
double r4841763 = r4841760 - r4841762;
double r4841764 = r4841761 / r4841763;
return r4841764;
}
double f(double alpha, double beta, double i) {
double r4841765 = beta;
double r4841766 = 2.8543550670306893e+154;
bool r4841767 = r4841765 <= r4841766;
double r4841768 = 1.0;
double r4841769 = alpha;
double r4841770 = r4841765 + r4841769;
double r4841771 = 2.0;
double r4841772 = i;
double r4841773 = r4841771 * r4841772;
double r4841774 = r4841770 + r4841773;
double r4841775 = 1.0;
double r4841776 = sqrt(r4841775);
double r4841777 = r4841774 - r4841776;
double r4841778 = r4841769 * r4841765;
double r4841779 = r4841772 + r4841770;
double r4841780 = r4841772 * r4841779;
double r4841781 = r4841778 + r4841780;
double r4841782 = r4841781 / r4841774;
double r4841783 = r4841777 / r4841782;
double r4841784 = r4841768 / r4841783;
double r4841785 = r4841780 / r4841774;
double r4841786 = r4841776 + r4841774;
double r4841787 = r4841785 / r4841786;
double r4841788 = r4841784 * r4841787;
double r4841789 = sqrt(r4841788);
double r4841790 = r4841789 * r4841789;
double r4841791 = 0.0;
double r4841792 = r4841767 ? r4841790 : r4841791;
return r4841792;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if beta < 2.8543550670306893e+154Initial program 50.3
rmApplied add-sqr-sqrt50.3
Applied difference-of-squares50.3
Applied times-frac35.5
Applied times-frac34.2
rmApplied *-un-lft-identity34.2
Applied associate-/l*34.2
rmApplied add-sqr-sqrt34.2
if 2.8543550670306893e+154 < beta Initial program 62.5
Taylor expanded around inf 47.2
Final simplification36.4
herbie shell --seed 2019138
(FPCore (alpha beta i)
:name "Octave 3.8, jcobi/4"
: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.0)))