\frac{n \cdot x - y \cdot y}{n \cdot n - n}\frac{n \cdot x - y \cdot y}{n \cdot n - n}double code(double n, double x, double y) {
return ((double) (((double) (((double) (n * x)) - ((double) (y * y)))) / ((double) (((double) (n * n)) - n))));
}
double code(double n, double x, double y) {
return ((double) (((double) (((double) (n * x)) - ((double) (y * y)))) / ((double) (((double) (n * n)) - n))));
}



Bits error versus n



Bits error versus x



Bits error versus y
Results
Initial program 22.2
Final simplification22.2
herbie shell --seed 2020152
(FPCore (n x y)
:name "(/ (- (* n x) (* y y)) (- (* n n) n))"
:precision binary64
(/ (- (* n x) (* y y)) (- (* n n) n)))