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



Bits error versus x
Results
if x < -10771490727067154 or 268163.125619068218 < x Initial program Error: 59.9 bits
rmApplied frac-subError: 62.2 bits
SimplifiedError: 62.2 bits
Taylor expanded around inf Error: 0.3 bits
SimplifiedError: 0.0 bits
if -10771490727067154 < x < 268163.125619068218Initial program Error: 0.6 bits
rmApplied frac-subError: 0.6 bits
SimplifiedError: 0.6 bits
rmApplied distribute-lft-inError: 0.6 bits
Applied associate--r+Error: 0.6 bits
SimplifiedError: 0.0 bits
Final simplificationError: 0.0 bits
herbie shell --seed 2020203
(FPCore (x)
:name "Asymptote C"
:precision binary64
(- (/ x (+ x 1.0)) (/ (+ x 1.0) (- x 1.0))))