\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.0} + 1.0}{2.0}\begin{array}{l}
\mathbf{if}\;\alpha \le 2.0760474421781283 \cdot 10^{+39}:\\
\;\;\;\;\frac{1.0 + \frac{\frac{\left(\beta + \alpha\right) - 2 \cdot i}{\sqrt{\left(2 \cdot i + \left(\beta + \alpha\right)\right) + 2.0}}}{\sqrt{\left(2 \cdot i + \left(\beta + \alpha\right)\right) + 2.0}} \cdot \left(\frac{\beta + \alpha}{2 \cdot i + \left(\beta + \alpha\right)} \cdot \frac{\beta - \alpha}{\left(\beta + \alpha\right) - 2 \cdot i}\right)}{2.0}\\
\mathbf{elif}\;\alpha \le 3.6345013768301078 \cdot 10^{+118}:\\
\;\;\;\;\frac{\frac{2.0}{\alpha} + \left(\frac{8.0}{\left(\alpha \cdot \alpha\right) \cdot \alpha} - \frac{4.0}{\alpha \cdot \alpha}\right)}{2.0}\\
\mathbf{elif}\;\alpha \le 1.480588010249348 \cdot 10^{+179}:\\
\;\;\;\;\frac{\frac{{1.0}^{3} + {\left(\left(\left(\frac{\beta + \alpha}{2 \cdot i + \left(\beta + \alpha\right)} \cdot \left(\beta - \alpha\right)\right) \cdot \frac{1}{\left(\beta + \alpha\right) - 2 \cdot i}\right) \cdot \frac{\left(\beta + \alpha\right) - 2 \cdot i}{\left(2 \cdot i + \left(\beta + \alpha\right)\right) + 2.0}\right)}^{3}}{\left(\left(\left(\frac{\beta + \alpha}{2 \cdot i + \left(\beta + \alpha\right)} \cdot \left(\beta - \alpha\right)\right) \cdot \frac{1}{\left(\beta + \alpha\right) - 2 \cdot i}\right) \cdot \frac{\left(\beta + \alpha\right) - 2 \cdot i}{\left(2 \cdot i + \left(\beta + \alpha\right)\right) + 2.0}\right) \cdot \left(\left(\left(\frac{\beta + \alpha}{2 \cdot i + \left(\beta + \alpha\right)} \cdot \left(\beta - \alpha\right)\right) \cdot \frac{1}{\left(\beta + \alpha\right) - 2 \cdot i}\right) \cdot \frac{\left(\beta + \alpha\right) - 2 \cdot i}{\left(2 \cdot i + \left(\beta + \alpha\right)\right) + 2.0}\right) + \left(1.0 \cdot 1.0 - 1.0 \cdot \left(\left(\left(\frac{\beta + \alpha}{2 \cdot i + \left(\beta + \alpha\right)} \cdot \left(\beta - \alpha\right)\right) \cdot \frac{1}{\left(\beta + \alpha\right) - 2 \cdot i}\right) \cdot \frac{\left(\beta + \alpha\right) - 2 \cdot i}{\left(2 \cdot i + \left(\beta + \alpha\right)\right) + 2.0}\right)\right)}}{2.0}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2.0}{\alpha} + \left(\frac{8.0}{\left(\alpha \cdot \alpha\right) \cdot \alpha} - \frac{4.0}{\alpha \cdot \alpha}\right)}{2.0}\\
\end{array}double f(double alpha, double beta, double i) {
double r2055739 = alpha;
double r2055740 = beta;
double r2055741 = r2055739 + r2055740;
double r2055742 = r2055740 - r2055739;
double r2055743 = r2055741 * r2055742;
double r2055744 = 2.0;
double r2055745 = i;
double r2055746 = r2055744 * r2055745;
double r2055747 = r2055741 + r2055746;
double r2055748 = r2055743 / r2055747;
double r2055749 = 2.0;
double r2055750 = r2055747 + r2055749;
double r2055751 = r2055748 / r2055750;
double r2055752 = 1.0;
double r2055753 = r2055751 + r2055752;
double r2055754 = r2055753 / r2055749;
return r2055754;
}
double f(double alpha, double beta, double i) {
double r2055755 = alpha;
double r2055756 = 2.0760474421781283e+39;
bool r2055757 = r2055755 <= r2055756;
double r2055758 = 1.0;
double r2055759 = beta;
double r2055760 = r2055759 + r2055755;
double r2055761 = 2.0;
double r2055762 = i;
double r2055763 = r2055761 * r2055762;
double r2055764 = r2055760 - r2055763;
double r2055765 = r2055763 + r2055760;
double r2055766 = 2.0;
double r2055767 = r2055765 + r2055766;
double r2055768 = sqrt(r2055767);
double r2055769 = r2055764 / r2055768;
double r2055770 = r2055769 / r2055768;
double r2055771 = r2055760 / r2055765;
double r2055772 = r2055759 - r2055755;
double r2055773 = r2055772 / r2055764;
double r2055774 = r2055771 * r2055773;
double r2055775 = r2055770 * r2055774;
double r2055776 = r2055758 + r2055775;
double r2055777 = r2055776 / r2055766;
double r2055778 = 3.6345013768301078e+118;
bool r2055779 = r2055755 <= r2055778;
double r2055780 = r2055766 / r2055755;
double r2055781 = 8.0;
double r2055782 = r2055755 * r2055755;
double r2055783 = r2055782 * r2055755;
double r2055784 = r2055781 / r2055783;
double r2055785 = 4.0;
double r2055786 = r2055785 / r2055782;
double r2055787 = r2055784 - r2055786;
double r2055788 = r2055780 + r2055787;
double r2055789 = r2055788 / r2055766;
double r2055790 = 1.480588010249348e+179;
bool r2055791 = r2055755 <= r2055790;
double r2055792 = 3.0;
double r2055793 = pow(r2055758, r2055792);
double r2055794 = r2055771 * r2055772;
double r2055795 = 1.0;
double r2055796 = r2055795 / r2055764;
double r2055797 = r2055794 * r2055796;
double r2055798 = r2055764 / r2055767;
double r2055799 = r2055797 * r2055798;
double r2055800 = pow(r2055799, r2055792);
double r2055801 = r2055793 + r2055800;
double r2055802 = r2055799 * r2055799;
double r2055803 = r2055758 * r2055758;
double r2055804 = r2055758 * r2055799;
double r2055805 = r2055803 - r2055804;
double r2055806 = r2055802 + r2055805;
double r2055807 = r2055801 / r2055806;
double r2055808 = r2055807 / r2055766;
double r2055809 = r2055791 ? r2055808 : r2055789;
double r2055810 = r2055779 ? r2055789 : r2055809;
double r2055811 = r2055757 ? r2055777 : r2055810;
return r2055811;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if alpha < 2.0760474421781283e+39Initial program 11.3
rmApplied *-un-lft-identity11.3
Applied flip-+13.5
Applied associate-/r/13.5
Applied times-frac13.5
Simplified1.0
rmApplied add-sqr-sqrt1.1
Applied associate-/r*1.1
if 2.0760474421781283e+39 < alpha < 3.6345013768301078e+118 or 1.480588010249348e+179 < alpha Initial program 54.5
rmApplied *-un-lft-identity54.5
Applied flip-+54.1
Applied associate-/r/54.1
Applied times-frac54.1
Simplified41.7
Taylor expanded around inf 40.8
Simplified40.8
if 3.6345013768301078e+118 < alpha < 1.480588010249348e+179Initial program 51.7
rmApplied *-un-lft-identity51.7
Applied flip-+51.4
Applied associate-/r/51.4
Applied times-frac51.4
Simplified37.4
rmApplied div-inv37.3
Applied associate-*r*37.3
rmApplied flip3-+37.3
Final simplification12.5
herbie shell --seed 2019155
(FPCore (alpha beta i)
:name "Octave 3.8, jcobi/2"
:pre (and (> alpha -1) (> beta -1) (> i 0))
(/ (+ (/ (/ (* (+ alpha beta) (- beta alpha)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) 2.0)) 1.0) 2.0))