double code(double x) {
return (x - sin(x)) / (x - tan(x));
}
double code(double x) {
double tmp;
if ((x <= -0.026395687002885145) || !(x <= 0.027677485715175896)) {
tmp = 1.0 / ((x - tan(x)) / (x - sin(x)));
} else {
tmp = ((x * x) * 0.225) - (0.5 + (0.009642857142857142 * pow(x, 4.0)));
}
return tmp;
}



Bits error versus x
Results
if x < -0.026395687002885145 or 0.0276774857151758959 < x Initial program 0.0
rmApplied clear-num0.1
if -0.026395687002885145 < x < 0.0276774857151758959Initial program 63.2
Taylor expanded around 0 0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2020338
(FPCore (x)
:name "sintan (problem 3.4.5)"
:precision binary64
(/ (- x (sin x)) (- x (tan x))))