{\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(\sqrt[3]{\cos \left(\left(\pi \cdot \left(\sqrt[3]{angle \cdot 0.005555555555555556} \cdot \sqrt[3]{angle \cdot 0.005555555555555556}\right)\right) \cdot \left({\left(\sqrt{180}\right)}^{-0.3333333333333333} \cdot \sqrt[3]{\frac{angle}{\sqrt{180}}}\right)\right)} \cdot \left(\sqrt[3]{\cos \left(\left(\pi \cdot \left(\sqrt[3]{angle \cdot 0.005555555555555556} \cdot \sqrt[3]{angle \cdot 0.005555555555555556}\right)\right) \cdot \left({\left(\sqrt{180}\right)}^{-0.3333333333333333} \cdot \sqrt[3]{\frac{angle}{\sqrt{180}}}\right)\right)} \cdot \sqrt[3]{\cos \left(\left(\pi \cdot \left(\sqrt[3]{angle \cdot 0.005555555555555556} \cdot \sqrt[3]{angle \cdot 0.005555555555555556}\right)\right) \cdot \left({\left(\sqrt{180}\right)}^{-0.3333333333333333} \cdot \sqrt[3]{\frac{angle}{\sqrt{180}}}\right)\right)}\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(0.005555555555555556 \cdot \left(\pi \cdot angle\right)\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
(*
(*
PI
(*
(cbrt (* angle 0.005555555555555556))
(cbrt (* angle 0.005555555555555556))))
(*
(pow (sqrt 180.0) -0.3333333333333333)
(cbrt (/ angle (sqrt 180.0)))))))
(*
(cbrt
(cos
(*
(*
PI
(*
(cbrt (* angle 0.005555555555555556))
(cbrt (* angle 0.005555555555555556))))
(*
(pow (sqrt 180.0) -0.3333333333333333)
(cbrt (/ angle (sqrt 180.0)))))))
(cbrt
(cos
(*
(*
PI
(*
(cbrt (* angle 0.005555555555555556))
(cbrt (* angle 0.005555555555555556))))
(*
(pow (sqrt 180.0) -0.3333333333333333)
(cbrt (/ angle (sqrt 180.0))))))))))
2.0)
(pow (* b (sin (* 0.005555555555555556 (* PI angle)))) 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((((double) M_PI) * (cbrt(angle * 0.005555555555555556) * cbrt(angle * 0.005555555555555556))) * (pow(sqrt(180.0), -0.3333333333333333) * cbrt(angle / sqrt(180.0))))) * (cbrt(cos((((double) M_PI) * (cbrt(angle * 0.005555555555555556) * cbrt(angle * 0.005555555555555556))) * (pow(sqrt(180.0), -0.3333333333333333) * cbrt(angle / sqrt(180.0))))) * cbrt(cos((((double) M_PI) * (cbrt(angle * 0.005555555555555556) * cbrt(angle * 0.005555555555555556))) * (pow(sqrt(180.0), -0.3333333333333333) * cbrt(angle / sqrt(180.0)))))))), 2.0) + pow((b * sin(0.005555555555555556 * (((double) M_PI) * angle))), 2.0);
}



Bits error versus a



Bits error versus b



Bits error versus angle
Results
Initial program 20.2
Taylor expanded around inf 20.2
Simplified20.2
rmApplied add-cube-cbrt_binary64_11320.2
Applied associate-*r*_binary64_1820.2
Simplified20.2
rmApplied add-sqr-sqrt_binary64_10020.2
Applied *-un-lft-identity_binary64_7820.2
Applied times-frac_binary64_8420.2
Applied cbrt-prod_binary64_10920.2
rmApplied add-cube-cbrt_binary64_11320.2
Simplified20.2
Simplified20.2
Final simplification20.2
herbie shell --seed 2020346
(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)))