\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2.0} + 1.0}{2.0}\begin{array}{l}
\mathbf{if}\;\alpha \le 139510554198006.0:\\
\;\;\;\;{e}^{\left(\log \left(\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0\right)}{2.0}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\sqrt[3]{\beta} \cdot \sqrt[3]{\beta}}{\sqrt[3]{\left(\alpha + \beta\right) + 2.0} \cdot \sqrt[3]{\left(\alpha + \beta\right) + 2.0}} \cdot \frac{\sqrt[3]{\beta}}{\sqrt[3]{\left(\alpha + \beta\right) + 2.0}} - \left(\frac{4.0}{\alpha \cdot \alpha} - \left(\frac{2.0}{\alpha} + \frac{\frac{8.0}{\alpha \cdot \alpha}}{\alpha}\right)\right)}{2.0}\\
\end{array}double f(double alpha, double beta) {
double r1576913 = beta;
double r1576914 = alpha;
double r1576915 = r1576913 - r1576914;
double r1576916 = r1576914 + r1576913;
double r1576917 = 2.0;
double r1576918 = r1576916 + r1576917;
double r1576919 = r1576915 / r1576918;
double r1576920 = 1.0;
double r1576921 = r1576919 + r1576920;
double r1576922 = r1576921 / r1576917;
return r1576922;
}
double f(double alpha, double beta) {
double r1576923 = alpha;
double r1576924 = 139510554198006.0;
bool r1576925 = r1576923 <= r1576924;
double r1576926 = exp(1.0);
double r1576927 = beta;
double r1576928 = r1576923 + r1576927;
double r1576929 = 2.0;
double r1576930 = r1576928 + r1576929;
double r1576931 = r1576927 / r1576930;
double r1576932 = r1576923 / r1576930;
double r1576933 = 1.0;
double r1576934 = r1576932 - r1576933;
double r1576935 = r1576931 - r1576934;
double r1576936 = r1576935 / r1576929;
double r1576937 = log(r1576936);
double r1576938 = pow(r1576926, r1576937);
double r1576939 = cbrt(r1576927);
double r1576940 = r1576939 * r1576939;
double r1576941 = cbrt(r1576930);
double r1576942 = r1576941 * r1576941;
double r1576943 = r1576940 / r1576942;
double r1576944 = r1576939 / r1576941;
double r1576945 = r1576943 * r1576944;
double r1576946 = 4.0;
double r1576947 = r1576923 * r1576923;
double r1576948 = r1576946 / r1576947;
double r1576949 = r1576929 / r1576923;
double r1576950 = 8.0;
double r1576951 = r1576950 / r1576947;
double r1576952 = r1576951 / r1576923;
double r1576953 = r1576949 + r1576952;
double r1576954 = r1576948 - r1576953;
double r1576955 = r1576945 - r1576954;
double r1576956 = r1576955 / r1576929;
double r1576957 = r1576925 ? r1576938 : r1576956;
return r1576957;
}



Bits error versus alpha



Bits error versus beta
Results
if alpha < 139510554198006.0Initial program 0.3
rmApplied div-sub0.3
Applied associate-+l-0.3
rmApplied add-exp-log0.3
rmApplied pow10.3
Applied log-pow0.3
Applied exp-prod0.3
Simplified0.3
if 139510554198006.0 < alpha Initial program 50.2
rmApplied div-sub50.2
Applied associate-+l-48.6
rmApplied add-cube-cbrt48.7
Applied add-cube-cbrt48.6
Applied times-frac48.6
Taylor expanded around inf 18.4
Simplified18.4
Final simplification6.0
herbie shell --seed 2019154 +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))