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) cos(a)) / ((double) tan(b)))) - ((double) sin(a))))));
}



Bits error versus r



Bits error versus a



Bits error versus b
Results
Initial program 15.5
rmApplied cos-sum0.3
rmApplied associate-*r/0.4
Simplified0.4
Taylor expanded around inf 0.4
Simplified0.4
rmApplied quot-tan0.4
Final simplification0.4
herbie shell --seed 2020190
(FPCore (r a b)
:name "rsin B"
:precision binary64
(* r (/ (sin b) (cos (+ a b)))))