\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 6.529674909757050223649546242132507716114 \cdot 10^{195}:\\
\;\;\;\;\frac{\frac{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \beta \cdot \alpha}{\left(\alpha + \beta\right) + 2 \cdot i} \cdot \frac{i \cdot \left(i + \left(\alpha + \beta\right)\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\sqrt{1} + \left(\left(\alpha + \beta\right) + 2 \cdot i\right)} \cdot \frac{1}{\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 r4453281 = i;
double r4453282 = alpha;
double r4453283 = beta;
double r4453284 = r4453282 + r4453283;
double r4453285 = r4453284 + r4453281;
double r4453286 = r4453281 * r4453285;
double r4453287 = r4453283 * r4453282;
double r4453288 = r4453287 + r4453286;
double r4453289 = r4453286 * r4453288;
double r4453290 = 2.0;
double r4453291 = r4453290 * r4453281;
double r4453292 = r4453284 + r4453291;
double r4453293 = r4453292 * r4453292;
double r4453294 = r4453289 / r4453293;
double r4453295 = 1.0;
double r4453296 = r4453293 - r4453295;
double r4453297 = r4453294 / r4453296;
return r4453297;
}
double f(double alpha, double beta, double i) {
double r4453298 = alpha;
double r4453299 = 6.52967490975705e+195;
bool r4453300 = r4453298 <= r4453299;
double r4453301 = i;
double r4453302 = beta;
double r4453303 = r4453298 + r4453302;
double r4453304 = r4453301 + r4453303;
double r4453305 = r4453301 * r4453304;
double r4453306 = r4453302 * r4453298;
double r4453307 = r4453305 + r4453306;
double r4453308 = 2.0;
double r4453309 = r4453308 * r4453301;
double r4453310 = r4453303 + r4453309;
double r4453311 = r4453307 / r4453310;
double r4453312 = r4453305 / r4453310;
double r4453313 = r4453311 * r4453312;
double r4453314 = 1.0;
double r4453315 = sqrt(r4453314);
double r4453316 = r4453315 + r4453310;
double r4453317 = r4453313 / r4453316;
double r4453318 = 1.0;
double r4453319 = r4453310 - r4453315;
double r4453320 = r4453318 / r4453319;
double r4453321 = r4453317 * r4453320;
double r4453322 = 0.0;
double r4453323 = r4453300 ? r4453321 : r4453322;
return r4453323;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if alpha < 6.52967490975705e+195Initial program 52.5
rmApplied add-sqr-sqrt52.5
Applied difference-of-squares52.5
Applied times-frac37.2
Applied times-frac35.0
rmApplied div-inv35.0
Applied associate-*r*35.0
rmApplied associate-*l/35.0
if 6.52967490975705e+195 < alpha Initial program 64.0
Taylor expanded around inf 44.1
Final simplification36.1
herbie shell --seed 2019170
(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)))