\frac{\frac{\frac{\left(\alpha + \beta\right) \cdot \left(\beta - \alpha\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + 1}{2}\begin{array}{l}
\mathbf{if}\;\alpha \le 9.76982660809475993097637758871956020598 \cdot 10^{62}:\\
\;\;\;\;\frac{1 + \sqrt[3]{\left(\left(\beta + \alpha\right) \cdot \frac{\sqrt[3]{\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)} \cdot \left(\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)} \cdot \frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)}\right)}}{2 + \left(i \cdot 2 + \left(\beta + \alpha\right)\right)}\right) \cdot \left(\left(\left(\beta + \alpha\right) \cdot \frac{\sqrt[3]{\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)} \cdot \left(\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)} \cdot \frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)}\right)}}{2 + \left(i \cdot 2 + \left(\beta + \alpha\right)\right)}\right) \cdot \left(\left(\beta + \alpha\right) \cdot \frac{\sqrt[3]{\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)} \cdot \left(\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)} \cdot \frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)}\right)}}{2 + \left(i \cdot 2 + \left(\beta + \alpha\right)\right)}\right)\right)}}{2}\\
\mathbf{elif}\;\alpha \le 2.266980258015373608283608969920590572301 \cdot 10^{166}:\\
\;\;\;\;\frac{\left(\frac{2}{\alpha} - \frac{4}{\alpha \cdot \alpha}\right) + \frac{8}{\alpha \cdot \left(\alpha \cdot \alpha\right)}}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{{\left(\left(\beta + \alpha\right) \cdot \frac{\sqrt[3]{\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)} \cdot \left(\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)} \cdot \frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)}\right)}}{2 + \left(i \cdot 2 + \left(\beta + \alpha\right)\right)}\right)}^{3} + {1}^{3}}{\left(\left(\beta + \alpha\right) \cdot \frac{\sqrt[3]{\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)} \cdot \left(\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)} \cdot \frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)}\right)}}{2 + \left(i \cdot 2 + \left(\beta + \alpha\right)\right)}\right) \cdot \left(\left(\beta + \alpha\right) \cdot \frac{\sqrt[3]{\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)} \cdot \left(\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)} \cdot \frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)}\right)}}{2 + \left(i \cdot 2 + \left(\beta + \alpha\right)\right)}\right) + \left(1 \cdot 1 - \left(\left(\beta + \alpha\right) \cdot \frac{\sqrt[3]{\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)} \cdot \left(\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)} \cdot \frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)}\right)}}{2 + \left(i \cdot 2 + \left(\beta + \alpha\right)\right)}\right) \cdot 1\right)}}{2}\\
\end{array}double f(double alpha, double beta, double i) {
double r4385099 = alpha;
double r4385100 = beta;
double r4385101 = r4385099 + r4385100;
double r4385102 = r4385100 - r4385099;
double r4385103 = r4385101 * r4385102;
double r4385104 = 2.0;
double r4385105 = i;
double r4385106 = r4385104 * r4385105;
double r4385107 = r4385101 + r4385106;
double r4385108 = r4385103 / r4385107;
double r4385109 = r4385107 + r4385104;
double r4385110 = r4385108 / r4385109;
double r4385111 = 1.0;
double r4385112 = r4385110 + r4385111;
double r4385113 = r4385112 / r4385104;
return r4385113;
}
double f(double alpha, double beta, double i) {
double r4385114 = alpha;
double r4385115 = 9.76982660809476e+62;
bool r4385116 = r4385114 <= r4385115;
double r4385117 = 1.0;
double r4385118 = beta;
double r4385119 = r4385118 + r4385114;
double r4385120 = r4385118 - r4385114;
double r4385121 = i;
double r4385122 = 2.0;
double r4385123 = r4385121 * r4385122;
double r4385124 = r4385123 + r4385119;
double r4385125 = r4385120 / r4385124;
double r4385126 = r4385125 * r4385125;
double r4385127 = r4385125 * r4385126;
double r4385128 = cbrt(r4385127);
double r4385129 = r4385122 + r4385124;
double r4385130 = r4385128 / r4385129;
double r4385131 = r4385119 * r4385130;
double r4385132 = r4385131 * r4385131;
double r4385133 = r4385131 * r4385132;
double r4385134 = cbrt(r4385133);
double r4385135 = r4385117 + r4385134;
double r4385136 = r4385135 / r4385122;
double r4385137 = 2.2669802580153736e+166;
bool r4385138 = r4385114 <= r4385137;
double r4385139 = r4385122 / r4385114;
double r4385140 = 4.0;
double r4385141 = r4385114 * r4385114;
double r4385142 = r4385140 / r4385141;
double r4385143 = r4385139 - r4385142;
double r4385144 = 8.0;
double r4385145 = r4385114 * r4385141;
double r4385146 = r4385144 / r4385145;
double r4385147 = r4385143 + r4385146;
double r4385148 = r4385147 / r4385122;
double r4385149 = 3.0;
double r4385150 = pow(r4385131, r4385149);
double r4385151 = pow(r4385117, r4385149);
double r4385152 = r4385150 + r4385151;
double r4385153 = r4385117 * r4385117;
double r4385154 = r4385131 * r4385117;
double r4385155 = r4385153 - r4385154;
double r4385156 = r4385132 + r4385155;
double r4385157 = r4385152 / r4385156;
double r4385158 = r4385157 / r4385122;
double r4385159 = r4385138 ? r4385148 : r4385158;
double r4385160 = r4385116 ? r4385136 : r4385159;
return r4385160;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if alpha < 9.76982660809476e+62Initial program 12.1
rmApplied *-un-lft-identity12.1
Applied *-un-lft-identity12.1
Applied times-frac1.6
Applied times-frac1.6
Simplified1.6
rmApplied add-cbrt-cube1.6
rmApplied add-cbrt-cube1.6
if 9.76982660809476e+62 < alpha < 2.2669802580153736e+166Initial program 45.7
rmApplied *-un-lft-identity45.7
Applied *-un-lft-identity45.7
Applied times-frac32.3
Applied times-frac32.2
Simplified32.2
rmApplied add-log-exp32.2
Taylor expanded around inf 39.9
Simplified39.9
if 2.2669802580153736e+166 < alpha Initial program 64.0
rmApplied *-un-lft-identity64.0
Applied *-un-lft-identity64.0
Applied times-frac47.4
Applied times-frac47.4
Simplified47.4
rmApplied add-cbrt-cube47.4
rmApplied flip3-+47.5
Final simplification13.2
herbie shell --seed 2019179
(FPCore (alpha beta i)
:name "Octave 3.8, jcobi/2"
:pre (and (> alpha -1.0) (> beta -1.0) (> i 0.0))
(/ (+ (/ (/ (* (+ alpha beta) (- beta alpha)) (+ (+ alpha beta) (* 2.0 i))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)) 1.0) 2.0))