double f(double x) {
double r13751295 = x;
double r13751296 = 1.0;
double r13751297 = r13751295 + r13751296;
double r13751298 = r13751295 / r13751297;
double r13751299 = r13751295 - r13751296;
double r13751300 = r13751297 / r13751299;
double r13751301 = r13751298 - r13751300;
return r13751301;
}
double f(double x) {
double r13751302 = x;
double r13751303 = -0.9912853724188615;
bool r13751304 = r13751302 <= r13751303;
double r13751305 = -1.0;
double r13751306 = r13751302 * r13751302;
double r13751307 = r13751305 / r13751306;
double r13751308 = -3.0;
double r13751309 = r13751306 * r13751302;
double r13751310 = r13751308 / r13751309;
double r13751311 = r13751308 / r13751302;
double r13751312 = r13751310 + r13751311;
double r13751313 = r13751307 + r13751312;
double r13751314 = 12935.179381933913;
bool r13751315 = r13751302 <= r13751314;
double r13751316 = 1.0;
double r13751317 = r13751316 + r13751302;
double r13751318 = r13751302 / r13751317;
double r13751319 = sqrt(r13751317);
double r13751320 = r13751302 - r13751316;
double r13751321 = r13751316 / r13751320;
double r13751322 = r13751319 * r13751321;
double r13751323 = r13751322 * r13751319;
double r13751324 = r13751318 - r13751323;
double r13751325 = r13751315 ? r13751324 : r13751313;
double r13751326 = r13751304 ? r13751313 : r13751325;
return r13751326;
}
\frac{x}{x + 1} - \frac{x + 1}{x - 1}\begin{array}{l}
\mathbf{if}\;x \le -0.9912853724188615:\\
\;\;\;\;\frac{-1}{x \cdot x} + \left(\frac{-3}{\left(x \cdot x\right) \cdot x} + \frac{-3}{x}\right)\\
\mathbf{elif}\;x \le 12935.179381933913:\\
\;\;\;\;\frac{x}{1 + x} - \left(\sqrt{1 + x} \cdot \frac{1}{x - 1}\right) \cdot \sqrt{1 + x}\\
\mathbf{else}:\\
\;\;\;\;\frac{-1}{x \cdot x} + \left(\frac{-3}{\left(x \cdot x\right) \cdot x} + \frac{-3}{x}\right)\\
\end{array}


Bits error versus x
if x < -0.9912853724188615 or 12935.179381933913 < x Initial program 59.0
Taylor expanded around -inf 0.5
Simplified0.5
Taylor expanded around -inf 0.5
Simplified0.2
if -0.9912853724188615 < x < 12935.179381933913Initial program 0.0
rmApplied div-inv0.0
rmApplied add-sqr-sqrt0.1
Applied associate-*l*0.1
Final simplification0.1
herbie shell --seed 2019101
(FPCore (x)
:name "Asymptote C"
(- (/ x (+ x 1)) (/ (+ x 1) (- x 1))))