{\left(a \cdot \sin \left(\frac{angle}{180} \cdot \pi\right)\right)}^{2} + {\left(b \cdot \cos \left(\frac{angle}{180} \cdot \pi\right)\right)}^{2}{\left(a \cdot \sin \left(\frac{1}{\sqrt{180}} \cdot \left(\pi \cdot \frac{angle}{\sqrt{180}}\right)\right)\right)}^{2} + {\left(b \cdot \sqrt[3]{{\left(\sqrt[3]{{\cos \left(\pi \cdot \frac{angle}{180}\right)}^{3}}\right)}^{3}}\right)}^{2}(FPCore (a b angle) :precision binary64 (+ (pow (* a (sin (* (/ angle 180.0) PI))) 2.0) (pow (* b (cos (* (/ angle 180.0) PI))) 2.0)))
(FPCore (a b angle) :precision binary64 (+ (pow (* a (sin (* (/ 1.0 (sqrt 180.0)) (* PI (/ angle (sqrt 180.0)))))) 2.0) (pow (* b (cbrt (pow (cbrt (pow (cos (* PI (/ angle 180.0))) 3.0)) 3.0))) 2.0)))
double code(double a, double b, double angle) {
return pow((a * sin((angle / 180.0) * ((double) M_PI))), 2.0) + pow((b * cos((angle / 180.0) * ((double) M_PI))), 2.0);
}
double code(double a, double b, double angle) {
return pow((a * sin((1.0 / sqrt(180.0)) * (((double) M_PI) * (angle / sqrt(180.0))))), 2.0) + pow((b * cbrt(pow(cbrt(pow(cos(((double) M_PI) * (angle / 180.0)), 3.0)), 3.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_44120.7
Applied *-un-lft-identity_binary64_41920.7
Applied times-frac_binary64_42520.5
Applied associate-*l*_binary64_36020.5
Simplified20.5
rmApplied add-cbrt-cube_binary64_45520.5
Simplified20.5
rmApplied add-cbrt-cube_binary64_45520.5
Simplified20.5
Final simplification20.5
herbie shell --seed 2020349
(FPCore (a b angle)
:name "ab-angle->ABCF A"
:precision binary64
(+ (pow (* a (sin (* (/ angle 180.0) PI))) 2.0) (pow (* b (cos (* (/ angle 180.0) PI))) 2.0)))