\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 r24977 = r;
double r24978 = b;
double r24979 = sin(r24978);
double r24980 = r24977 * r24979;
double r24981 = a;
double r24982 = r24981 + r24978;
double r24983 = cos(r24982);
double r24984 = r24980 / r24983;
return r24984;
}
double f(double r, double a, double b) {
double r24985 = r;
double r24986 = b;
double r24987 = sin(r24986);
double r24988 = a;
double r24989 = cos(r24988);
double r24990 = cos(r24986);
double r24991 = r24989 * r24990;
double r24992 = 3.0;
double r24993 = pow(r24991, r24992);
double r24994 = sin(r24988);
double r24995 = r24994 * r24987;
double r24996 = pow(r24995, r24992);
double r24997 = r24993 - r24996;
double r24998 = r24987 / r24997;
double r24999 = r24985 * r24998;
double r25000 = r24991 * r24991;
double r25001 = r24995 * r24995;
double r25002 = r24991 * r24995;
double r25003 = r25001 + r25002;
double r25004 = r25000 + r25003;
double r25005 = r24999 * r25004;
return r25005;
}



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))))