\cos \left(x + \varepsilon\right) - \cos x
\begin{array}{l}
\mathbf{if}\;\varepsilon \le -0.008151886891035409199446348793571814894676:\\
\;\;\;\;\cos x \cdot \cos \varepsilon - \left(\sin x \cdot \sin \varepsilon + \cos x\right)\\
\mathbf{elif}\;\varepsilon \le 0.04140103919524167758181576459719508420676:\\
\;\;\;\;\left(-2 \cdot \sin \left(\frac{\varepsilon}{2}\right)\right) \cdot \sin \left(\frac{\left(x + \varepsilon\right) + x}{2}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right) - \cos x\\
\end{array}double f(double x, double eps) {
double r30979 = x;
double r30980 = eps;
double r30981 = r30979 + r30980;
double r30982 = cos(r30981);
double r30983 = cos(r30979);
double r30984 = r30982 - r30983;
return r30984;
}
double f(double x, double eps) {
double r30985 = eps;
double r30986 = -0.00815188689103541;
bool r30987 = r30985 <= r30986;
double r30988 = x;
double r30989 = cos(r30988);
double r30990 = cos(r30985);
double r30991 = r30989 * r30990;
double r30992 = sin(r30988);
double r30993 = sin(r30985);
double r30994 = r30992 * r30993;
double r30995 = r30994 + r30989;
double r30996 = r30991 - r30995;
double r30997 = 0.04140103919524168;
bool r30998 = r30985 <= r30997;
double r30999 = -2.0;
double r31000 = 2.0;
double r31001 = r30985 / r31000;
double r31002 = sin(r31001);
double r31003 = r30999 * r31002;
double r31004 = r30988 + r30985;
double r31005 = r31004 + r30988;
double r31006 = r31005 / r31000;
double r31007 = sin(r31006);
double r31008 = r31003 * r31007;
double r31009 = r30991 - r30994;
double r31010 = r31009 - r30989;
double r31011 = r30998 ? r31008 : r31010;
double r31012 = r30987 ? r30996 : r31011;
return r31012;
}



Bits error versus x



Bits error versus eps
Results
if eps < -0.00815188689103541Initial program 30.7
rmApplied cos-sum0.8
Applied associate--l-0.8
if -0.00815188689103541 < eps < 0.04140103919524168Initial program 48.8
rmApplied diff-cos37.0
Simplified0.6
rmApplied associate-*r*0.6
if 0.04140103919524168 < eps Initial program 30.9
rmApplied cos-sum0.8
Final simplification0.7
herbie shell --seed 2019326
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
:precision binary64
(- (cos (+ x eps)) (cos x)))