\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \frac{1}{x - 1}\frac{2}{\left(\left(x - 1\right) \cdot \left(x + 1\right)\right) \cdot x}double f(double x) {
double r4175424 = 1.0;
double r4175425 = x;
double r4175426 = r4175425 + r4175424;
double r4175427 = r4175424 / r4175426;
double r4175428 = 2.0;
double r4175429 = r4175428 / r4175425;
double r4175430 = r4175427 - r4175429;
double r4175431 = r4175425 - r4175424;
double r4175432 = r4175424 / r4175431;
double r4175433 = r4175430 + r4175432;
return r4175433;
}
double f(double x) {
double r4175434 = 2.0;
double r4175435 = x;
double r4175436 = 1.0;
double r4175437 = r4175435 - r4175436;
double r4175438 = r4175435 + r4175436;
double r4175439 = r4175437 * r4175438;
double r4175440 = r4175439 * r4175435;
double r4175441 = r4175434 / r4175440;
return r4175441;
}




Bits error versus x
Results
| Original | 9.9 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 9.9
rmApplied +-commutative9.9
rmApplied associate-+r-10.0
rmApplied frac-add26.3
Applied frac-sub26.0
Taylor expanded around 0 0.2
Final simplification0.2
herbie shell --seed 2019170 +o rules:numerics
(FPCore (x)
:name "3frac (problem 3.3.3)"
:herbie-target
(/ 2.0 (* x (- (* x x) 1.0)))
(+ (- (/ 1.0 (+ x 1.0)) (/ 2.0 x)) (/ 1.0 (- x 1.0))))