\frac{1}{x} - \frac{1}{\tan x}{x}^{5} \cdot \frac{2}{945} + \left(\left(\left(\frac{1}{45} \cdot x\right) \cdot x + \frac{1}{3}\right) \cdot \left(\left(\frac{1}{45} \cdot x\right) \cdot x - \frac{1}{3}\right)\right) \cdot \frac{x}{\left(\frac{1}{45} \cdot x\right) \cdot x - \frac{1}{3}}double f(double x) {
double r2242151 = 1.0;
double r2242152 = x;
double r2242153 = r2242151 / r2242152;
double r2242154 = tan(r2242152);
double r2242155 = r2242151 / r2242154;
double r2242156 = r2242153 - r2242155;
return r2242156;
}
double f(double x) {
double r2242157 = x;
double r2242158 = 5.0;
double r2242159 = pow(r2242157, r2242158);
double r2242160 = 0.0021164021164021165;
double r2242161 = r2242159 * r2242160;
double r2242162 = 0.022222222222222223;
double r2242163 = r2242162 * r2242157;
double r2242164 = r2242163 * r2242157;
double r2242165 = 0.3333333333333333;
double r2242166 = r2242164 + r2242165;
double r2242167 = r2242164 - r2242165;
double r2242168 = r2242166 * r2242167;
double r2242169 = r2242157 / r2242167;
double r2242170 = r2242168 * r2242169;
double r2242171 = r2242161 + r2242170;
return r2242171;
}




Bits error versus x
Results
| Original | 60.0 |
|---|---|
| Target | 0.1 |
| Herbie | 0.2 |
Initial program 60.0
Taylor expanded around 0 0.3
Simplified0.3
rmApplied flip-+0.3
Applied associate-*l/0.3
rmApplied *-un-lft-identity0.3
Applied times-frac0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2019134
(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))))