{\left(a \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2}{\left(a \cdot \left(\left(\sqrt[3]{\cos \left(\frac{\sqrt[3]{angle}}{\sqrt{180}} \cdot \left(\pi \cdot \left(\sqrt[3]{angle} \cdot \frac{\sqrt[3]{angle}}{\sqrt{180}}\right)\right)\right)} \cdot \sqrt[3]{\cos \left(\frac{\pi \cdot angle}{{\left(\sqrt{180}\right)}^{2}}\right)}\right) \cdot \sqrt[3]{\cos \left(\pi \cdot \frac{angle}{180}\right)}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2}(FPCore (a b angle) :precision binary64 (+ (pow (* a (cos (* PI (/ angle 180.0)))) 2.0) (pow (* b (sin (* PI (/ angle 180.0)))) 2.0)))
(FPCore (a b angle)
:precision binary64
(+
(pow
(*
a
(*
(*
(cbrt
(cos
(*
(/ (cbrt angle) (sqrt 180.0))
(* PI (* (cbrt angle) (/ (cbrt angle) (sqrt 180.0)))))))
(cbrt (cos (/ (* PI angle) (pow (sqrt 180.0) 2.0)))))
(cbrt (cos (* PI (/ angle 180.0))))))
2.0)
(pow (* b (sin (* PI (/ angle 180.0)))) 2.0)))double code(double a, double b, double angle) {
return pow((a * cos(((double) M_PI) * (angle / 180.0))), 2.0) + pow((b * sin(((double) M_PI) * (angle / 180.0))), 2.0);
}
double code(double a, double b, double angle) {
return pow((a * ((cbrt(cos((cbrt(angle) / sqrt(180.0)) * (((double) M_PI) * (cbrt(angle) * (cbrt(angle) / sqrt(180.0)))))) * cbrt(cos((((double) M_PI) * angle) / pow(sqrt(180.0), 2.0)))) * cbrt(cos(((double) M_PI) * (angle / 180.0))))), 2.0) + pow((b * sin(((double) M_PI) * (angle / 180.0))), 2.0);
}



Bits error versus a



Bits error versus b



Bits error versus angle
Results
Initial program 20.5
rmApplied add-sqr-sqrt_binary64_10020.6
Applied add-cube-cbrt_binary64_11320.6
Applied times-frac_binary64_8420.6
Applied associate-*r*_binary64_1820.6
Simplified20.6
rmApplied add-cube-cbrt_binary64_11320.6
Simplified20.6
Taylor expanded around inf 20.6
Final simplification20.6
herbie shell --seed 2021176
(FPCore (a b angle)
:name "ab-angle->ABCF C"
:precision binary64
(+ (pow (* a (cos (* PI (/ angle 180.0)))) 2.0) (pow (* b (sin (* PI (/ angle 180.0)))) 2.0)))