\frac{c + \sqrt{c \cdot c + 4}}{2} - x\frac{c + \sqrt{c \cdot c + 4}}{2} - xdouble code(double c, double x) {
return ((double) (((double) (((double) (c + ((double) sqrt(((double) (((double) (c * c)) + 4.0)))))) / 2.0)) - x));
}
double code(double c, double x) {
return ((double) (((double) (((double) (c + ((double) sqrt(((double) (((double) (c * c)) + 4.0)))))) / 2.0)) - x));
}



Bits error versus c



Bits error versus x
Results
Initial program 18.3
Final simplification18.3
herbie shell --seed 2020153
(FPCore (c x)
:name "(- (/ (+ c (sqrt (+ (* c c) 4))) 2) x)"
:precision binary64
(- (/ (+ c (sqrt (+ (* c c) 4.0))) 2.0) x))