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



Bits error versus r



Bits error versus a



Bits error versus b
Results
Initial program 15.3
rmApplied cos-sum0.3
Taylor expanded around inf 0.3
Simplified0.4
rmApplied clear-num0.4
Final simplification0.4
herbie shell --seed 2020182
(FPCore (r a b)
:name "rsin B"
:precision binary64
(* r (/ (sin b) (cos (+ a b)))))