\frac{x}{x + 1} - \frac{x + 1}{x - 1}\begin{array}{l}
\mathbf{if}\;x \le -15417.9220111292761 \lor \neg \left(x \le 14949.242297652785\right):\\
\;\;\;\;\frac{-1}{{x}^{2}} - \mathsf{fma}\left(3, \frac{1}{x}, 3 \cdot \frac{1}{{x}^{3}}\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt[3]{{\left(\frac{{\left(\frac{x}{x + 1}\right)}^{3} - {\left(\frac{x + 1}{x - 1}\right)}^{3}}{\mathsf{fma}\left(\frac{x + 1}{x - 1}, \frac{x}{x + 1} + \frac{x + 1}{x - 1}, \frac{x}{x + 1} \cdot \frac{x}{x + 1}\right)}\right)}^{3}}\\
\end{array}double code(double x) {
return ((double) (((double) (x / ((double) (x + 1.0)))) - ((double) (((double) (x + 1.0)) / ((double) (x - 1.0))))));
}
double code(double x) {
double VAR;
if (((x <= -15417.922011129276) || !(x <= 14949.242297652785))) {
VAR = ((double) (((double) (((double) -(1.0)) / ((double) pow(x, 2.0)))) - ((double) fma(3.0, ((double) (1.0 / x)), ((double) (3.0 * ((double) (1.0 / ((double) pow(x, 3.0))))))))));
} else {
VAR = ((double) cbrt(((double) pow(((double) (((double) (((double) pow(((double) (x / ((double) (x + 1.0)))), 3.0)) - ((double) pow(((double) (((double) (x + 1.0)) / ((double) (x - 1.0)))), 3.0)))) / ((double) fma(((double) (((double) (x + 1.0)) / ((double) (x - 1.0)))), ((double) (((double) (x / ((double) (x + 1.0)))) + ((double) (((double) (x + 1.0)) / ((double) (x - 1.0)))))), ((double) (((double) (x / ((double) (x + 1.0)))) * ((double) (x / ((double) (x + 1.0)))))))))), 3.0))));
}
return VAR;
}



Bits error versus x
Results
if x < -15417.922011129276 or 14949.242297652785 < x Initial program 59.5
Taylor expanded around inf 0.3
Simplified0.3
if -15417.922011129276 < x < 14949.242297652785Initial program 0.1
rmApplied flip3--0.1
Simplified0.1
rmApplied add-cbrt-cube0.1
Applied add-cbrt-cube0.1
Applied cbrt-undiv0.1
Simplified0.1
Final simplification0.2
herbie shell --seed 2020121 +o rules:numerics
(FPCore (x)
:name "Asymptote C"
:precision binary64
(- (/ x (+ x 1)) (/ (+ x 1) (- x 1))))