\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 6.301800243409153329403242826949981708813 \cdot 10^{165}:\\
\;\;\;\;\frac{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\alpha + \beta\right) + 2 \cdot i} \cdot \sqrt{\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}} \cdot \frac{\frac{\sqrt{\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 r105036 = i;
double r105037 = alpha;
double r105038 = beta;
double r105039 = r105037 + r105038;
double r105040 = r105039 + r105036;
double r105041 = r105036 * r105040;
double r105042 = r105038 * r105037;
double r105043 = r105042 + r105041;
double r105044 = r105041 * r105043;
double r105045 = 2.0;
double r105046 = r105045 * r105036;
double r105047 = r105039 + r105046;
double r105048 = r105047 * r105047;
double r105049 = r105044 / r105048;
double r105050 = 1.0;
double r105051 = r105048 - r105050;
double r105052 = r105049 / r105051;
return r105052;
}
double f(double alpha, double beta, double i) {
double r105053 = alpha;
double r105054 = 6.301800243409153e+165;
bool r105055 = r105053 <= r105054;
double r105056 = i;
double r105057 = beta;
double r105058 = r105053 + r105057;
double r105059 = r105058 + r105056;
double r105060 = r105056 * r105059;
double r105061 = 2.0;
double r105062 = r105061 * r105056;
double r105063 = r105058 + r105062;
double r105064 = r105060 / r105063;
double r105065 = r105057 * r105053;
double r105066 = r105065 + r105060;
double r105067 = sqrt(r105066);
double r105068 = r105064 * r105067;
double r105069 = 1.0;
double r105070 = sqrt(r105069);
double r105071 = r105063 + r105070;
double r105072 = r105068 / r105071;
double r105073 = r105067 / r105063;
double r105074 = r105063 - r105070;
double r105075 = r105073 / r105074;
double r105076 = r105072 * r105075;
double r105077 = 0.0;
double r105078 = r105055 ? r105076 : r105077;
return r105078;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if alpha < 6.301800243409153e+165Initial program 52.0
rmApplied add-sqr-sqrt52.0
Applied difference-of-squares52.0
Applied times-frac36.8
Applied times-frac35.0
rmApplied *-un-lft-identity35.0
Applied *-un-lft-identity35.0
Applied add-sqr-sqrt35.0
Applied times-frac35.0
Applied times-frac35.1
Applied associate-*r*35.1
Simplified35.1
if 6.301800243409153e+165 < alpha Initial program 64.0
Taylor expanded around inf 47.6
Final simplification36.9
herbie shell --seed 2020002
(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)))