\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \frac{1}{x - 1}\begin{array}{l}
\mathbf{if}\;x \le -92.034279482346449 \lor \neg \left(x \le 95.3045887785727217\right):\\
\;\;\;\;2 \cdot \left(\frac{1}{{x}^{7}} + \left(\frac{1}{{x}^{5}} + \frac{1}{{x}^{3}}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \left(\sqrt[3]{\frac{1}{x - 1}} \cdot \sqrt[3]{\frac{1}{x - 1}}\right) \cdot \sqrt[3]{\frac{1}{x - 1}}\\
\end{array}double code(double x) {
return ((double) (((double) (((double) (1.0 / ((double) (x + 1.0)))) - ((double) (2.0 / x)))) + ((double) (1.0 / ((double) (x - 1.0))))));
}
double code(double x) {
double VAR;
if (((x <= -92.03427948234645) || !(x <= 95.30458877857272))) {
VAR = ((double) (2.0 * ((double) (((double) (1.0 / ((double) pow(x, 7.0)))) + ((double) (((double) (1.0 / ((double) pow(x, 5.0)))) + ((double) (1.0 / ((double) pow(x, 3.0))))))))));
} else {
VAR = ((double) (((double) (((double) (1.0 / ((double) (x + 1.0)))) - ((double) (2.0 / x)))) + ((double) (((double) (((double) cbrt(((double) (1.0 / ((double) (x - 1.0)))))) * ((double) cbrt(((double) (1.0 / ((double) (x - 1.0)))))))) * ((double) cbrt(((double) (1.0 / ((double) (x - 1.0))))))))));
}
return VAR;
}




Bits error versus x
Results
| Original | 9.7 |
|---|---|
| Target | 0.2 |
| Herbie | 0.3 |
if x < -92.034279482346449 or 95.3045887785727217 < x Initial program 19.7
Taylor expanded around inf 0.5
Simplified0.5
if -92.034279482346449 < x < 95.3045887785727217Initial program 0.0
rmApplied add-cube-cbrt0.1
Final simplification0.3
herbie shell --seed 2020161
(FPCore (x)
:name "3frac (problem 3.3.3)"
:precision binary64
:herbie-target
(/ 2.0 (* x (- (* x x) 1.0)))
(+ (- (/ 1.0 (+ x 1.0)) (/ 2.0 x)) (/ 1.0 (- x 1.0))))