\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}\;\alpha \le 9.866752982876504 \cdot 10^{+203}:\\
\;\;\;\;\frac{\frac{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \beta \cdot \alpha}{\left(\alpha + \beta\right) + 2 \cdot i} \cdot \frac{\frac{i \cdot \left(i + \left(\alpha + \beta\right)\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\sqrt{1.0} + \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - \sqrt{1.0}}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}double f(double alpha, double beta, double i) {
double r5339108 = i;
double r5339109 = alpha;
double r5339110 = beta;
double r5339111 = r5339109 + r5339110;
double r5339112 = r5339111 + r5339108;
double r5339113 = r5339108 * r5339112;
double r5339114 = r5339110 * r5339109;
double r5339115 = r5339114 + r5339113;
double r5339116 = r5339113 * r5339115;
double r5339117 = 2.0;
double r5339118 = r5339117 * r5339108;
double r5339119 = r5339111 + r5339118;
double r5339120 = r5339119 * r5339119;
double r5339121 = r5339116 / r5339120;
double r5339122 = 1.0;
double r5339123 = r5339120 - r5339122;
double r5339124 = r5339121 / r5339123;
return r5339124;
}
double f(double alpha, double beta, double i) {
double r5339125 = alpha;
double r5339126 = 9.866752982876504e+203;
bool r5339127 = r5339125 <= r5339126;
double r5339128 = i;
double r5339129 = beta;
double r5339130 = r5339125 + r5339129;
double r5339131 = r5339128 + r5339130;
double r5339132 = r5339128 * r5339131;
double r5339133 = r5339129 * r5339125;
double r5339134 = r5339132 + r5339133;
double r5339135 = 2.0;
double r5339136 = r5339135 * r5339128;
double r5339137 = r5339130 + r5339136;
double r5339138 = r5339134 / r5339137;
double r5339139 = r5339132 / r5339137;
double r5339140 = 1.0;
double r5339141 = sqrt(r5339140);
double r5339142 = r5339141 + r5339137;
double r5339143 = r5339139 / r5339142;
double r5339144 = r5339138 * r5339143;
double r5339145 = r5339137 - r5339141;
double r5339146 = r5339144 / r5339145;
double r5339147 = 0.0;
double r5339148 = r5339127 ? r5339146 : r5339147;
return r5339148;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if alpha < 9.866752982876504e+203Initial program 51.0
rmApplied add-sqr-sqrt51.0
Applied difference-of-squares51.0
Applied times-frac36.3
Applied times-frac34.3
rmApplied associate-*r/34.3
if 9.866752982876504e+203 < alpha Initial program 62.6
Taylor expanded around inf 45.5
Final simplification35.5
herbie shell --seed 2019164
(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)))