\frac{x - \sin x}{x - \tan x}\begin{array}{l}
\mathbf{if}\;x \leq -0.029873394671195968 \lor \neg \left(x \leq 0.02727798752832581\right):\\
\;\;\;\;\frac{\frac{x}{x - \tan x} \cdot \frac{x}{x - \tan x} - \frac{{\left(\sin x\right)}^{2}}{\left(x - \tan x\right) \cdot \left(x - \tan x\right)}}{\frac{x}{x - \tan x} + \frac{\sin x}{x - \tan x}}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot 0.225\right) + \left(-0.5 + {x}^{4} \cdot -0.009642857142857142\right)\\
\end{array}(FPCore (x) :precision binary64 (/ (- x (sin x)) (- x (tan x))))
(FPCore (x)
:precision binary64
(if (or (<= x -0.029873394671195968) (not (<= x 0.02727798752832581)))
(/
(-
(* (/ x (- x (tan x))) (/ x (- x (tan x))))
(/ (pow (sin x) 2.0) (* (- x (tan x)) (- x (tan x)))))
(+ (/ x (- x (tan x))) (/ (sin x) (- x (tan x)))))
(+ (* x (* x 0.225)) (+ -0.5 (* (pow x 4.0) -0.009642857142857142)))))double code(double x) {
return (((double) (x - ((double) sin(x)))) / ((double) (x - ((double) tan(x)))));
}
double code(double x) {
double tmp;
if (((x <= -0.029873394671195968) || !(x <= 0.02727798752832581))) {
tmp = (((double) (((double) ((x / ((double) (x - ((double) tan(x))))) * (x / ((double) (x - ((double) tan(x))))))) - (((double) pow(((double) sin(x)), 2.0)) / ((double) (((double) (x - ((double) tan(x)))) * ((double) (x - ((double) tan(x))))))))) / ((double) ((x / ((double) (x - ((double) tan(x))))) + (((double) sin(x)) / ((double) (x - ((double) tan(x))))))));
} else {
tmp = ((double) (((double) (x * ((double) (x * 0.225)))) + ((double) (-0.5 + ((double) (((double) pow(x, 4.0)) * -0.009642857142857142))))));
}
return tmp;
}



Bits error versus x
Results
if x < -0.0298733946711959679 or 0.02727798752832581 < x Initial program Error: 0.0 bits
rmApplied div-subError: 0.0 bits
rmApplied flip--Error: 0.1 bits
SimplifiedError: 0.1 bits
if -0.0298733946711959679 < x < 0.02727798752832581Initial program Error: 63.2 bits
Taylor expanded around 0 Error: 0.0 bits
SimplifiedError: 0.0 bits
Final simplificationError: 0.0 bits
herbie shell --seed 2020204
(FPCore (x)
:name "sintan (problem 3.4.5)"
:precision binary64
(/ (- x (sin x)) (- x (tan x))))