\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(\beta + 2 \cdot i\right) + \alpha}{\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 r189316 = i;
double r189317 = alpha;
double r189318 = beta;
double r189319 = r189317 + r189318;
double r189320 = r189319 + r189316;
double r189321 = r189316 * r189320;
double r189322 = r189318 * r189317;
double r189323 = r189322 + r189321;
double r189324 = r189321 * r189323;
double r189325 = 2.0;
double r189326 = r189325 * r189316;
double r189327 = r189319 + r189326;
double r189328 = r189327 * r189327;
double r189329 = r189324 / r189328;
double r189330 = 1.0;
double r189331 = r189328 - r189330;
double r189332 = r189329 / r189331;
return r189332;
}
double f(double alpha, double beta, double i) {
double r189333 = i;
double r189334 = 3.7225611126883413e+87;
bool r189335 = r189333 <= r189334;
double r189336 = beta;
double r189337 = 2.0;
double r189338 = r189337 * r189333;
double r189339 = r189336 + r189338;
double r189340 = alpha;
double r189341 = r189339 + r189340;
double r189342 = r189340 + r189336;
double r189343 = r189342 + r189333;
double r189344 = r189341 / r189343;
double r189345 = r189333 / r189344;
double r189346 = r189342 + r189338;
double r189347 = 1.0;
double r189348 = sqrt(r189347);
double r189349 = r189346 + r189348;
double r189350 = r189345 / r189349;
double r189351 = r189336 * r189340;
double r189352 = r189333 * r189343;
double r189353 = r189351 + r189352;
double r189354 = sqrt(r189353);
double r189355 = r189346 / r189354;
double r189356 = r189354 / r189355;
double r189357 = r189350 * r189356;
double r189358 = r189346 - r189348;
double r189359 = r189357 / r189358;
double r189360 = 6.726192858084654e+153;
bool r189361 = r189333 <= r189360;
double r189362 = 0.25;
double r189363 = 2.0;
double r189364 = pow(r189333, r189363);
double r189365 = r189362 * r189364;
double r189366 = r189346 * r189346;
double r189367 = r189366 - r189347;
double r189368 = r189365 / r189367;
double r189369 = 0.0;
double r189370 = r189361 ? r189368 : r189369;
double r189371 = r189335 ? r189359 : r189370;
return r189371;
}



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
Simplified7.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)))