\frac{\frac{\frac{\left(\alpha + \beta\right) \cdot \left(\beta - \alpha\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + 1}{2}\begin{array}{l}
\mathbf{if}\;\alpha \le 19393672272838702021867389529627543207940:\\
\;\;\;\;\frac{\frac{\log \left(e^{{\left(\frac{\frac{1}{\sqrt{1}}}{\frac{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}}{\frac{\alpha + \beta}{\sqrt{1}}}}\right)}^{3} + {1}^{3}}\right)}{\frac{\frac{1}{\sqrt{1}}}{\frac{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}}{\frac{\alpha + \beta}{\sqrt{1}}}} \cdot \frac{\frac{1}{\sqrt{1}}}{\frac{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}}{\frac{\alpha + \beta}{\sqrt{1}}}} + \left(1 \cdot 1 - \frac{\frac{1}{\sqrt{1}}}{\frac{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}}{\frac{\alpha + \beta}{\sqrt{1}}}} \cdot 1\right)}}{2}\\
\mathbf{elif}\;\alpha \le 1.262557062960945982227949818049794246774 \cdot 10^{165}:\\
\;\;\;\;\frac{\left(2 \cdot \frac{1}{\alpha} + 8 \cdot \frac{1}{{\alpha}^{3}}\right) - 4 \cdot \frac{1}{{\alpha}^{2}}}{2}\\
\mathbf{elif}\;\alpha \le 1.468201707068765789213963828469773650324 \cdot 10^{209}:\\
\;\;\;\;\frac{\frac{\frac{\frac{1}{\sqrt{1}}}{\frac{\frac{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}}{1}}{\alpha + \beta}}}{\frac{\frac{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}}{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}}{\frac{1}{\sqrt{1}}}} + 1}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(2 \cdot \frac{1}{\alpha} + 8 \cdot \frac{1}{{\alpha}^{3}}\right) - 4 \cdot \frac{1}{{\alpha}^{2}}}{2}\\
\end{array}double f(double alpha, double beta, double i) {
double r794286 = alpha;
double r794287 = beta;
double r794288 = r794286 + r794287;
double r794289 = r794287 - r794286;
double r794290 = r794288 * r794289;
double r794291 = 2.0;
double r794292 = i;
double r794293 = r794291 * r794292;
double r794294 = r794288 + r794293;
double r794295 = r794290 / r794294;
double r794296 = r794294 + r794291;
double r794297 = r794295 / r794296;
double r794298 = 1.0;
double r794299 = r794297 + r794298;
double r794300 = r794299 / r794291;
return r794300;
}
double f(double alpha, double beta, double i) {
double r794301 = alpha;
double r794302 = 1.9393672272838702e+40;
bool r794303 = r794301 <= r794302;
double r794304 = 1.0;
double r794305 = sqrt(r794304);
double r794306 = r794304 / r794305;
double r794307 = beta;
double r794308 = r794301 + r794307;
double r794309 = 2.0;
double r794310 = i;
double r794311 = r794309 * r794310;
double r794312 = r794308 + r794311;
double r794313 = r794312 + r794309;
double r794314 = r794307 - r794301;
double r794315 = r794314 / r794312;
double r794316 = r794313 / r794315;
double r794317 = r794308 / r794305;
double r794318 = r794316 / r794317;
double r794319 = r794306 / r794318;
double r794320 = 3.0;
double r794321 = pow(r794319, r794320);
double r794322 = 1.0;
double r794323 = pow(r794322, r794320);
double r794324 = r794321 + r794323;
double r794325 = exp(r794324);
double r794326 = log(r794325);
double r794327 = r794319 * r794319;
double r794328 = r794322 * r794322;
double r794329 = r794319 * r794322;
double r794330 = r794328 - r794329;
double r794331 = r794327 + r794330;
double r794332 = r794326 / r794331;
double r794333 = r794332 / r794309;
double r794334 = 1.262557062960946e+165;
bool r794335 = r794301 <= r794334;
double r794336 = r794304 / r794301;
double r794337 = r794309 * r794336;
double r794338 = 8.0;
double r794339 = pow(r794301, r794320);
double r794340 = r794304 / r794339;
double r794341 = r794338 * r794340;
double r794342 = r794337 + r794341;
double r794343 = 4.0;
double r794344 = 2.0;
double r794345 = pow(r794301, r794344);
double r794346 = r794304 / r794345;
double r794347 = r794343 * r794346;
double r794348 = r794342 - r794347;
double r794349 = r794348 / r794309;
double r794350 = 1.4682017070687658e+209;
bool r794351 = r794301 <= r794350;
double r794352 = sqrt(r794313);
double r794353 = r794352 / r794304;
double r794354 = r794353 / r794308;
double r794355 = r794306 / r794354;
double r794356 = r794352 / r794315;
double r794357 = r794356 / r794306;
double r794358 = r794355 / r794357;
double r794359 = r794358 + r794322;
double r794360 = r794359 / r794309;
double r794361 = r794351 ? r794360 : r794349;
double r794362 = r794335 ? r794349 : r794361;
double r794363 = r794303 ? r794333 : r794362;
return r794363;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if alpha < 1.9393672272838702e+40Initial program 11.0
rmApplied *-un-lft-identity11.0
Applied times-frac1.1
Applied associate-/l*1.1
rmApplied add-sqr-sqrt1.1
Applied *-un-lft-identity1.1
Applied times-frac1.1
Applied associate-/l*1.1
rmApplied flip3-+1.1
rmApplied add-log-exp1.1
Applied add-log-exp1.1
Applied sum-log1.1
Simplified1.1
if 1.9393672272838702e+40 < alpha < 1.262557062960946e+165 or 1.4682017070687658e+209 < alpha Initial program 51.9
Taylor expanded around inf 41.5
if 1.262557062960946e+165 < alpha < 1.4682017070687658e+209Initial program 64.0
rmApplied *-un-lft-identity64.0
Applied times-frac41.8
Applied associate-/l*41.8
rmApplied add-sqr-sqrt41.8
Applied *-un-lft-identity41.8
Applied times-frac41.8
Applied associate-/l*41.8
rmApplied div-inv41.8
Applied *-un-lft-identity41.8
Applied *-un-lft-identity41.8
Applied times-frac41.8
Applied add-sqr-sqrt42.0
Applied times-frac42.0
Applied times-frac42.1
Applied associate-/r*42.0
Simplified42.0
Final simplification13.1
herbie shell --seed 2019298
(FPCore (alpha beta i)
:name "Octave 3.8, jcobi/2"
:precision binary64
:pre (and (> alpha -1) (> beta -1) (> i 0.0))
(/ (+ (/ (/ (* (+ alpha beta) (- beta alpha)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) 2)) 1) 2))