\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}\;x \le -3.899795259514831596491496327990136660459 \cdot 10^{-188}:\\
\;\;\;\;\frac{\frac{\cos \left(2 \cdot x\right)}{\left|{\left({cos}^{1} \cdot {sin}^{1}\right)}^{1} \cdot x\right|}}{{\left(\sqrt{\left|{\left({cos}^{1} \cdot {sin}^{1}\right)}^{1} \cdot x\right|}\right)}^{2}}\\
\mathbf{elif}\;x \le 3.099208061336091548327514698736633602822 \cdot 10^{-150}:\\
\;\;\;\;\frac{1}{\left|{cos}^{\left(\frac{2}{2}\right)} \cdot \left(x \cdot {sin}^{\left(\frac{2}{2}\right)}\right)\right|} \cdot \frac{\cos \left(2 \cdot x\right)}{\left|{cos}^{\left(\frac{2}{2}\right)} \cdot \left(x \cdot {sin}^{\left(\frac{2}{2}\right)}\right)\right|}\\
\mathbf{else}:\\
\;\;\;\;\cos \left(2 \cdot x\right) \cdot \frac{1}{{\left(\left|{\left({cos}^{1} \cdot {sin}^{1}\right)}^{1} \cdot x\right|\right)}^{2}}\\
\end{array}double f(double x, double cos, double sin) {
double r52875 = 2.0;
double r52876 = x;
double r52877 = r52875 * r52876;
double r52878 = cos(r52877);
double r52879 = cos;
double r52880 = pow(r52879, r52875);
double r52881 = sin;
double r52882 = pow(r52881, r52875);
double r52883 = r52876 * r52882;
double r52884 = r52883 * r52876;
double r52885 = r52880 * r52884;
double r52886 = r52878 / r52885;
return r52886;
}
double f(double x, double cos, double sin) {
double r52887 = x;
double r52888 = -3.8997952595148316e-188;
bool r52889 = r52887 <= r52888;
double r52890 = 2.0;
double r52891 = r52890 * r52887;
double r52892 = cos(r52891);
double r52893 = cos;
double r52894 = 1.0;
double r52895 = pow(r52893, r52894);
double r52896 = sin;
double r52897 = pow(r52896, r52894);
double r52898 = r52895 * r52897;
double r52899 = pow(r52898, r52894);
double r52900 = r52899 * r52887;
double r52901 = fabs(r52900);
double r52902 = r52892 / r52901;
double r52903 = sqrt(r52901);
double r52904 = 2.0;
double r52905 = pow(r52903, r52904);
double r52906 = r52902 / r52905;
double r52907 = 3.0992080613360915e-150;
bool r52908 = r52887 <= r52907;
double r52909 = 1.0;
double r52910 = r52890 / r52904;
double r52911 = pow(r52893, r52910);
double r52912 = pow(r52896, r52910);
double r52913 = r52887 * r52912;
double r52914 = r52911 * r52913;
double r52915 = fabs(r52914);
double r52916 = r52909 / r52915;
double r52917 = r52892 / r52915;
double r52918 = r52916 * r52917;
double r52919 = pow(r52901, r52904);
double r52920 = r52909 / r52919;
double r52921 = r52892 * r52920;
double r52922 = r52908 ? r52918 : r52921;
double r52923 = r52889 ? r52906 : r52922;
return r52923;
}



Bits error versus x



Bits error versus cos



Bits error versus sin
Results
if x < -3.8997952595148316e-188Initial program 26.2
rmApplied sqr-pow26.2
Applied associate-*r*21.0
rmApplied add-sqr-sqrt21.1
Simplified21.0
Simplified2.6
Taylor expanded around inf 2.1
rmApplied add-sqr-sqrt2.3
Applied unpow-prod-down2.3
Applied associate-/r*1.9
Simplified1.8
if -3.8997952595148316e-188 < x < 3.0992080613360915e-150Initial program 44.2
rmApplied sqr-pow44.2
Applied associate-*r*30.4
rmApplied add-sqr-sqrt30.5
Simplified30.4
Simplified4.8
rmApplied *-un-lft-identity4.8
Applied times-frac4.6
if 3.0992080613360915e-150 < x Initial program 25.7
rmApplied sqr-pow25.7
Applied associate-*r*20.3
rmApplied add-sqr-sqrt20.4
Simplified20.4
Simplified2.4
Taylor expanded around inf 1.8
rmApplied div-inv1.8
Final simplification2.2
herbie shell --seed 2019352
(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))))