\frac{x}{x + 1} - \frac{x + 1}{x - 1}\begin{array}{l}
\mathbf{if}\;x \leq -13944.835699645984 \lor \neg \left(x \leq 13458.95778492134\right):\\
\;\;\;\;\frac{-1}{x \cdot x} - \left(\frac{3}{x} + \frac{3}{{x}^{3}}\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{1}{x + 1} - \frac{x + 1}{x + -1}\\
\end{array}(FPCore (x) :precision binary64 (- (/ x (+ x 1.0)) (/ (+ x 1.0) (- x 1.0))))
(FPCore (x) :precision binary64 (if (or (<= x -13944.835699645984) (not (<= x 13458.95778492134))) (- (/ -1.0 (* x x)) (+ (/ 3.0 x) (/ 3.0 (pow x 3.0)))) (- (* x (/ 1.0 (+ x 1.0))) (/ (+ x 1.0) (+ x -1.0)))))
double code(double x) {
return (x / (x + 1.0)) - ((x + 1.0) / (x - 1.0));
}
double code(double x) {
double tmp;
if ((x <= -13944.835699645984) || !(x <= 13458.95778492134)) {
tmp = (-1.0 / (x * x)) - ((3.0 / x) + (3.0 / pow(x, 3.0)));
} else {
tmp = (x * (1.0 / (x + 1.0))) - ((x + 1.0) / (x + -1.0));
}
return tmp;
}



Bits error versus x
Results
if x < -13944.8 or 13459 < x Initial program 59.2
Taylor expanded around inf 0.3
Simplified0.0
if -13944.8 < x < 13459Initial program 0.1
rmApplied div-inv_binary64_5750.1
Final simplification0.0
herbie shell --seed 2020231
(FPCore (x)
:name "Asymptote C"
:precision binary64
(- (/ x (+ x 1.0)) (/ (+ x 1.0) (- x 1.0))))