\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}\;i \le 1.337520642759350918795226553730261611952 \cdot 10^{154}:\\
\;\;\;\;\sqrt{\frac{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right) + \alpha \cdot \beta}{\left(\alpha + \beta\right) + i \cdot 2}}{\left(\left(\alpha + \beta\right) + i \cdot 2\right) - \sqrt{1}}} \cdot \left(\frac{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\alpha + \beta\right) + i \cdot 2}}{\sqrt{1} + \left(\left(\alpha + \beta\right) + i \cdot 2\right)} \cdot \sqrt{\frac{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right) + \alpha \cdot \beta}{\left(\alpha + \beta\right) + i \cdot 2}}{\left(\left(\alpha + \beta\right) + i \cdot 2\right) - \sqrt{1}}}\right)\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}double f(double alpha, double beta, double i) {
double r4104882 = i;
double r4104883 = alpha;
double r4104884 = beta;
double r4104885 = r4104883 + r4104884;
double r4104886 = r4104885 + r4104882;
double r4104887 = r4104882 * r4104886;
double r4104888 = r4104884 * r4104883;
double r4104889 = r4104888 + r4104887;
double r4104890 = r4104887 * r4104889;
double r4104891 = 2.0;
double r4104892 = r4104891 * r4104882;
double r4104893 = r4104885 + r4104892;
double r4104894 = r4104893 * r4104893;
double r4104895 = r4104890 / r4104894;
double r4104896 = 1.0;
double r4104897 = r4104894 - r4104896;
double r4104898 = r4104895 / r4104897;
return r4104898;
}
double f(double alpha, double beta, double i) {
double r4104899 = i;
double r4104900 = 1.3375206427593509e+154;
bool r4104901 = r4104899 <= r4104900;
double r4104902 = alpha;
double r4104903 = beta;
double r4104904 = r4104902 + r4104903;
double r4104905 = r4104904 + r4104899;
double r4104906 = r4104899 * r4104905;
double r4104907 = r4104902 * r4104903;
double r4104908 = r4104906 + r4104907;
double r4104909 = 2.0;
double r4104910 = r4104899 * r4104909;
double r4104911 = r4104904 + r4104910;
double r4104912 = r4104908 / r4104911;
double r4104913 = 1.0;
double r4104914 = sqrt(r4104913);
double r4104915 = r4104911 - r4104914;
double r4104916 = r4104912 / r4104915;
double r4104917 = sqrt(r4104916);
double r4104918 = r4104906 / r4104911;
double r4104919 = r4104914 + r4104911;
double r4104920 = r4104918 / r4104919;
double r4104921 = r4104920 * r4104917;
double r4104922 = r4104917 * r4104921;
double r4104923 = 0.0;
double r4104924 = r4104901 ? r4104922 : r4104923;
return r4104924;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if i < 1.3375206427593509e+154Initial program 44.1
rmApplied add-sqr-sqrt44.1
Applied difference-of-squares44.1
Applied times-frac15.6
Applied times-frac10.9
rmApplied add-sqr-sqrt11.0
Applied associate-*r*11.0
if 1.3375206427593509e+154 < i Initial program 64.0
Taylor expanded around inf 61.9
Final simplification36.8
herbie shell --seed 2019174 +o rules:numerics
(FPCore (alpha beta i)
:name "Octave 3.8, jcobi/4"
:pre (and (> alpha -1.0) (> beta -1.0) (> i 1.0))
(/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i)))) (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0)))