\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 9.5093000760638009 \cdot 10^{51}:\\
\;\;\;\;\frac{e^{\log \left(\mathsf{fma}\left(\frac{\frac{\alpha + \beta}{1}}{1}, \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}, 1\right)\right)}}{2}\\
\mathbf{elif}\;\alpha \le 3.9744006379666212 \cdot 10^{102}:\\
\;\;\;\;\frac{\mathsf{fma}\left(-4, \frac{1}{{\alpha}^{2}}, \mathsf{fma}\left(8, \frac{1}{{\alpha}^{3}}, \frac{2}{\alpha}\right)\right)}{2}\\
\mathbf{elif}\;\alpha \le 4.4633300930535925 \cdot 10^{121}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\frac{\frac{\alpha + \beta}{1}}{1}, \frac{\frac{\frac{\beta - \alpha}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}}}{\mathsf{fma}\left(i, 2, \alpha + \beta\right)}}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}}, 1\right)}{2}\\
\mathbf{elif}\;\alpha \le 4.22071431264922 \cdot 10^{169}:\\
\;\;\;\;\frac{\mathsf{fma}\left(-4, \frac{1}{{\alpha}^{2}}, \mathsf{fma}\left(8, \frac{1}{{\alpha}^{3}}, \frac{2}{\alpha}\right)\right)}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(\frac{\frac{\beta}{\alpha + \left(\beta + \mathsf{fma}\left(i, 2, 2\right)\right)}}{\mathsf{fma}\left(i, 2, \alpha + \beta\right)} - \frac{\frac{\alpha}{\alpha + \left(\beta + \mathsf{fma}\left(i, 2, 2\right)\right)}}{\mathsf{fma}\left(i, 2, \alpha + \beta\right)}\right) \cdot \left(\alpha + \beta\right) + 1}{2}\\
\end{array}double code(double alpha, double beta, double i) {
return ((((((alpha + beta) * (beta - alpha)) / ((alpha + beta) + (2.0 * i))) / (((alpha + beta) + (2.0 * i)) + 2.0)) + 1.0) / 2.0);
}
double code(double alpha, double beta, double i) {
double VAR;
if ((alpha <= 9.509300076063801e+51)) {
VAR = (exp(log(fma((((alpha + beta) / 1.0) / 1.0), (((beta - alpha) / ((alpha + beta) + (2.0 * i))) / (((alpha + beta) + (2.0 * i)) + 2.0)), 1.0))) / 2.0);
} else {
double VAR_1;
if ((alpha <= 3.974400637966621e+102)) {
VAR_1 = (fma(-4.0, (1.0 / pow(alpha, 2.0)), fma(8.0, (1.0 / pow(alpha, 3.0)), (2.0 / alpha))) / 2.0);
} else {
double VAR_2;
if ((alpha <= 4.4633300930535925e+121)) {
VAR_2 = (fma((((alpha + beta) / 1.0) / 1.0), ((((beta - alpha) / sqrt((((alpha + beta) + (2.0 * i)) + 2.0))) / fma(i, 2.0, (alpha + beta))) / sqrt((((alpha + beta) + (2.0 * i)) + 2.0))), 1.0) / 2.0);
} else {
double VAR_3;
if ((alpha <= 4.22071431264922e+169)) {
VAR_3 = (fma(-4.0, (1.0 / pow(alpha, 2.0)), fma(8.0, (1.0 / pow(alpha, 3.0)), (2.0 / alpha))) / 2.0);
} else {
VAR_3 = ((((((beta / (alpha + (beta + fma(i, 2.0, 2.0)))) / fma(i, 2.0, (alpha + beta))) - ((alpha / (alpha + (beta + fma(i, 2.0, 2.0)))) / fma(i, 2.0, (alpha + beta)))) * (alpha + beta)) + 1.0) / 2.0);
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if alpha < 9.509300076063801e+51Initial program 12.0
rmApplied *-un-lft-identity12.0
Applied *-un-lft-identity12.0
Applied times-frac1.2
Applied times-frac1.2
Applied fma-def1.2
rmApplied add-exp-log1.2
if 9.509300076063801e+51 < alpha < 3.974400637966621e+102 or 4.4633300930535925e+121 < alpha < 4.22071431264922e+169Initial program 44.4
rmApplied *-un-lft-identity44.4
Applied *-un-lft-identity44.4
Applied times-frac31.1
Applied times-frac30.9
Applied fma-def31.0
Taylor expanded around inf 40.4
Simplified40.4
if 3.974400637966621e+102 < alpha < 4.4633300930535925e+121Initial program 43.4
rmApplied *-un-lft-identity43.4
Applied *-un-lft-identity43.4
Applied times-frac33.0
Applied times-frac32.8
Applied fma-def32.7
rmApplied add-sqr-sqrt32.8
Applied associate-/r*32.8
Simplified32.8
if 4.22071431264922e+169 < alpha Initial program 64.0
rmApplied *-un-lft-identity64.0
Applied *-un-lft-identity64.0
Applied times-frac48.1
Applied times-frac48.2
Applied fma-def49.1
rmApplied div-sub49.1
Applied div-sub49.1
Simplified49.1
Simplified49.1
rmApplied fma-udef48.2
Simplified48.2
Final simplification13.2
herbie shell --seed 2020102 +o rules:numerics
(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))