\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}\;\alpha \le 3.204294196618706018722370553408534967361 \cdot 10^{184}:\\
\;\;\;\;\left(\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}\right) \cdot \frac{1}{\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 r7829964 = i;
double r7829965 = alpha;
double r7829966 = beta;
double r7829967 = r7829965 + r7829966;
double r7829968 = r7829967 + r7829964;
double r7829969 = r7829964 * r7829968;
double r7829970 = r7829966 * r7829965;
double r7829971 = r7829970 + r7829969;
double r7829972 = r7829969 * r7829971;
double r7829973 = 2.0;
double r7829974 = r7829973 * r7829964;
double r7829975 = r7829967 + r7829974;
double r7829976 = r7829975 * r7829975;
double r7829977 = r7829972 / r7829976;
double r7829978 = 1.0;
double r7829979 = r7829976 - r7829978;
double r7829980 = r7829977 / r7829979;
return r7829980;
}
double f(double alpha, double beta, double i) {
double r7829981 = alpha;
double r7829982 = 3.204294196618706e+184;
bool r7829983 = r7829981 <= r7829982;
double r7829984 = i;
double r7829985 = beta;
double r7829986 = r7829981 + r7829985;
double r7829987 = r7829986 + r7829984;
double r7829988 = r7829984 * r7829987;
double r7829989 = 2.0;
double r7829990 = r7829989 * r7829984;
double r7829991 = r7829986 + r7829990;
double r7829992 = r7829988 / r7829991;
double r7829993 = 1.0;
double r7829994 = sqrt(r7829993);
double r7829995 = r7829991 + r7829994;
double r7829996 = r7829992 / r7829995;
double r7829997 = r7829985 * r7829981;
double r7829998 = r7829997 + r7829988;
double r7829999 = r7829998 / r7829991;
double r7830000 = r7829996 * r7829999;
double r7830001 = 1.0;
double r7830002 = r7829991 - r7829994;
double r7830003 = r7830001 / r7830002;
double r7830004 = r7830000 * r7830003;
double r7830005 = 0.0;
double r7830006 = r7829983 ? r7830004 : r7830005;
return r7830006;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if alpha < 3.204294196618706e+184Initial program 52.6
rmApplied add-sqr-sqrt52.6
Applied difference-of-squares52.6
Applied times-frac37.5
Applied times-frac35.4
rmApplied div-inv35.5
Applied associate-*r*35.5
if 3.204294196618706e+184 < alpha Initial program 64.0
Taylor expanded around inf 45.8
Final simplification36.8
herbie shell --seed 2019173
(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)))