{\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 \cos \left(\left(\pi \cdot \frac{\sqrt[3]{angle} \cdot \sqrt[3]{angle}}{\sqrt{180}}\right) \cdot \frac{\sqrt[3]{\sqrt[3]{angle}} \cdot \left(\sqrt[3]{\sqrt[3]{angle}} \cdot \sqrt[3]{\sqrt[3]{angle}}\right)}{\sqrt{180}}\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
(cos
(*
(* PI (/ (* (cbrt angle) (cbrt angle)) (sqrt 180.0)))
(/
(* (cbrt (cbrt angle)) (* (cbrt (cbrt angle)) (cbrt (cbrt angle))))
(sqrt 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 * cos((((double) M_PI) * ((cbrt(angle) * cbrt(angle)) / sqrt(180.0))) * ((cbrt(cbrt(angle)) * (cbrt(cbrt(angle)) * cbrt(cbrt(angle)))) / sqrt(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.4
rmApplied add-sqr-sqrt_binary64_10020.4
Applied add-cube-cbrt_binary64_11320.4
Applied times-frac_binary64_8420.4
Applied associate-*r*_binary64_1820.4
rmApplied add-cube-cbrt_binary64_11320.4
Final simplification20.4
herbie shell --seed 2020343
(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)))