\sqrt{x + \frac{1}{x}} + \sqrt{x - \frac{1}{x}}\sqrt{x + \frac{1}{x}} + \sqrt{x - \frac{1}{x}}double code(double x) {
return ((double) (((double) sqrt(((double) (x + ((double) (1.0 / x)))))) + ((double) sqrt(((double) (x - ((double) (1.0 / x))))))));
}
double code(double x) {
return ((double) (((double) sqrt(((double) (x + ((double) (1.0 / x)))))) + ((double) sqrt(((double) (x - ((double) (1.0 / x))))))));
}



Bits error versus x
Results
Initial program 0.0
Final simplification0.0
herbie shell --seed 2020153
(FPCore (x)
:name "(+ (sqrt (+ x (/ 1 x))) (sqrt (- x (/ 1 x))))"
:precision binary64
(+ (sqrt (+ x (/ 1.0 x))) (sqrt (- x (/ 1.0 x)))))