\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}\;\beta \le 2.565406916826700694059349673172180602042 \cdot 10^{215}:\\
\;\;\;\;\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{1}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - \sqrt{1}}{\frac{\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}double f(double alpha, double beta, double i) {
double r95793 = i;
double r95794 = alpha;
double r95795 = beta;
double r95796 = r95794 + r95795;
double r95797 = r95796 + r95793;
double r95798 = r95793 * r95797;
double r95799 = r95795 * r95794;
double r95800 = r95799 + r95798;
double r95801 = r95798 * r95800;
double r95802 = 2.0;
double r95803 = r95802 * r95793;
double r95804 = r95796 + r95803;
double r95805 = r95804 * r95804;
double r95806 = r95801 / r95805;
double r95807 = 1.0;
double r95808 = r95805 - r95807;
double r95809 = r95806 / r95808;
return r95809;
}
double f(double alpha, double beta, double i) {
double r95810 = beta;
double r95811 = 2.5654069168267007e+215;
bool r95812 = r95810 <= r95811;
double r95813 = i;
double r95814 = alpha;
double r95815 = r95814 + r95810;
double r95816 = r95815 + r95813;
double r95817 = r95813 * r95816;
double r95818 = 2.0;
double r95819 = r95818 * r95813;
double r95820 = r95815 + r95819;
double r95821 = r95817 / r95820;
double r95822 = 1.0;
double r95823 = sqrt(r95822);
double r95824 = r95820 + r95823;
double r95825 = r95821 / r95824;
double r95826 = 1.0;
double r95827 = r95820 - r95823;
double r95828 = r95810 * r95814;
double r95829 = r95828 + r95817;
double r95830 = r95829 / r95820;
double r95831 = r95827 / r95830;
double r95832 = r95826 / r95831;
double r95833 = r95825 * r95832;
double r95834 = 0.0;
double r95835 = r95812 ? r95833 : r95834;
return r95835;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if beta < 2.5654069168267007e+215Initial program 53.0
rmApplied add-sqr-sqrt53.0
Applied difference-of-squares53.0
Applied times-frac37.8
Applied times-frac35.5
rmApplied *-un-lft-identity35.5
Applied *-un-lft-identity35.5
Applied times-frac35.5
Applied associate-/l*35.5
if 2.5654069168267007e+215 < beta Initial program 64.0
Taylor expanded around inf 43.0
Final simplification36.2
herbie shell --seed 2020001
(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)))