\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}\;\beta \le 2.35404421953347644 \cdot 10^{209}:\\
\;\;\;\;1 \cdot \frac{\frac{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}} \cdot \frac{\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\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 r208838 = i;
double r208839 = alpha;
double r208840 = beta;
double r208841 = r208839 + r208840;
double r208842 = r208841 + r208838;
double r208843 = r208838 * r208842;
double r208844 = r208840 * r208839;
double r208845 = r208844 + r208843;
double r208846 = r208843 * r208845;
double r208847 = 2.0;
double r208848 = r208847 * r208838;
double r208849 = r208841 + r208848;
double r208850 = r208849 * r208849;
double r208851 = r208846 / r208850;
double r208852 = 1.0;
double r208853 = r208850 - r208852;
double r208854 = r208851 / r208853;
return r208854;
}
double f(double alpha, double beta, double i) {
double r208855 = beta;
double r208856 = 2.3540442195334764e+209;
bool r208857 = r208855 <= r208856;
double r208858 = 1.0;
double r208859 = i;
double r208860 = alpha;
double r208861 = r208860 + r208855;
double r208862 = r208861 + r208859;
double r208863 = r208859 * r208862;
double r208864 = 2.0;
double r208865 = r208864 * r208859;
double r208866 = r208861 + r208865;
double r208867 = r208863 / r208866;
double r208868 = 1.0;
double r208869 = sqrt(r208868);
double r208870 = r208866 + r208869;
double r208871 = r208867 / r208870;
double r208872 = r208855 * r208860;
double r208873 = r208872 + r208863;
double r208874 = r208873 / r208866;
double r208875 = r208871 * r208874;
double r208876 = r208866 - r208869;
double r208877 = r208875 / r208876;
double r208878 = r208858 * r208877;
double r208879 = 0.0;
double r208880 = r208857 ? r208878 : r208879;
return r208880;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if beta < 2.3540442195334764e+209Initial program 53.0
rmApplied add-sqr-sqrt53.0
Applied difference-of-squares53.0
Applied times-frac37.9
Applied times-frac35.5
rmApplied associate-*r/35.4
rmApplied *-un-lft-identity35.4
if 2.3540442195334764e+209 < beta Initial program 64.0
rmApplied add-sqr-sqrt64.0
Applied difference-of-squares64.0
Applied times-frac57.0
Applied times-frac56.1
rmApplied associate-*r/56.1
rmApplied *-un-lft-identity56.1
Taylor expanded around inf 43.2
Final simplification36.3
herbie shell --seed 2020018
(FPCore (alpha beta i)
:name "Octave 3.8, jcobi/4"
:precision binary64
:pre (and (> alpha -1) (> beta -1) (> i 1))
(/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1)))