\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 \leq 5.5059367458921626 \cdot 10^{+131}:\\
\;\;\;\;\frac{\frac{i \cdot \left(\alpha + \left(i + \beta\right)\right)}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \frac{\frac{i \cdot \left(\alpha + \left(i + \beta\right)\right) + \alpha \cdot \beta}{\alpha + \left(\beta + i \cdot 2\right)}}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) - 1}}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) + 1}\\
\mathbf{else}:\\
\;\;\;\;0.0625\\
\end{array}(FPCore (alpha beta i) :precision binary64 (/ (/ (* (* 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)))
(FPCore (alpha beta i)
:precision binary64
(if (<= i 5.5059367458921626e+131)
(/
(*
(/ (* i (+ alpha (+ i beta))) (+ alpha (+ beta (* i 2.0))))
(/
(/
(+ (* i (+ alpha (+ i beta))) (* alpha beta))
(+ alpha (+ beta (* i 2.0))))
(- (+ alpha (+ beta (* i 2.0))) 1.0)))
(+ (+ alpha (+ beta (* i 2.0))) 1.0))
0.0625))double code(double alpha, double beta, double i) {
return (((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);
}
double code(double alpha, double beta, double i) {
double tmp;
if (i <= 5.5059367458921626e+131) {
tmp = (((i * (alpha + (i + beta))) / (alpha + (beta + (i * 2.0)))) * ((((i * (alpha + (i + beta))) + (alpha * beta)) / (alpha + (beta + (i * 2.0)))) / ((alpha + (beta + (i * 2.0))) - 1.0))) / ((alpha + (beta + (i * 2.0))) + 1.0);
} else {
tmp = 0.0625;
}
return tmp;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if i < 5.5059367458921626e131Initial program 40.1
rmApplied difference-of-sqr-1_binary6440.1
Applied times-frac_binary6414.5
Applied times-frac_binary6410.0
Simplified10.0
Simplified10.0
rmApplied associate-*l/_binary6410.0
if 5.5059367458921626e131 < i Initial program 64.0
Taylor expanded around inf 11.2
Final simplification10.7
herbie shell --seed 2021173
(FPCore (alpha beta i)
:name "Octave 3.8, jcobi/4"
:precision binary64
: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)))