\frac{x - \sin x}{x - \tan x}\begin{array}{l}
\mathbf{if}\;x \le -0.02805121074086329088781432972155016614124:\\
\;\;\;\;\frac{x}{x - \tan x} - \frac{\sin x}{x - \tan x}\\
\mathbf{elif}\;x \le 1.586202501676462706114989487105049192905:\\
\;\;\;\;\frac{9}{40} \cdot {x}^{2} - \mathsf{fma}\left(\frac{27}{2800}, {x}^{4}, \frac{1}{2}\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\frac{x - \sin x}{x - \tan x}} \cdot \sqrt{\frac{x - \sin x}{x - \tan x}}\\
\end{array}double f(double x) {
double r29626 = x;
double r29627 = sin(r29626);
double r29628 = r29626 - r29627;
double r29629 = tan(r29626);
double r29630 = r29626 - r29629;
double r29631 = r29628 / r29630;
return r29631;
}
double f(double x) {
double r29632 = x;
double r29633 = -0.02805121074086329;
bool r29634 = r29632 <= r29633;
double r29635 = tan(r29632);
double r29636 = r29632 - r29635;
double r29637 = r29632 / r29636;
double r29638 = sin(r29632);
double r29639 = r29638 / r29636;
double r29640 = r29637 - r29639;
double r29641 = 1.5862025016764627;
bool r29642 = r29632 <= r29641;
double r29643 = 0.225;
double r29644 = 2.0;
double r29645 = pow(r29632, r29644);
double r29646 = r29643 * r29645;
double r29647 = 0.009642857142857142;
double r29648 = 4.0;
double r29649 = pow(r29632, r29648);
double r29650 = 0.5;
double r29651 = fma(r29647, r29649, r29650);
double r29652 = r29646 - r29651;
double r29653 = r29632 - r29638;
double r29654 = r29653 / r29636;
double r29655 = sqrt(r29654);
double r29656 = r29655 * r29655;
double r29657 = r29642 ? r29652 : r29656;
double r29658 = r29634 ? r29640 : r29657;
return r29658;
}



Bits error versus x
if x < -0.02805121074086329Initial program 0.0
rmApplied add-cube-cbrt1.4
Applied fma-neg1.4
rmApplied div-sub1.4
Simplified0.1
Simplified0.0
if -0.02805121074086329 < x < 1.5862025016764627Initial program 62.9
Taylor expanded around 0 0.1
Simplified0.1
if 1.5862025016764627 < x Initial program 0.0
rmApplied add-sqr-sqrt0.0
Final simplification0.1
herbie shell --seed 2019306 +o rules:numerics
(FPCore (x)
:name "sintan (problem 3.4.5)"
:precision binary64
(/ (- x (sin x)) (- x (tan x))))