\frac{1}{x} - \frac{1}{\tan x}0.02222222222222222307030925492199457949027 \cdot {x}^{3} + \left(0.002116402116402116544841005563171165704262 \cdot {x}^{5} + 0.3333333333333333148296162562473909929395 \cdot x\right)double f(double x) {
double r120740 = 1.0;
double r120741 = x;
double r120742 = r120740 / r120741;
double r120743 = tan(r120741);
double r120744 = r120740 / r120743;
double r120745 = r120742 - r120744;
return r120745;
}
double f(double x) {
double r120746 = 0.022222222222222223;
double r120747 = x;
double r120748 = 3.0;
double r120749 = pow(r120747, r120748);
double r120750 = r120746 * r120749;
double r120751 = 0.0021164021164021165;
double r120752 = 5.0;
double r120753 = pow(r120747, r120752);
double r120754 = r120751 * r120753;
double r120755 = 0.3333333333333333;
double r120756 = r120755 * r120747;
double r120757 = r120754 + r120756;
double r120758 = r120750 + r120757;
return r120758;
}




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 2019297
(FPCore (x)
:name "invcot (example 3.9)"
:precision binary64
:pre (and (< -0.0259999999999999988 x) (< x 0.0259999999999999988))
:herbie-target
(if (< (fabs x) 0.0259999999999999988) (* (/ x 3) (+ 1 (/ (* x x) 15))) (- (/ 1 x) (/ 1 (tan x))))
(- (/ 1 x) (/ 1 (tan x))))