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



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