\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.5553661537212927 \cdot 10^{137}:\\
\;\;\;\;\frac{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\alpha + \beta\right) + 2 \cdot i} \cdot \frac{\frac{\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}}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}}\\
\mathbf{elif}\;i \le 5.66991920974387326 \cdot 10^{155}:\\
\;\;\;\;\frac{\frac{i}{\sqrt{\left(\alpha + \beta\right) + 2 \cdot i}} \cdot \frac{\left(\alpha + \beta\right) + i}{\sqrt{\left(\alpha + \beta\right) + 2 \cdot i}}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}} \cdot \frac{i}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - \sqrt{1}}\\
\mathbf{elif}\;i \le 2.8985313523571793 \cdot 10^{168}:\\
\;\;\;\;\frac{\frac{i}{\sqrt{\left(\alpha + \beta\right) + 2 \cdot i}} \cdot \frac{\left(\alpha + \beta\right) + i}{\sqrt{\left(\alpha + \beta\right) + 2 \cdot i}}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}} \cdot \frac{0.25 \cdot \alpha + \left(0.5 \cdot i + 0.25 \cdot \beta\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - \sqrt{1}}\\
\mathbf{elif}\;i \le 5.81785064838450142 \cdot 10^{178}:\\
\;\;\;\;\frac{\frac{i}{\sqrt{\left(\alpha + \beta\right) + 2 \cdot i}} \cdot \frac{\left(\alpha + \beta\right) + i}{\sqrt{\left(\alpha + \beta\right) + 2 \cdot i}}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}} \cdot \frac{i}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - \sqrt{1}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{i}{\sqrt{\left(\alpha + \beta\right) + 2 \cdot i}} \cdot \frac{\left(\alpha + \beta\right) + i}{\sqrt{\left(\alpha + \beta\right) + 2 \cdot i}}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}} \cdot \frac{0.25 \cdot \alpha + \left(0.5 \cdot i + 0.25 \cdot \beta\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - \sqrt{1}}\\
\end{array}double f(double alpha, double beta, double i) {
double r185171 = i;
double r185172 = alpha;
double r185173 = beta;
double r185174 = r185172 + r185173;
double r185175 = r185174 + r185171;
double r185176 = r185171 * r185175;
double r185177 = r185173 * r185172;
double r185178 = r185177 + r185176;
double r185179 = r185176 * r185178;
double r185180 = 2.0;
double r185181 = r185180 * r185171;
double r185182 = r185174 + r185181;
double r185183 = r185182 * r185182;
double r185184 = r185179 / r185183;
double r185185 = 1.0;
double r185186 = r185183 - r185185;
double r185187 = r185184 / r185186;
return r185187;
}
double f(double alpha, double beta, double i) {
double r185188 = i;
double r185189 = 1.5553661537212927e+137;
bool r185190 = r185188 <= r185189;
double r185191 = alpha;
double r185192 = beta;
double r185193 = r185191 + r185192;
double r185194 = r185193 + r185188;
double r185195 = r185188 * r185194;
double r185196 = 2.0;
double r185197 = r185196 * r185188;
double r185198 = r185193 + r185197;
double r185199 = r185195 / r185198;
double r185200 = r185192 * r185191;
double r185201 = r185200 + r185195;
double r185202 = r185201 / r185198;
double r185203 = 1.0;
double r185204 = sqrt(r185203);
double r185205 = r185198 - r185204;
double r185206 = r185202 / r185205;
double r185207 = r185199 * r185206;
double r185208 = r185198 + r185204;
double r185209 = r185207 / r185208;
double r185210 = 5.669919209743873e+155;
bool r185211 = r185188 <= r185210;
double r185212 = sqrt(r185198);
double r185213 = r185188 / r185212;
double r185214 = r185194 / r185212;
double r185215 = r185213 * r185214;
double r185216 = r185215 / r185208;
double r185217 = r185188 / r185205;
double r185218 = r185216 * r185217;
double r185219 = 2.8985313523571793e+168;
bool r185220 = r185188 <= r185219;
double r185221 = 0.25;
double r185222 = r185221 * r185191;
double r185223 = 0.5;
double r185224 = r185223 * r185188;
double r185225 = r185221 * r185192;
double r185226 = r185224 + r185225;
double r185227 = r185222 + r185226;
double r185228 = r185227 / r185205;
double r185229 = r185216 * r185228;
double r185230 = 5.8178506483845014e+178;
bool r185231 = r185188 <= r185230;
double r185232 = r185231 ? r185218 : r185229;
double r185233 = r185220 ? r185229 : r185232;
double r185234 = r185211 ? r185218 : r185233;
double r185235 = r185190 ? r185209 : r185234;
return r185235;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if i < 1.5553661537212927e+137Initial program 41.6
rmApplied add-sqr-sqrt41.6
Applied difference-of-squares41.6
Applied times-frac15.0
Applied times-frac10.7
rmApplied associate-*l/10.7
if 1.5553661537212927e+137 < i < 5.669919209743873e+155 or 2.8985313523571793e+168 < i < 5.8178506483845014e+178Initial program 64.0
rmApplied add-sqr-sqrt64.0
Applied difference-of-squares64.0
Applied times-frac37.6
Applied times-frac37.0
rmApplied add-sqr-sqrt37.2
Applied times-frac37.2
Taylor expanded around inf 39.2
if 5.669919209743873e+155 < i < 2.8985313523571793e+168 or 5.8178506483845014e+178 < i Initial program 64.0
rmApplied add-sqr-sqrt64.0
Applied difference-of-squares64.0
Applied times-frac64.0
Applied times-frac64.0
rmApplied add-sqr-sqrt64.0
Applied times-frac64.0
Taylor expanded around 0 9.6
Final simplification12.9
herbie shell --seed 2020049
(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)))