\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2.0} + 1.0}{2.0}\begin{array}{l}
\mathbf{if}\;\alpha \le 5.669431488829545 \cdot 10^{+27}:\\
\;\;\;\;{e}^{\left(\log \left(\frac{\frac{\beta}{\left(2.0 + \alpha\right) + \beta} - \left(\frac{\alpha}{\left(2.0 + \alpha\right) + \beta} - 1.0\right)}{2.0}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(\frac{4.0}{\alpha \cdot \alpha} - \left(\frac{2.0}{\alpha} + \frac{\frac{8.0}{\alpha}}{\alpha \cdot \alpha}\right)\right)}{2.0}\\
\end{array}double f(double alpha, double beta) {
double r3566998 = beta;
double r3566999 = alpha;
double r3567000 = r3566998 - r3566999;
double r3567001 = r3566999 + r3566998;
double r3567002 = 2.0;
double r3567003 = r3567001 + r3567002;
double r3567004 = r3567000 / r3567003;
double r3567005 = 1.0;
double r3567006 = r3567004 + r3567005;
double r3567007 = r3567006 / r3567002;
return r3567007;
}
double f(double alpha, double beta) {
double r3567008 = alpha;
double r3567009 = 5.669431488829545e+27;
bool r3567010 = r3567008 <= r3567009;
double r3567011 = exp(1.0);
double r3567012 = beta;
double r3567013 = 2.0;
double r3567014 = r3567013 + r3567008;
double r3567015 = r3567014 + r3567012;
double r3567016 = r3567012 / r3567015;
double r3567017 = r3567008 / r3567015;
double r3567018 = 1.0;
double r3567019 = r3567017 - r3567018;
double r3567020 = r3567016 - r3567019;
double r3567021 = r3567020 / r3567013;
double r3567022 = log(r3567021);
double r3567023 = pow(r3567011, r3567022);
double r3567024 = r3567008 + r3567012;
double r3567025 = r3567024 + r3567013;
double r3567026 = r3567012 / r3567025;
double r3567027 = 4.0;
double r3567028 = r3567008 * r3567008;
double r3567029 = r3567027 / r3567028;
double r3567030 = r3567013 / r3567008;
double r3567031 = 8.0;
double r3567032 = r3567031 / r3567008;
double r3567033 = r3567032 / r3567028;
double r3567034 = r3567030 + r3567033;
double r3567035 = r3567029 - r3567034;
double r3567036 = r3567026 - r3567035;
double r3567037 = r3567036 / r3567013;
double r3567038 = r3567010 ? r3567023 : r3567037;
return r3567038;
}



Bits error versus alpha



Bits error versus beta
Results
if alpha < 5.669431488829545e+27Initial program 1.2
rmApplied div-sub1.2
Applied associate-+l-1.2
rmApplied add-exp-log1.2
Applied add-exp-log1.2
Applied div-exp1.2
Simplified1.2
rmApplied pow11.2
Applied log-pow1.2
Applied exp-prod1.2
Simplified1.2
rmApplied *-un-lft-identity1.2
Applied *-un-lft-identity1.2
Applied times-frac1.2
Simplified1.2
if 5.669431488829545e+27 < alpha Initial program 50.6
rmApplied div-sub50.6
Applied associate-+l-49.0
Taylor expanded around inf 17.8
Simplified17.8
Final simplification6.2
herbie shell --seed 2019165 +o rules:numerics
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/1"
:pre (and (> alpha -1) (> beta -1))
(/ (+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0) 2.0))