\cos \left(x + \varepsilon\right) - \cos x
\begin{array}{l}
\mathbf{if}\;\varepsilon \le -1.153154815626698153605163400370998364602 \cdot 10^{-10} \lor \neg \left(\varepsilon \le 1.131984740243919275020376483120720934039 \cdot 10^{-6}\right):\\
\;\;\;\;\frac{{\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right)}^{3} - {\left(\cos x\right)}^{2} \cdot \cos x}{\left(\cos \varepsilon \cdot \cos x - \sin x \cdot \sin \varepsilon\right) \cdot \left(\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right) + \cos x\right) + \cos x \cdot \cos x}\\
\mathbf{else}:\\
\;\;\;\;\varepsilon \cdot \left(\left(\frac{1}{6} \cdot {x}^{3} - x\right) - \varepsilon \cdot \frac{1}{2}\right)\\
\end{array}double f(double x, double eps) {
double r67374 = x;
double r67375 = eps;
double r67376 = r67374 + r67375;
double r67377 = cos(r67376);
double r67378 = cos(r67374);
double r67379 = r67377 - r67378;
return r67379;
}
double f(double x, double eps) {
double r67380 = eps;
double r67381 = -1.1531548156266982e-10;
bool r67382 = r67380 <= r67381;
double r67383 = 1.1319847402439193e-06;
bool r67384 = r67380 <= r67383;
double r67385 = !r67384;
bool r67386 = r67382 || r67385;
double r67387 = x;
double r67388 = cos(r67387);
double r67389 = cos(r67380);
double r67390 = r67388 * r67389;
double r67391 = sin(r67387);
double r67392 = sin(r67380);
double r67393 = r67391 * r67392;
double r67394 = r67390 - r67393;
double r67395 = 3.0;
double r67396 = pow(r67394, r67395);
double r67397 = 2.0;
double r67398 = pow(r67388, r67397);
double r67399 = r67398 * r67388;
double r67400 = r67396 - r67399;
double r67401 = r67389 * r67388;
double r67402 = r67401 - r67393;
double r67403 = r67394 + r67388;
double r67404 = r67402 * r67403;
double r67405 = r67388 * r67388;
double r67406 = r67404 + r67405;
double r67407 = r67400 / r67406;
double r67408 = 0.16666666666666666;
double r67409 = pow(r67387, r67395);
double r67410 = r67408 * r67409;
double r67411 = r67410 - r67387;
double r67412 = 0.5;
double r67413 = r67380 * r67412;
double r67414 = r67411 - r67413;
double r67415 = r67380 * r67414;
double r67416 = r67386 ? r67407 : r67415;
return r67416;
}



Bits error versus x



Bits error versus eps
Results
if eps < -1.1531548156266982e-10 or 1.1319847402439193e-06 < eps Initial program 30.5
rmApplied cos-sum1.2
rmApplied flip3--1.4
Simplified1.4
rmApplied add-cube-cbrt2.0
Applied unpow-prod-down2.0
Simplified1.6
Simplified1.4
if -1.1531548156266982e-10 < eps < 1.1319847402439193e-06Initial program 49.7
Taylor expanded around 0 32.7
Simplified32.7
Final simplification16.8
herbie shell --seed 2019353
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
:precision binary64
(- (cos (+ x eps)) (cos x)))