{\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 \sqrt[3]{{\cos \left(\frac{\pi \cdot angle}{180}\right)}^{3}}\right)}^{2} + {\left(b \cdot \sin \left(\frac{\pi \cdot 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 (pow (cos (/ (* PI angle) 180.0)) 3.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(pow(cos((((double) M_PI) * angle) / 180.0), 3.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.3
rmApplied associate-*r/_binary64_2020.4
rmApplied add-cbrt-cube_binary64_11420.4
Simplified20.4
Final simplification20.4
herbie shell --seed 2020355
(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)))