\frac{\cos \left(2 \cdot x\right)}{{cos}^{2} \cdot \left(\left(x \cdot {sin}^{2}\right) \cdot x\right)}\begin{array}{l}
\mathbf{if}\;cos \le -4.65962850443893997 \cdot 10^{-117}:\\
\;\;\;\;\frac{\frac{\frac{\cos \left(2 \cdot x\right)}{\sqrt{\left|{cos}^{\left(\frac{2}{2}\right)} \cdot \left(x \cdot {sin}^{\left(\frac{2}{2}\right)}\right)\right|}}}{\sqrt{\left|{cos}^{\left(\frac{2}{2}\right)} \cdot \left(x \cdot {sin}^{\left(\frac{2}{2}\right)}\right)\right|}}}{\left|{cos}^{\left(\frac{2}{2}\right)} \cdot \left(x \cdot {sin}^{\left(\frac{2}{2}\right)}\right)\right|}\\
\mathbf{elif}\;cos \le 4.88319876901602 \cdot 10^{-111}:\\
\;\;\;\;\frac{\cos \left(2 \cdot x\right)}{{\left(\left|{\left({cos}^{1} \cdot {sin}^{1}\right)}^{1} \cdot x\right|\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\cos \left(2 \cdot x\right)}{\left|{cos}^{\left(\frac{\frac{2}{2}}{2}\right)} \cdot \left({cos}^{\left(\frac{\frac{2}{2}}{2}\right)} \cdot \left(x \cdot {sin}^{\left(\frac{2}{2}\right)}\right)\right)\right|}}{\left|{cos}^{\left(\frac{2}{2}\right)} \cdot \left(x \cdot {sin}^{\left(\frac{2}{2}\right)}\right)\right|}\\
\end{array}double f(double x, double cos, double sin) {
double r68200 = 2.0;
double r68201 = x;
double r68202 = r68200 * r68201;
double r68203 = cos(r68202);
double r68204 = cos;
double r68205 = pow(r68204, r68200);
double r68206 = sin;
double r68207 = pow(r68206, r68200);
double r68208 = r68201 * r68207;
double r68209 = r68208 * r68201;
double r68210 = r68205 * r68209;
double r68211 = r68203 / r68210;
return r68211;
}
double f(double x, double cos, double sin) {
double r68212 = cos;
double r68213 = -4.65962850443894e-117;
bool r68214 = r68212 <= r68213;
double r68215 = 2.0;
double r68216 = x;
double r68217 = r68215 * r68216;
double r68218 = cos(r68217);
double r68219 = 2.0;
double r68220 = r68215 / r68219;
double r68221 = pow(r68212, r68220);
double r68222 = sin;
double r68223 = pow(r68222, r68220);
double r68224 = r68216 * r68223;
double r68225 = r68221 * r68224;
double r68226 = fabs(r68225);
double r68227 = sqrt(r68226);
double r68228 = r68218 / r68227;
double r68229 = r68228 / r68227;
double r68230 = r68229 / r68226;
double r68231 = 4.88319876901602e-111;
bool r68232 = r68212 <= r68231;
double r68233 = 1.0;
double r68234 = pow(r68212, r68233);
double r68235 = pow(r68222, r68233);
double r68236 = r68234 * r68235;
double r68237 = pow(r68236, r68233);
double r68238 = r68237 * r68216;
double r68239 = fabs(r68238);
double r68240 = pow(r68239, r68219);
double r68241 = r68218 / r68240;
double r68242 = r68220 / r68219;
double r68243 = pow(r68212, r68242);
double r68244 = r68243 * r68224;
double r68245 = r68243 * r68244;
double r68246 = fabs(r68245);
double r68247 = r68218 / r68246;
double r68248 = r68247 / r68226;
double r68249 = r68232 ? r68241 : r68248;
double r68250 = r68214 ? r68230 : r68249;
return r68250;
}



Bits error versus x



Bits error versus cos



Bits error versus sin
Results
if cos < -4.65962850443894e-117Initial program 22.0
rmApplied sqr-pow22.0
Applied associate-*r*14.8
rmApplied add-sqr-sqrt14.9
Simplified14.8
Simplified1.8
rmApplied associate-/r*1.6
rmApplied add-sqr-sqrt1.6
Applied associate-/r*1.6
if -4.65962850443894e-117 < cos < 4.88319876901602e-111Initial program 52.6
rmApplied sqr-pow52.6
Applied associate-*r*50.4
rmApplied add-sqr-sqrt50.4
Simplified50.4
Simplified5.4
Taylor expanded around inf 4.8
if 4.88319876901602e-111 < cos Initial program 22.9
rmApplied sqr-pow22.9
Applied associate-*r*16.3
rmApplied add-sqr-sqrt16.4
Simplified16.3
Simplified2.0
rmApplied associate-/r*1.6
rmApplied sqr-pow1.7
Applied associate-*l*1.7
Final simplification2.3
herbie shell --seed 2020089 +o rules:numerics
(FPCore (x cos sin)
:name "cos(2*x)/(cos^2(x)*sin^2(x))"
:precision binary64
(/ (cos (* 2 x)) (* (pow cos 2) (* (* x (pow sin 2)) x))))