r \cdot \frac{\sin b}{\cos \left(a + b\right)}r \cdot \left(\frac{\sin b}{\left(\cos a \cdot \cos b\right) \cdot \left(\cos a \cdot \cos b\right) - \left(\sin a \cdot \sin b\right) \cdot \left(\sin a \cdot \sin b\right)} \cdot \left(\cos a \cdot \cos b + \sin a \cdot \sin b\right)\right)double code(double r, double a, double b) {
return ((double) (r * ((double) (((double) sin(b)) / ((double) cos(((double) (a + b))))))));
}
double code(double r, double a, double b) {
return ((double) (r * ((double) (((double) (((double) sin(b)) / ((double) (((double) (((double) (((double) cos(a)) * ((double) cos(b)))) * ((double) (((double) cos(a)) * ((double) cos(b)))))) - ((double) (((double) (((double) sin(a)) * ((double) sin(b)))) * ((double) (((double) sin(a)) * ((double) sin(b)))))))))) * ((double) (((double) (((double) cos(a)) * ((double) cos(b)))) + ((double) (((double) sin(a)) * ((double) sin(b))))))))));
}



Bits error versus r



Bits error versus a



Bits error versus b
Results
Initial program 15.6
rmApplied cos-sum0.3
rmApplied flip--0.4
Applied associate-/r/0.4
Final simplification0.4
herbie shell --seed 2020126
(FPCore (r a b)
:name "r*sin(b)/cos(a+b), B"
:precision binary64
(* r (/ (sin b) (cos (+ a b)))))