\cos \left(x + \varepsilon\right) - \cos x
\begin{array}{l}
\mathbf{if}\;\varepsilon \le -8.7495484524256561 \cdot 10^{-5}:\\
\;\;\;\;\frac{\left(\cos x \cdot \cos \varepsilon\right) \cdot \left(\cos x \cdot \cos \varepsilon\right) - \left(\sin x \cdot \sin \varepsilon\right) \cdot \left(\sin x \cdot \sin \varepsilon\right)}{\cos x \cdot \cos \varepsilon + \sin x \cdot \sin \varepsilon} - \cos x\\
\mathbf{elif}\;\varepsilon \le 1.87956082054209138 \cdot 10^{-15}:\\
\;\;\;\;\frac{1}{24} \cdot {\varepsilon}^{4} - \left(x \cdot \varepsilon + \frac{1}{2} \cdot {\varepsilon}^{2}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{{\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right)}^{3} - {\left(\cos x\right)}^{3}}{\left(\cos \varepsilon \cdot \cos x - \sin x \cdot \sin \varepsilon\right) \cdot \frac{{\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right)}^{2} - {\left(\cos x\right)}^{2}}{\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right) - \cos x} + \cos x \cdot \cos x}\\
\end{array}double code(double x, double eps) {
return ((double) (((double) cos(((double) (x + eps)))) - ((double) cos(x))));
}
double code(double x, double eps) {
double VAR;
if ((eps <= -8.749548452425656e-05)) {
VAR = ((double) (((double) (((double) (((double) (((double) (((double) cos(x)) * ((double) cos(eps)))) * ((double) (((double) cos(x)) * ((double) cos(eps)))))) - ((double) (((double) (((double) sin(x)) * ((double) sin(eps)))) * ((double) (((double) sin(x)) * ((double) sin(eps)))))))) / ((double) (((double) (((double) cos(x)) * ((double) cos(eps)))) + ((double) (((double) sin(x)) * ((double) sin(eps)))))))) - ((double) cos(x))));
} else {
double VAR_1;
if ((eps <= 1.8795608205420914e-15)) {
VAR_1 = ((double) (((double) (0.041666666666666664 * ((double) pow(eps, 4.0)))) - ((double) (((double) (x * eps)) + ((double) (0.5 * ((double) pow(eps, 2.0))))))));
} else {
VAR_1 = ((double) (((double) (((double) pow(((double) (((double) (((double) cos(x)) * ((double) cos(eps)))) - ((double) (((double) sin(x)) * ((double) sin(eps)))))), 3.0)) - ((double) pow(((double) cos(x)), 3.0)))) / ((double) (((double) (((double) (((double) (((double) cos(eps)) * ((double) cos(x)))) - ((double) (((double) sin(x)) * ((double) sin(eps)))))) * ((double) (((double) (((double) pow(((double) (((double) (((double) cos(x)) * ((double) cos(eps)))) - ((double) (((double) sin(x)) * ((double) sin(eps)))))), 2.0)) - ((double) pow(((double) cos(x)), 2.0)))) / ((double) (((double) (((double) (((double) cos(x)) * ((double) cos(eps)))) - ((double) (((double) sin(x)) * ((double) sin(eps)))))) - ((double) cos(x)))))))) + ((double) (((double) cos(x)) * ((double) cos(x))))))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus x



Bits error versus eps
Results
if eps < -8.7495484524256561e-5Initial program 30.8
rmApplied cos-sum0.9
rmApplied flip--1.0
if -8.7495484524256561e-5 < eps < 1.87956082054209138e-15Initial program 49.6
rmApplied cos-sum49.3
rmApplied add-cbrt-cube49.3
Simplified49.3
Taylor expanded around 0 30.9
if 1.87956082054209138e-15 < eps Initial program 30.8
rmApplied cos-sum2.3
rmApplied flip3--2.4
Simplified2.4
rmApplied flip-+2.4
Simplified2.4
Final simplification16.0
herbie shell --seed 2020157
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
:precision binary64
(- (cos (+ x eps)) (cos x)))