\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 r130464 = 1.0;
double r130465 = x;
double r130466 = r130464 / r130465;
double r130467 = tan(r130465);
double r130468 = r130464 / r130467;
double r130469 = r130466 - r130468;
return r130469;
}
double f(double x) {
double r130470 = 0.022222222222222223;
double r130471 = x;
double r130472 = 3.0;
double r130473 = pow(r130471, r130472);
double r130474 = r130470 * r130473;
double r130475 = 0.0021164021164021165;
double r130476 = 5.0;
double r130477 = pow(r130471, r130476);
double r130478 = r130475 * r130477;
double r130479 = 0.3333333333333333;
double r130480 = r130479 * r130471;
double r130481 = r130478 + r130480;
double r130482 = r130474 + r130481;
return r130482;
}




Bits error versus x
Results
| Original | 59.8 |
|---|---|
| Target | 0.1 |
| Herbie | 0.3 |
Initial program 59.8
Taylor expanded around 0 0.3
Final simplification0.3
herbie shell --seed 2020021
(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))))