\frac{x}{x + 1} - \frac{x + 1}{x - 1}\begin{array}{l}
\mathbf{if}\;x \leq -14614.459261386228 \lor \neg \left(x \leq 13111.01440601656\right):\\
\;\;\;\;\frac{-1}{x \cdot x} - \left(\frac{3}{x} + \frac{3}{{x}^{3}}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{{x}^{3} \cdot {\left(x + -1\right)}^{3} - {\left(x + 1\right)}^{6}}{{\left(-1 + x \cdot x\right)}^{3}}}{\frac{x}{x + 1} \cdot \frac{x}{x + 1} + \frac{x + \left(x + 1\right) \cdot \frac{x + 1}{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 -14614.459261386228) (not (<= x 13111.01440601656)))
(- (/ -1.0 (* x x)) (+ (/ 3.0 x) (/ 3.0 (pow x 3.0))))
(/
(/
(- (* (pow x 3.0) (pow (+ x -1.0) 3.0)) (pow (+ x 1.0) 6.0))
(pow (+ -1.0 (* x x)) 3.0))
(+
(* (/ x (+ x 1.0)) (/ x (+ x 1.0)))
(/ (+ x (* (+ 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 <= -14614.459261386228) || !(x <= 13111.01440601656)) {
tmp = (-1.0 / (x * x)) - ((3.0 / x) + (3.0 / pow(x, 3.0)));
} else {
tmp = (((pow(x, 3.0) * pow((x + -1.0), 3.0)) - pow((x + 1.0), 6.0)) / pow((-1.0 + (x * x)), 3.0)) / (((x / (x + 1.0)) * (x / (x + 1.0))) + ((x + ((x + 1.0) * ((x + 1.0) / (x + -1.0)))) / (x + -1.0)));
}
return tmp;
}



Bits error versus x
Results
if x < -14614.4592613862278 or 13111.0144060165603 < x Initial program 59.4
Taylor expanded around inf 0.3
Simplified0.0
if -14614.4592613862278 < x < 13111.0144060165603Initial program 0.1
rmApplied flip3--_binary64_14590.1
Simplified0.1
rmApplied cube-div_binary64_14840.1
Applied cube-div_binary64_14840.1
Applied frac-sub_binary64_14640.1
Simplified0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2020292
(FPCore (x)
:name "Asymptote C"
:precision binary64
(- (/ x (+ x 1.0)) (/ (+ x 1.0) (- x 1.0))))