\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}\;i \le 3.7225611126883413 \cdot 10^{87}:\\
\;\;\;\;\frac{\frac{\frac{i}{\frac{\left(\alpha + \beta\right) + 2 \cdot i}{\left(\alpha + \beta\right) + i}}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}} \cdot \frac{\sqrt{\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)}}{\frac{\left(\alpha + \beta\right) + 2 \cdot i}{\sqrt{\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)}}}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - \sqrt{1}}\\
\mathbf{elif}\;i \le 6.72619285808465442 \cdot 10^{153}:\\
\;\;\;\;\frac{0.25 \cdot {i}^{2}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}double f(double alpha, double beta, double i) {
double r133304 = i;
double r133305 = alpha;
double r133306 = beta;
double r133307 = r133305 + r133306;
double r133308 = r133307 + r133304;
double r133309 = r133304 * r133308;
double r133310 = r133306 * r133305;
double r133311 = r133310 + r133309;
double r133312 = r133309 * r133311;
double r133313 = 2.0;
double r133314 = r133313 * r133304;
double r133315 = r133307 + r133314;
double r133316 = r133315 * r133315;
double r133317 = r133312 / r133316;
double r133318 = 1.0;
double r133319 = r133316 - r133318;
double r133320 = r133317 / r133319;
return r133320;
}
double f(double alpha, double beta, double i) {
double r133321 = i;
double r133322 = 3.7225611126883413e+87;
bool r133323 = r133321 <= r133322;
double r133324 = alpha;
double r133325 = beta;
double r133326 = r133324 + r133325;
double r133327 = 2.0;
double r133328 = r133327 * r133321;
double r133329 = r133326 + r133328;
double r133330 = r133326 + r133321;
double r133331 = r133329 / r133330;
double r133332 = r133321 / r133331;
double r133333 = 1.0;
double r133334 = sqrt(r133333);
double r133335 = r133329 + r133334;
double r133336 = r133332 / r133335;
double r133337 = r133325 * r133324;
double r133338 = r133321 * r133330;
double r133339 = r133337 + r133338;
double r133340 = sqrt(r133339);
double r133341 = r133329 / r133340;
double r133342 = r133340 / r133341;
double r133343 = r133336 * r133342;
double r133344 = r133329 - r133334;
double r133345 = r133343 / r133344;
double r133346 = 6.726192858084654e+153;
bool r133347 = r133321 <= r133346;
double r133348 = 0.25;
double r133349 = 2.0;
double r133350 = pow(r133321, r133349);
double r133351 = r133348 * r133350;
double r133352 = r133329 * r133329;
double r133353 = r133352 - r133333;
double r133354 = r133351 / r133353;
double r133355 = 0.0;
double r133356 = r133347 ? r133354 : r133355;
double r133357 = r133323 ? r133345 : r133356;
return r133357;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if i < 3.7225611126883413e+87Initial program 29.1
rmApplied add-sqr-sqrt29.1
Applied difference-of-squares29.1
Applied times-frac12.0
Applied times-frac7.6
rmApplied associate-*r/7.6
rmApplied add-sqr-sqrt7.6
Applied associate-/l*7.6
rmApplied associate-/l*7.6
if 3.7225611126883413e+87 < i < 6.726192858084654e+153Initial program 64.0
Taylor expanded around inf 18.4
if 6.726192858084654e+153 < i Initial program 64.0
Taylor expanded around inf 61.9
Final simplification37.2
herbie shell --seed 2020045
(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)))