\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.0} + 1.0}{2.0}\begin{array}{l}
\mathbf{if}\;\alpha \le 2.778527467618903 \cdot 10^{+23}:\\
\;\;\;\;\frac{\sqrt[3]{\left(\left(\beta + \alpha\right) \cdot \frac{\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)}}{2.0 + \left(i \cdot 2 + \left(\beta + \alpha\right)\right)} + 1.0\right) \cdot \left(\left(\left(\beta + \alpha\right) \cdot \frac{\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)}}{2.0 + \left(i \cdot 2 + \left(\beta + \alpha\right)\right)} + 1.0\right) \cdot \left(\left(\beta + \alpha\right) \cdot \frac{\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)}}{2.0 + \left(i \cdot 2 + \left(\beta + \alpha\right)\right)} + 1.0\right)\right)}}{2.0}\\
\mathbf{elif}\;\alpha \le 7.90075049213558 \cdot 10^{+53}:\\
\;\;\;\;\frac{\frac{\frac{\frac{8.0}{\alpha} - 4.0}{\alpha}}{\alpha} + \frac{2.0}{\alpha}}{2.0}\\
\mathbf{elif}\;\alpha \le 4.755359132529547 \cdot 10^{+94}:\\
\;\;\;\;\frac{\left(\beta + \alpha\right) \cdot \frac{\frac{\frac{\beta - \alpha}{i \cdot 2 + \left(\beta + \alpha\right)}}{\sqrt{2.0 + \left(i \cdot 2 + \left(\beta + \alpha\right)\right)}}}{\sqrt{2.0 + \left(i \cdot 2 + \left(\beta + \alpha\right)\right)}} + 1.0}{2.0}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{\frac{8.0}{\alpha} - 4.0}{\alpha}}{\alpha} + \frac{2.0}{\alpha}}{2.0}\\
\end{array}double f(double alpha, double beta, double i) {
double r32217685 = alpha;
double r32217686 = beta;
double r32217687 = r32217685 + r32217686;
double r32217688 = r32217686 - r32217685;
double r32217689 = r32217687 * r32217688;
double r32217690 = 2.0;
double r32217691 = i;
double r32217692 = r32217690 * r32217691;
double r32217693 = r32217687 + r32217692;
double r32217694 = r32217689 / r32217693;
double r32217695 = 2.0;
double r32217696 = r32217693 + r32217695;
double r32217697 = r32217694 / r32217696;
double r32217698 = 1.0;
double r32217699 = r32217697 + r32217698;
double r32217700 = r32217699 / r32217695;
return r32217700;
}
double f(double alpha, double beta, double i) {
double r32217701 = alpha;
double r32217702 = 2.778527467618903e+23;
bool r32217703 = r32217701 <= r32217702;
double r32217704 = beta;
double r32217705 = r32217704 + r32217701;
double r32217706 = r32217704 - r32217701;
double r32217707 = i;
double r32217708 = 2.0;
double r32217709 = r32217707 * r32217708;
double r32217710 = r32217709 + r32217705;
double r32217711 = r32217706 / r32217710;
double r32217712 = 2.0;
double r32217713 = r32217712 + r32217710;
double r32217714 = r32217711 / r32217713;
double r32217715 = r32217705 * r32217714;
double r32217716 = 1.0;
double r32217717 = r32217715 + r32217716;
double r32217718 = r32217717 * r32217717;
double r32217719 = r32217717 * r32217718;
double r32217720 = cbrt(r32217719);
double r32217721 = r32217720 / r32217712;
double r32217722 = 7.90075049213558e+53;
bool r32217723 = r32217701 <= r32217722;
double r32217724 = 8.0;
double r32217725 = r32217724 / r32217701;
double r32217726 = 4.0;
double r32217727 = r32217725 - r32217726;
double r32217728 = r32217727 / r32217701;
double r32217729 = r32217728 / r32217701;
double r32217730 = r32217712 / r32217701;
double r32217731 = r32217729 + r32217730;
double r32217732 = r32217731 / r32217712;
double r32217733 = 4.755359132529547e+94;
bool r32217734 = r32217701 <= r32217733;
double r32217735 = sqrt(r32217713);
double r32217736 = r32217711 / r32217735;
double r32217737 = r32217736 / r32217735;
double r32217738 = r32217705 * r32217737;
double r32217739 = r32217738 + r32217716;
double r32217740 = r32217739 / r32217712;
double r32217741 = r32217734 ? r32217740 : r32217732;
double r32217742 = r32217723 ? r32217732 : r32217741;
double r32217743 = r32217703 ? r32217721 : r32217742;
return r32217743;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if alpha < 2.778527467618903e+23Initial program 11.2
rmApplied *-un-lft-identity11.2
Applied times-frac0.4
Simplified0.4
rmApplied *-un-lft-identity0.4
Applied times-frac0.4
Simplified0.4
rmApplied add-cbrt-cube0.4
if 2.778527467618903e+23 < alpha < 7.90075049213558e+53 or 4.755359132529547e+94 < alpha Initial program 53.7
rmApplied *-un-lft-identity53.7
Applied times-frac39.8
Simplified39.8
rmApplied *-un-lft-identity39.8
Applied times-frac39.8
Simplified39.8
Taylor expanded around -inf 42.0
Simplified42.0
if 7.90075049213558e+53 < alpha < 4.755359132529547e+94Initial program 37.6
rmApplied *-un-lft-identity37.6
Applied times-frac26.6
Simplified26.6
rmApplied *-un-lft-identity26.6
Applied times-frac26.5
Simplified26.5
rmApplied add-sqr-sqrt26.5
Applied associate-/r*26.5
Final simplification12.4
herbie shell --seed 2019107
(FPCore (alpha beta i)
:name "Octave 3.8, jcobi/2"
:pre (and (> alpha -1) (> beta -1) (> i 0))
(/ (+ (/ (/ (* (+ alpha beta) (- beta alpha)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) 2.0)) 1.0) 2.0))