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



Bits error versus x
Results
if x < -0.0302754788051601Initial program 1.0
rmApplied associate-/r*0.5
if -0.0302754788051601 < x < 0.031101198683564729Initial program 62.4
Taylor expanded around 0 0.0
Simplified0.0
if 0.031101198683564729 < x Initial program 1.1
rmApplied clear-num1.1
Simplified1.1
Final simplification0.4
herbie shell --seed 2020192
(FPCore (x)
:name "cos2 (problem 3.4.1)"
:precision binary64
(/ (- 1.0 (cos x)) (* x x)))