\cos \left(x + \varepsilon\right) - \cos x
\begin{array}{l}
\mathbf{if}\;\varepsilon \le -5.4894795427244131 \cdot 10^{-7}:\\
\;\;\;\;\cos x \cdot \cos \varepsilon - \log \left(e^{\sin x \cdot \sin \varepsilon + \cos x}\right)\\
\mathbf{elif}\;\varepsilon \le 3.3962311967673259 \cdot 10^{-13}:\\
\;\;\;\;\varepsilon \cdot \left(\left(\frac{1}{6} \cdot {x}^{3} - x\right) - \varepsilon \cdot \frac{1}{2}\right)\\
\mathbf{else}:\\
\;\;\;\;\cos x \cdot \cos \varepsilon - \left(\sqrt[3]{{\left(\sin x \cdot \sin \varepsilon\right)}^{3}} + \cos x\right)\\
\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 <= -5.489479542724413e-07)) {
VAR = ((double) (((double) (((double) cos(x)) * ((double) cos(eps)))) - ((double) log(((double) exp(((double) (((double) (((double) sin(x)) * ((double) sin(eps)))) + ((double) cos(x))))))))));
} else {
double VAR_1;
if ((eps <= 3.396231196767326e-13)) {
VAR_1 = ((double) (eps * ((double) (((double) (((double) (0.16666666666666666 * ((double) pow(x, 3.0)))) - x)) - ((double) (eps * 0.5))))));
} else {
VAR_1 = ((double) (((double) (((double) cos(x)) * ((double) cos(eps)))) - ((double) (((double) cbrt(((double) pow(((double) (((double) sin(x)) * ((double) sin(eps)))), 3.0)))) + ((double) cos(x))))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus x



Bits error versus eps
Results
if eps < -5.489479542724413e-07Initial program 29.3
rmApplied cos-sum1.0
Applied associate--l-1.0
rmApplied add-log-exp1.1
Applied add-log-exp1.2
Applied sum-log1.2
Simplified1.1
if -5.489479542724413e-07 < eps < 3.396231196767326e-13Initial program 49.4
Taylor expanded around 0 31.5
Simplified31.5
if 3.396231196767326e-13 < eps Initial program 31.7
rmApplied cos-sum2.2
Applied associate--l-2.2
rmApplied add-cbrt-cube2.3
Applied add-cbrt-cube2.3
Applied cbrt-unprod2.3
Simplified2.3
Final simplification15.9
herbie shell --seed 2020122
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
:precision binary64
(- (cos (+ x eps)) (cos x)))