\cos \left(x + \varepsilon\right) - \cos x
\begin{array}{l}
\mathbf{if}\;\varepsilon \le -2.9500988603657079 \cdot 10^{-11}:\\
\;\;\;\;\log \left(e^{\cos x \cdot \cos \varepsilon - \left(\sin x \cdot \sin \varepsilon + \cos x\right)}\right)\\
\mathbf{elif}\;\varepsilon \le 4.8007364015852064 \cdot 10^{-11}:\\
\;\;\;\;\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 - \sqrt[3]{{\left(\sin x \cdot \sin \varepsilon + \cos x\right)}^{3}}\\
\end{array}double code(double x, double eps) {
return (cos((x + eps)) - cos(x));
}
double code(double x, double eps) {
double VAR;
if ((eps <= -2.950098860365708e-11)) {
VAR = log(exp(((cos(x) * cos(eps)) - ((sin(x) * sin(eps)) + cos(x)))));
} else {
double VAR_1;
if ((eps <= 4.8007364015852064e-11)) {
VAR_1 = (eps * (((0.16666666666666666 * pow(x, 3.0)) - x) - (eps * 0.5)));
} else {
VAR_1 = ((cos(x) * cos(eps)) - cbrt(pow(((sin(x) * sin(eps)) + cos(x)), 3.0)));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus x



Bits error versus eps
Results
if eps < -2.950098860365708e-11Initial program 30.6
rmApplied cos-sum1.5
Applied associate--l-1.6
rmApplied add-log-exp1.7
Applied add-log-exp1.7
Applied sum-log1.7
Applied add-log-exp1.9
Applied diff-log2.0
Simplified1.7
if -2.950098860365708e-11 < eps < 4.8007364015852064e-11Initial program 49.6
Taylor expanded around 0 31.1
Simplified31.1
if 4.8007364015852064e-11 < eps Initial program 29.9
rmApplied cos-sum1.5
Applied associate--l-1.5
rmApplied add-cbrt-cube1.7
Simplified1.7
Final simplification15.7
herbie shell --seed 2020092
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
:precision binary64
(- (cos (+ x eps)) (cos x)))