\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}\;\alpha \le 1.794632074747508277370817417090895285519 \cdot 10^{178}:\\
\;\;\;\;\frac{\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 \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}}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}double f(double alpha, double beta, double i) {
double r206335 = i;
double r206336 = alpha;
double r206337 = beta;
double r206338 = r206336 + r206337;
double r206339 = r206338 + r206335;
double r206340 = r206335 * r206339;
double r206341 = r206337 * r206336;
double r206342 = r206341 + r206340;
double r206343 = r206340 * r206342;
double r206344 = 2.0;
double r206345 = r206344 * r206335;
double r206346 = r206338 + r206345;
double r206347 = r206346 * r206346;
double r206348 = r206343 / r206347;
double r206349 = 1.0;
double r206350 = r206347 - r206349;
double r206351 = r206348 / r206350;
return r206351;
}
double f(double alpha, double beta, double i) {
double r206352 = alpha;
double r206353 = 1.7946320747475083e+178;
bool r206354 = r206352 <= r206353;
double r206355 = i;
double r206356 = beta;
double r206357 = r206352 + r206356;
double r206358 = r206357 + r206355;
double r206359 = r206355 * r206358;
double r206360 = 2.0;
double r206361 = r206360 * r206355;
double r206362 = r206357 + r206361;
double r206363 = r206359 / r206362;
double r206364 = 1.0;
double r206365 = sqrt(r206364);
double r206366 = r206362 + r206365;
double r206367 = r206363 / r206366;
double r206368 = r206356 * r206352;
double r206369 = r206368 + r206359;
double r206370 = r206369 / r206362;
double r206371 = r206367 * r206370;
double r206372 = r206362 - r206365;
double r206373 = r206371 / r206372;
double r206374 = 0.0;
double r206375 = r206354 ? r206373 : r206374;
return r206375;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if alpha < 1.7946320747475083e+178Initial program 52.6
rmApplied add-sqr-sqrt52.6
Applied difference-of-squares52.6
Applied times-frac37.6
Applied times-frac35.6
rmApplied associate-*r/35.6
if 1.7946320747475083e+178 < alpha Initial program 64.0
rmApplied add-sqr-sqrt64.0
Applied difference-of-squares64.0
Applied times-frac56.9
Applied times-frac53.5
Taylor expanded around inf 46.9
Final simplification37.2
herbie shell --seed 2019362
(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)))