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



Bits error versus x
Results
if x < -13286.958638175067 or 11564.8798436448142 < x Initial program 59.3
Taylor expanded around inf 0.3
Simplified0.3
if -13286.958638175067 < x < 11564.8798436448142Initial program 0.1
rmApplied frac-sub0.1
Simplified0.1
Final simplification0.2
herbie shell --seed 2020198
(FPCore (x)
:name "Asymptote C"
:precision binary64
(- (/ x (+ x 1.0)) (/ (+ x 1.0) (- x 1.0))))