\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 6.1912244088489427 \cdot 10^{177}:\\
\;\;\;\;\frac{\left(\alpha + \beta\right) \cdot \frac{\frac{\beta}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + \sqrt[3]{{\left(\left(-\frac{\frac{\alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}\right) \cdot \left(\alpha + \beta\right) + 1\right)}^{3}}}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(\alpha + \beta\right) \cdot \frac{\frac{\beta}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + \left(\left(2 \cdot \frac{1}{\alpha} + 8 \cdot \frac{1}{{\alpha}^{3}}\right) - 4 \cdot \frac{1}{{\alpha}^{2}}\right)}{2}\\
\end{array}double f(double alpha, double beta, double i) {
double r115810 = alpha;
double r115811 = beta;
double r115812 = r115810 + r115811;
double r115813 = r115811 - r115810;
double r115814 = r115812 * r115813;
double r115815 = 2.0;
double r115816 = i;
double r115817 = r115815 * r115816;
double r115818 = r115812 + r115817;
double r115819 = r115814 / r115818;
double r115820 = r115818 + r115815;
double r115821 = r115819 / r115820;
double r115822 = 1.0;
double r115823 = r115821 + r115822;
double r115824 = r115823 / r115815;
return r115824;
}
double f(double alpha, double beta, double i) {
double r115825 = alpha;
double r115826 = 6.191224408848943e+177;
bool r115827 = r115825 <= r115826;
double r115828 = beta;
double r115829 = r115825 + r115828;
double r115830 = 2.0;
double r115831 = i;
double r115832 = r115830 * r115831;
double r115833 = r115829 + r115832;
double r115834 = r115828 / r115833;
double r115835 = r115833 + r115830;
double r115836 = r115834 / r115835;
double r115837 = r115829 * r115836;
double r115838 = r115825 / r115833;
double r115839 = r115838 / r115835;
double r115840 = -r115839;
double r115841 = r115840 * r115829;
double r115842 = 1.0;
double r115843 = r115841 + r115842;
double r115844 = 3.0;
double r115845 = pow(r115843, r115844);
double r115846 = cbrt(r115845);
double r115847 = r115837 + r115846;
double r115848 = r115847 / r115830;
double r115849 = 1.0;
double r115850 = r115849 / r115825;
double r115851 = r115830 * r115850;
double r115852 = 8.0;
double r115853 = pow(r115825, r115844);
double r115854 = r115849 / r115853;
double r115855 = r115852 * r115854;
double r115856 = r115851 + r115855;
double r115857 = 4.0;
double r115858 = 2.0;
double r115859 = pow(r115825, r115858);
double r115860 = r115849 / r115859;
double r115861 = r115857 * r115860;
double r115862 = r115856 - r115861;
double r115863 = r115837 + r115862;
double r115864 = r115863 / r115830;
double r115865 = r115827 ? r115848 : r115864;
return r115865;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if alpha < 6.191224408848943e+177Initial program 17.9
rmApplied *-un-lft-identity17.9
Applied *-un-lft-identity17.9
Applied times-frac6.7
Applied times-frac6.6
Simplified6.6
rmApplied div-sub6.6
Applied div-sub6.6
rmApplied sub-neg6.6
Applied distribute-lft-in6.6
Applied associate-+l+6.5
Simplified6.5
rmApplied add-cbrt-cube6.5
Simplified6.5
if 6.191224408848943e+177 < alpha Initial program 64.0
rmApplied *-un-lft-identity64.0
Applied *-un-lft-identity64.0
Applied times-frac49.0
Applied times-frac49.0
Simplified49.0
rmApplied div-sub49.0
Applied div-sub49.0
rmApplied sub-neg49.0
Applied distribute-lft-in49.0
Applied associate-+l+48.6
Simplified48.6
Taylor expanded around inf 39.3
Final simplification11.2
herbie shell --seed 2020049
(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))