\frac{1}{x} - \frac{1}{\tan x}\frac{x}{\frac{\frac{1}{3} - \frac{1}{45} \cdot \left(x \cdot x\right)}{\frac{1}{9} - \left(\frac{1}{45} \cdot \left(x \cdot x\right)\right) \cdot \left(\frac{1}{45} \cdot \left(x \cdot x\right)\right)}} + \frac{2}{945} \cdot {x}^{5}double f(double x) {
double r3956287 = 1.0;
double r3956288 = x;
double r3956289 = r3956287 / r3956288;
double r3956290 = tan(r3956288);
double r3956291 = r3956287 / r3956290;
double r3956292 = r3956289 - r3956291;
return r3956292;
}
double f(double x) {
double r3956293 = x;
double r3956294 = 0.3333333333333333;
double r3956295 = 0.022222222222222223;
double r3956296 = r3956293 * r3956293;
double r3956297 = r3956295 * r3956296;
double r3956298 = r3956294 - r3956297;
double r3956299 = 0.1111111111111111;
double r3956300 = r3956297 * r3956297;
double r3956301 = r3956299 - r3956300;
double r3956302 = r3956298 / r3956301;
double r3956303 = r3956293 / r3956302;
double r3956304 = 0.0021164021164021165;
double r3956305 = 5.0;
double r3956306 = pow(r3956293, r3956305);
double r3956307 = r3956304 * r3956306;
double r3956308 = r3956303 + r3956307;
return r3956308;
}




Bits error versus x
Results
| Original | 59.8 |
|---|---|
| Target | 0.1 |
| Herbie | 0.0 |
Initial program 59.8
Taylor expanded around 0 0.3
Simplified0.3
rmApplied flip-+0.3
Applied associate-*r/0.3
rmApplied associate-/l*0.0
Final simplification0.0
herbie shell --seed 2019112
(FPCore (x)
:name "invcot (example 3.9)"
:pre (and (< -0.026 x) (< x 0.026))
:herbie-target
(if (< (fabs x) 0.026) (* (/ x 3) (+ 1 (/ (* x x) 15))) (- (/ 1 x) (/ 1 (tan x))))
(- (/ 1 x) (/ 1 (tan x))))