\frac{\frac{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}
\begin{array}{l}
t_0 := \left(\alpha + \beta\right) + 2\\
t_1 := \frac{1}{t_0}\\
\frac{\frac{\left(\alpha + 1\right) \cdot \left(t_1 + \beta \cdot t_1\right)}{t_0}}{\alpha + \left(\beta + 3\right)}
\end{array}
(FPCore (alpha beta) :precision binary64 (/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2.0 1.0))) (+ (+ alpha beta) (* 2.0 1.0))) (+ (+ (+ alpha beta) (* 2.0 1.0)) 1.0)))
(FPCore (alpha beta) :precision binary64 (let* ((t_0 (+ (+ alpha beta) 2.0)) (t_1 (/ 1.0 t_0))) (/ (/ (* (+ alpha 1.0) (+ t_1 (* beta t_1))) t_0) (+ alpha (+ beta 3.0)))))
double code(double alpha, double beta) {
return (((((alpha + beta) + (beta * alpha)) + 1.0) / ((alpha + beta) + (2.0 * 1.0))) / ((alpha + beta) + (2.0 * 1.0))) / (((alpha + beta) + (2.0 * 1.0)) + 1.0);
}
double code(double alpha, double beta) {
double t_0 = (alpha + beta) + 2.0;
double t_1 = 1.0 / t_0;
return (((alpha + 1.0) * (t_1 + (beta * t_1))) / t_0) / (alpha + (beta + 3.0));
}



Bits error versus alpha



Bits error versus beta
Results
Initial program 3.8
Simplified2.3
Applied *-un-lft-identity_binary642.3
Applied times-frac_binary640.1
Applied associate-*r/_binary640.1
Applied associate-*r/_binary640.1
Applied distribute-lft-in_binary640.1
Final simplification0.1
herbie shell --seed 2021280
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/3"
:precision binary64
:pre (and (> alpha -1.0) (> beta -1.0))
(/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2.0 1.0))) (+ (+ alpha beta) (* 2.0 1.0))) (+ (+ (+ alpha beta) (* 2.0 1.0)) 1.0)))