\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 40.4
Final simplification40.4
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)))))