\frac{x}{x + 1} - \frac{x + 1}{x - 1}\begin{array}{l}
\mathbf{if}\;x \le -8404.501425772927:\\
\;\;\;\;\frac{-3}{x} + \left(\frac{-1}{x \cdot x} + \frac{\frac{-3}{x}}{x \cdot x}\right)\\
\mathbf{elif}\;x \le 6607.223075937772:\\
\;\;\;\;\frac{\frac{x}{\sqrt[3]{1 + x} \cdot \sqrt[3]{1 + x}}}{\sqrt[3]{1 + x}} - \frac{1 + x}{x - 1}\\
\mathbf{else}:\\
\;\;\;\;\frac{-3}{x} + \left(\frac{-1}{x \cdot x} + \frac{\frac{-3}{x}}{x \cdot x}\right)\\
\end{array}double f(double x) {
double r4644531 = x;
double r4644532 = 1.0;
double r4644533 = r4644531 + r4644532;
double r4644534 = r4644531 / r4644533;
double r4644535 = r4644531 - r4644532;
double r4644536 = r4644533 / r4644535;
double r4644537 = r4644534 - r4644536;
return r4644537;
}
double f(double x) {
double r4644538 = x;
double r4644539 = -8404.501425772927;
bool r4644540 = r4644538 <= r4644539;
double r4644541 = -3.0;
double r4644542 = r4644541 / r4644538;
double r4644543 = -1.0;
double r4644544 = r4644538 * r4644538;
double r4644545 = r4644543 / r4644544;
double r4644546 = r4644542 / r4644544;
double r4644547 = r4644545 + r4644546;
double r4644548 = r4644542 + r4644547;
double r4644549 = 6607.223075937772;
bool r4644550 = r4644538 <= r4644549;
double r4644551 = 1.0;
double r4644552 = r4644551 + r4644538;
double r4644553 = cbrt(r4644552);
double r4644554 = r4644553 * r4644553;
double r4644555 = r4644538 / r4644554;
double r4644556 = r4644555 / r4644553;
double r4644557 = r4644538 - r4644551;
double r4644558 = r4644552 / r4644557;
double r4644559 = r4644556 - r4644558;
double r4644560 = r4644550 ? r4644559 : r4644548;
double r4644561 = r4644540 ? r4644548 : r4644560;
return r4644561;
}



Bits error versus x
Results
if x < -8404.501425772927 or 6607.223075937772 < x Initial program 59.1
Taylor expanded around inf 0.3
Simplified0.0
Taylor expanded around -inf 0.3
Simplified0.0
if -8404.501425772927 < x < 6607.223075937772Initial program 0.1
rmApplied add-cube-cbrt0.1
Applied associate-/r*0.1
Final simplification0.1
herbie shell --seed 2019138
(FPCore (x)
:name "Asymptote C"
(- (/ x (+ x 1)) (/ (+ x 1) (- x 1))))