\cos \left(x + \varepsilon\right) - \cos x
\begin{array}{l}
\mathbf{if}\;\varepsilon \le -3.02333196694931366 \cdot 10^{-13}:\\
\;\;\;\;\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}\\
\mathbf{elif}\;\varepsilon \le 1.16717710983658106 \cdot 10^{-10}:\\
\;\;\;\;\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 <= -3.0233319669493137e-13)) {
VAR = ((pow(((cos(x) * cos(eps)) - (sin(x) * sin(eps))), 3.0) - pow(cos(x), 3.0)) / ((((cos(eps) * cos(x)) - (sin(x) * sin(eps))) * ((pow(((cos(x) * cos(eps)) - (sin(x) * sin(eps))), 2.0) - pow(cos(x), 2.0)) / (((cos(x) * cos(eps)) - (sin(x) * sin(eps))) - cos(x)))) + (cos(x) * cos(x))));
} else {
double VAR_1;
if ((eps <= 1.167177109836581e-10)) {
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 < -3.0233319669493137e-13Initial program 31.7
rmApplied cos-sum1.8
rmApplied flip3--2.0
Simplified2.0
rmApplied flip-+2.0
Simplified2.0
if -3.0233319669493137e-13 < eps < 1.167177109836581e-10Initial program 48.4
Taylor expanded around 0 30.9
Simplified30.9
if 1.167177109836581e-10 < eps Initial program 30.7
rmApplied cos-sum1.5
Applied associate--l-1.5
rmApplied add-cbrt-cube1.7
Simplified1.7
Final simplification16.1
herbie shell --seed 2020091
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
:precision binary64
(- (cos (+ x eps)) (cos x)))