r \cdot \frac{\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 r24725 = r;
double r24726 = b;
double r24727 = sin(r24726);
double r24728 = a;
double r24729 = r24728 + r24726;
double r24730 = cos(r24729);
double r24731 = r24727 / r24730;
double r24732 = r24725 * r24731;
return r24732;
}
double f(double r, double a, double b) {
double r24733 = r;
double r24734 = b;
double r24735 = sin(r24734);
double r24736 = a;
double r24737 = cos(r24736);
double r24738 = cos(r24734);
double r24739 = r24737 * r24738;
double r24740 = 3.0;
double r24741 = pow(r24739, r24740);
double r24742 = sin(r24736);
double r24743 = r24742 * r24735;
double r24744 = pow(r24743, r24740);
double r24745 = r24741 - r24744;
double r24746 = r24735 / r24745;
double r24747 = r24733 * r24746;
double r24748 = r24739 * r24739;
double r24749 = r24743 * r24743;
double r24750 = r24739 * r24743;
double r24751 = r24749 + r24750;
double r24752 = r24748 + r24751;
double r24753 = r24747 * r24752;
return r24753;
}



Bits error versus r



Bits error versus a



Bits error versus b
Results
Initial program 14.7
rmApplied cos-sum0.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), B"
:precision binary64
(* r (/ (sin b) (cos (+ a b)))))