\frac{a - \frac{b \cdot b}{n}}{n - 1}\frac{a - \frac{b \cdot b}{n}}{n - 1}double code(double a, double b, double n) {
return ((double) (((double) (a - ((double) (((double) (b * b)) / n)))) / ((double) (n - 1.0))));
}
double code(double a, double b, double n) {
return ((double) (((double) (a - ((double) (((double) (b * b)) / n)))) / ((double) (n - 1.0))));
}



Bits error versus a



Bits error versus b



Bits error versus n
Results
Initial program 8.4
Final simplification8.4
herbie shell --seed 2020153
(FPCore (a b n)
:name "(/ (- a (/ (* b b) n)) (- n 1))"
:precision binary64
(/ (- a (/ (* b b) n)) (- n 1.0)))