\frac{r \cdot \sin b}{\cos \left(a + b\right)}\frac{\sin b}{\cos a \cdot \cos b - \sqrt[3]{\left(\left(\sin a \cdot \sin b\right) \cdot \left(\sin a \cdot \sin b\right)\right) \cdot \left(\sin a \cdot \sin b\right)}} \cdot rdouble f(double r, double a, double b) {
double r1447660 = r;
double r1447661 = b;
double r1447662 = sin(r1447661);
double r1447663 = r1447660 * r1447662;
double r1447664 = a;
double r1447665 = r1447664 + r1447661;
double r1447666 = cos(r1447665);
double r1447667 = r1447663 / r1447666;
return r1447667;
}
double f(double r, double a, double b) {
double r1447668 = b;
double r1447669 = sin(r1447668);
double r1447670 = a;
double r1447671 = cos(r1447670);
double r1447672 = cos(r1447668);
double r1447673 = r1447671 * r1447672;
double r1447674 = sin(r1447670);
double r1447675 = r1447674 * r1447669;
double r1447676 = r1447675 * r1447675;
double r1447677 = r1447676 * r1447675;
double r1447678 = cbrt(r1447677);
double r1447679 = r1447673 - r1447678;
double r1447680 = r1447669 / r1447679;
double r1447681 = r;
double r1447682 = r1447680 * r1447681;
return r1447682;
}



Bits error versus r



Bits error versus a



Bits error versus b
Results
Initial program 15.0
rmApplied cos-sum0.3
rmApplied *-un-lft-identity0.3
Applied times-frac0.3
Simplified0.3
rmApplied add-cbrt-cube0.4
Applied add-cbrt-cube0.4
Applied cbrt-unprod0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2019120 +o rules:numerics
(FPCore (r a b)
:name "r*sin(b)/cos(a+b), A"
(/ (* r (sin b)) (cos (+ a b))))