\cos \left(x + \varepsilon\right) - \cos x
\begin{array}{l}
\mathbf{if}\;\varepsilon \le -1.184078211908822 \cdot 10^{-13}:\\
\;\;\;\;\log \left(e^{\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right) - \cos x}\right)\\
\mathbf{elif}\;\varepsilon \le 1.69858524052719363 \cdot 10^{-7}:\\
\;\;\;\;\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 - \frac{\left(\sin x \cdot \sin \varepsilon\right) \cdot \left(\sin x \cdot \sin \varepsilon\right) - \cos x \cdot \cos x}{\sin x \cdot \sin \varepsilon - \cos x}\\
\end{array}double code(double x, double eps) {
return (cos((x + eps)) - cos(x));
}
double code(double x, double eps) {
double VAR;
if ((eps <= -1.184078211908822e-13)) {
VAR = log(exp((((cos(x) * cos(eps)) - (sin(x) * sin(eps))) - cos(x))));
} else {
double VAR_1;
if ((eps <= 1.6985852405271936e-07)) {
VAR_1 = (eps * (((0.16666666666666666 * pow(x, 3.0)) - x) - (eps * 0.5)));
} else {
VAR_1 = ((cos(x) * cos(eps)) - ((((sin(x) * sin(eps)) * (sin(x) * sin(eps))) - (cos(x) * cos(x))) / ((sin(x) * sin(eps)) - cos(x))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus x



Bits error versus eps
Results
if eps < -1.184078211908822e-13Initial program 30.9
rmApplied cos-sum1.9
rmApplied add-log-exp2.0
Applied add-log-exp2.0
Applied add-log-exp2.2
Applied diff-log2.2
Applied diff-log2.4
Simplified2.0
if -1.184078211908822e-13 < eps < 1.6985852405271936e-07Initial program 48.8
Taylor expanded around 0 31.7
Simplified31.7
if 1.6985852405271936e-07 < eps Initial program 29.8
rmApplied cos-sum1.1
Applied associate--l-1.1
rmApplied flip-+1.2
Final simplification16.3
herbie shell --seed 2020079
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
:precision binary64
(- (cos (+ x eps)) (cos x)))