\log \left(\frac{1}{x} + \frac{\sqrt{1 - x \cdot x}}{x}\right)\log \left(\frac{1}{x} + \frac{\sqrt{1 - x \cdot x}}{x}\right)double code(double x) {
return ((double) log(((double) (((double) (1.0 / x)) + ((double) (((double) sqrt(((double) (1.0 - ((double) (x * x)))))) / x))))));
}
double code(double x) {
return ((double) log(((double) (((double) (1.0 / x)) + ((double) (((double) sqrt(((double) (1.0 - ((double) (x * x)))))) / x))))));
}



Bits error versus x
Results
Initial program 0.1
Final simplification0.1
herbie shell --seed 2020123 +o rules:numerics
(FPCore (x)
:name "Hyperbolic arc-(co)secant"
:precision binary64
(log (+ (/ 1 x) (/ (sqrt (- 1 (* x x))) x))))