\frac{1}{x} - \frac{1}{\tan x}0.0222222222222222231 \cdot {x}^{3} + \left(0.00211640211640211654 \cdot {x}^{5} + 0.333333333333333315 \cdot x\right)double f(double x) {
double r118183 = 1.0;
double r118184 = x;
double r118185 = r118183 / r118184;
double r118186 = tan(r118184);
double r118187 = r118183 / r118186;
double r118188 = r118185 - r118187;
return r118188;
}
double f(double x) {
double r118189 = 0.022222222222222223;
double r118190 = x;
double r118191 = 3.0;
double r118192 = pow(r118190, r118191);
double r118193 = r118189 * r118192;
double r118194 = 0.0021164021164021165;
double r118195 = 5.0;
double r118196 = pow(r118190, r118195);
double r118197 = r118194 * r118196;
double r118198 = 0.3333333333333333;
double r118199 = r118198 * r118190;
double r118200 = r118197 + r118199;
double r118201 = r118193 + r118200;
return r118201;
}




Bits error versus x
Results
| Original | 59.9 |
|---|---|
| Target | 0.1 |
| Herbie | 0.3 |
Initial program 59.9
Taylor expanded around 0 0.3
Final simplification0.3
herbie shell --seed 2020065
(FPCore (x)
:name "invcot (example 3.9)"
:precision binary64
: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))))