\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}\;\beta \le 4.07062621862513012 \cdot 10^{207}:\\
\;\;\;\;\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{\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\alpha + \beta\right) + 2 \cdot i}} \cdot \frac{\sqrt{\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 r88345 = i;
double r88346 = alpha;
double r88347 = beta;
double r88348 = r88346 + r88347;
double r88349 = r88348 + r88345;
double r88350 = r88345 * r88349;
double r88351 = r88347 * r88346;
double r88352 = r88351 + r88350;
double r88353 = r88350 * r88352;
double r88354 = 2.0;
double r88355 = r88354 * r88345;
double r88356 = r88348 + r88355;
double r88357 = r88356 * r88356;
double r88358 = r88353 / r88357;
double r88359 = 1.0;
double r88360 = r88357 - r88359;
double r88361 = r88358 / r88360;
return r88361;
}
double f(double alpha, double beta, double i) {
double r88362 = beta;
double r88363 = 4.07062621862513e+207;
bool r88364 = r88362 <= r88363;
double r88365 = i;
double r88366 = alpha;
double r88367 = r88366 + r88362;
double r88368 = r88367 + r88365;
double r88369 = r88365 * r88368;
double r88370 = 2.0;
double r88371 = r88370 * r88365;
double r88372 = r88367 + r88371;
double r88373 = r88369 / r88372;
double r88374 = 1.0;
double r88375 = sqrt(r88374);
double r88376 = r88372 + r88375;
double r88377 = r88373 / r88376;
double r88378 = r88362 * r88366;
double r88379 = r88378 + r88369;
double r88380 = r88379 / r88372;
double r88381 = sqrt(r88380);
double r88382 = r88372 - r88375;
double r88383 = r88381 / r88382;
double r88384 = r88381 * r88383;
double r88385 = r88377 * r88384;
double r88386 = 0.0;
double r88387 = r88364 ? r88385 : r88386;
return r88387;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if beta < 4.07062621862513e+207Initial program 52.8
rmApplied add-sqr-sqrt52.8
Applied difference-of-squares52.8
Applied times-frac37.6
Applied times-frac35.3
rmApplied *-un-lft-identity35.3
Applied add-sqr-sqrt35.4
Applied times-frac35.4
Simplified35.4
if 4.07062621862513e+207 < beta Initial program 64.0
Taylor expanded around inf 44.1
Final simplification36.4
herbie shell --seed 2019198
(FPCore (alpha beta i)
:name "Octave 3.8, jcobi/4"
:pre (and (> alpha -1.0) (> beta -1.0) (> i 1.0))
(/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i)))) (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0)))