double code(double x) {
return ((double) (((double) (1.0 - ((double) cos(x)))) / ((double) (x * x))));
}
double code(double x) {
double VAR;
if ((x <= -0.033404312087007934)) {
VAR = ((double) (((double) (1.0 / x)) * ((double) (((double) (((double) (1.0 * 1.0)) - ((double) (((double) cos(x)) * ((double) cos(x)))))) / ((double) (x * ((double) (1.0 + ((double) cos(x))))))))));
} else {
double VAR_1;
if ((x <= 0.028650006657888936)) {
VAR_1 = ((double) (((double) (((double) pow(x, 4.0)) * 0.001388888888888889)) + ((double) (0.5 + ((double) (x * ((double) (x * -0.041666666666666664))))))));
} else {
VAR_1 = ((double) (((double) (1.0 / x)) * ((double) (((double) (((double) pow(1.0, 3.0)) - ((double) pow(((double) cos(x)), 3.0)))) / ((double) (x * ((double) (((double) (1.0 * 1.0)) + ((double) (((double) cos(x)) * ((double) (1.0 + ((double) cos(x))))))))))))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus x
Results
if x < -0.033404312087007934Initial program 1.1
rmApplied *-un-lft-identity1.1
Applied times-frac0.5
rmApplied flip--0.7
Applied associate-/l/0.7
if -0.033404312087007934 < x < 0.028650006657888936Initial program 62.2
Taylor expanded around 0 0.0
Simplified0.0
if 0.028650006657888936 < x Initial program 1.1
rmApplied *-un-lft-identity1.1
Applied times-frac0.5
rmApplied flip3--0.6
Applied associate-/l/0.6
Simplified0.6
Final simplification0.3
herbie shell --seed 2020190
(FPCore (x)
:name "cos2 (problem 3.4.1)"
:precision binary64
(/ (- 1.0 (cos x)) (* x x)))