\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 1.8209498526474936 \cdot 10^{+220}:\\
\;\;\;\;\frac{\frac{i \cdot \left(i + \left(\beta + \alpha\right)\right) + \alpha \cdot \beta}{\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)}}{\left(\left(\beta + \alpha\right) + 2 \cdot i\right) - \sqrt{1.0}}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}double f(double alpha, double beta, double i) {
double r31607706 = i;
double r31607707 = alpha;
double r31607708 = beta;
double r31607709 = r31607707 + r31607708;
double r31607710 = r31607709 + r31607706;
double r31607711 = r31607706 * r31607710;
double r31607712 = r31607708 * r31607707;
double r31607713 = r31607712 + r31607711;
double r31607714 = r31607711 * r31607713;
double r31607715 = 2.0;
double r31607716 = r31607715 * r31607706;
double r31607717 = r31607709 + r31607716;
double r31607718 = r31607717 * r31607717;
double r31607719 = r31607714 / r31607718;
double r31607720 = 1.0;
double r31607721 = r31607718 - r31607720;
double r31607722 = r31607719 / r31607721;
return r31607722;
}
double f(double alpha, double beta, double i) {
double r31607723 = beta;
double r31607724 = 1.8209498526474936e+220;
bool r31607725 = r31607723 <= r31607724;
double r31607726 = i;
double r31607727 = alpha;
double r31607728 = r31607723 + r31607727;
double r31607729 = r31607726 + r31607728;
double r31607730 = r31607726 * r31607729;
double r31607731 = r31607727 * r31607723;
double r31607732 = r31607730 + r31607731;
double r31607733 = 2.0;
double r31607734 = r31607733 * r31607726;
double r31607735 = r31607728 + r31607734;
double r31607736 = r31607732 / r31607735;
double r31607737 = r31607730 / r31607735;
double r31607738 = 1.0;
double r31607739 = sqrt(r31607738);
double r31607740 = r31607739 + r31607735;
double r31607741 = r31607737 / r31607740;
double r31607742 = r31607736 * r31607741;
double r31607743 = r31607735 - r31607739;
double r31607744 = r31607742 / r31607743;
double r31607745 = 0.0;
double r31607746 = r31607725 ? r31607744 : r31607745;
return r31607746;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if beta < 1.8209498526474936e+220Initial program 51.6
rmApplied add-sqr-sqrt51.6
Applied difference-of-squares51.6
Applied times-frac36.8
Applied times-frac34.6
rmApplied associate-*r/34.5
if 1.8209498526474936e+220 < beta Initial program 62.6
rmApplied add-sqr-sqrt62.6
Applied difference-of-squares62.6
Applied times-frac56.6
Applied times-frac56.1
Taylor expanded around -inf 43.0
Final simplification35.4
herbie shell --seed 2019104
(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)))