\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 r5240168 = i;
double r5240169 = alpha;
double r5240170 = beta;
double r5240171 = r5240169 + r5240170;
double r5240172 = r5240171 + r5240168;
double r5240173 = r5240168 * r5240172;
double r5240174 = r5240170 * r5240169;
double r5240175 = r5240174 + r5240173;
double r5240176 = r5240173 * r5240175;
double r5240177 = 2.0;
double r5240178 = r5240177 * r5240168;
double r5240179 = r5240171 + r5240178;
double r5240180 = r5240179 * r5240179;
double r5240181 = r5240176 / r5240180;
double r5240182 = 1.0;
double r5240183 = r5240180 - r5240182;
double r5240184 = r5240181 / r5240183;
return r5240184;
}
double f(double alpha, double beta, double i) {
double r5240185 = i;
double r5240186 = 1.3375206427593509e+154;
bool r5240187 = r5240185 <= r5240186;
double r5240188 = alpha;
double r5240189 = beta;
double r5240190 = r5240188 + r5240189;
double r5240191 = r5240190 + r5240185;
double r5240192 = r5240185 * r5240191;
double r5240193 = r5240188 * r5240189;
double r5240194 = r5240192 + r5240193;
double r5240195 = 2.0;
double r5240196 = r5240185 * r5240195;
double r5240197 = r5240190 + r5240196;
double r5240198 = r5240194 / r5240197;
double r5240199 = 1.0;
double r5240200 = sqrt(r5240199);
double r5240201 = r5240197 - r5240200;
double r5240202 = r5240198 / r5240201;
double r5240203 = sqrt(r5240202);
double r5240204 = r5240192 / r5240197;
double r5240205 = r5240200 + r5240197;
double r5240206 = r5240204 / r5240205;
double r5240207 = r5240206 * r5240203;
double r5240208 = r5240203 * r5240207;
double r5240209 = 0.0;
double r5240210 = r5240187 ? r5240208 : r5240209;
return r5240210;
}



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
(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)))