\frac{r \cdot \sin b}{\cos \left(a + b\right)}\left(r \cdot \frac{\sin b}{{\left(\cos a \cdot \cos b\right)}^{3} - {\left(\sin a \cdot \sin b\right)}^{3}}\right) \cdot \left(\left(\cos a \cdot \cos b\right) \cdot \left(\cos a \cdot \cos b\right) + \left(\left(\sin a \cdot \sin b\right) \cdot \left(\sin a \cdot \sin b\right) + \left(\cos a \cdot \cos b\right) \cdot \left(\sin a \cdot \sin b\right)\right)\right)double f(double r, double a, double b) {
double r24976 = r;
double r24977 = b;
double r24978 = sin(r24977);
double r24979 = r24976 * r24978;
double r24980 = a;
double r24981 = r24980 + r24977;
double r24982 = cos(r24981);
double r24983 = r24979 / r24982;
return r24983;
}
double f(double r, double a, double b) {
double r24984 = r;
double r24985 = b;
double r24986 = sin(r24985);
double r24987 = a;
double r24988 = cos(r24987);
double r24989 = cos(r24985);
double r24990 = r24988 * r24989;
double r24991 = 3.0;
double r24992 = pow(r24990, r24991);
double r24993 = sin(r24987);
double r24994 = r24993 * r24986;
double r24995 = pow(r24994, r24991);
double r24996 = r24992 - r24995;
double r24997 = r24986 / r24996;
double r24998 = r24984 * r24997;
double r24999 = r24990 * r24990;
double r25000 = r24994 * r24994;
double r25001 = r24990 * r24994;
double r25002 = r25000 + r25001;
double r25003 = r24999 + r25002;
double r25004 = r24998 * r25003;
return r25004;
}



Bits error versus r



Bits error versus a



Bits error versus b
Results
Initial program 14.7
rmApplied cos-sum0.3
rmApplied *-un-lft-identity0.3
Applied times-frac0.3
Simplified0.3
rmApplied flip3--0.4
Applied associate-/r/0.4
Applied associate-*r*0.4
Final simplification0.4
herbie shell --seed 2019323
(FPCore (r a b)
:name "r*sin(b)/cos(a+b), A"
:precision binary64
(/ (* r (sin b)) (cos (+ a b))))