\frac{r \cdot \sin b}{\cos \left(a + b\right)}r \cdot \frac{\sin b}{\cos a \cdot \cos b - \sqrt[3]{{\left(\sin b \cdot \sin a\right)}^{3}}}double f(double r, double a, double b) {
double r24663 = r;
double r24664 = b;
double r24665 = sin(r24664);
double r24666 = r24663 * r24665;
double r24667 = a;
double r24668 = r24667 + r24664;
double r24669 = cos(r24668);
double r24670 = r24666 / r24669;
return r24670;
}
double f(double r, double a, double b) {
double r24671 = r;
double r24672 = b;
double r24673 = sin(r24672);
double r24674 = a;
double r24675 = cos(r24674);
double r24676 = cos(r24672);
double r24677 = r24675 * r24676;
double r24678 = sin(r24674);
double r24679 = r24673 * r24678;
double r24680 = 3.0;
double r24681 = pow(r24679, r24680);
double r24682 = cbrt(r24681);
double r24683 = r24677 - r24682;
double r24684 = r24673 / r24683;
double r24685 = r24671 * r24684;
return r24685;
}



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 add-cbrt-cube0.4
Applied add-cbrt-cube0.4
Applied cbrt-unprod0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2019325 +o rules:numerics
(FPCore (r a b)
:name "r*sin(b)/cos(a+b), A"
:precision binary64
(/ (* r (sin b)) (cos (+ a b))))