double code(double x, double eps) {
return ((double) (((double) cos(((double) (x + eps)))) - ((double) cos(x))));
}
double code(double x, double eps) {
return ((double) (-2.0 * ((double) (((double) (((double) sin(x)) * ((double) (((double) cos(((double) (eps * 0.5)))) * ((double) sin(((double) (eps / 2.0)))))))) + ((double) (((double) cos(x)) * ((double) (((double) sin(((double) (eps / 2.0)))) * ((double) sin(((double) (eps * 0.5))))))))))));
}



Bits error versus x



Bits error versus eps
Results
Initial program 39.1
rmApplied diff-cos33.6
Simplified15.0
Taylor expanded around inf 14.9
Simplified14.9
rmApplied sin-sum0.4
Applied distribute-lft-in0.4
Simplified0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2020191
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
:precision binary64
(- (cos (+ x eps)) (cos x)))