\frac{x - \sin x}{x - \tan x}\begin{array}{l}
\mathbf{if}\;x \leq -1.5604285908945796:\\
\;\;\;\;\sqrt{\frac{x - \sin x}{x - \tan x}} \cdot \sqrt{\frac{x - \sin x}{x - \tan x}}\\
\mathbf{elif}\;x \leq 0.0363653530648353:\\
\;\;\;\;\left(x \cdot x\right) \cdot 0.225 - \left(0.5 + 0.009642857142857142 \cdot {x}^{4}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x - \tan x} - \frac{\sin x}{x - \tan x}\\
\end{array}(FPCore (x) :precision binary64 (/ (- x (sin x)) (- x (tan x))))
(FPCore (x)
:precision binary64
(if (<= x -1.5604285908945796)
(*
(sqrt (/ (- x (sin x)) (- x (tan x))))
(sqrt (/ (- x (sin x)) (- x (tan x)))))
(if (<= x 0.0363653530648353)
(- (* (* x x) 0.225) (+ 0.5 (* 0.009642857142857142 (pow x 4.0))))
(- (/ x (- x (tan x))) (/ (sin x) (- x (tan x)))))))double code(double x) {
return (x - sin(x)) / (x - tan(x));
}
double code(double x) {
double tmp;
if (x <= -1.5604285908945796) {
tmp = sqrt((x - sin(x)) / (x - tan(x))) * sqrt((x - sin(x)) / (x - tan(x)));
} else if (x <= 0.0363653530648353) {
tmp = ((x * x) * 0.225) - (0.5 + (0.009642857142857142 * pow(x, 4.0)));
} else {
tmp = (x / (x - tan(x))) - (sin(x) / (x - tan(x)));
}
return tmp;
}



Bits error versus x
Results
if x < -1.56042859089457964Initial program 0.0
rmApplied add-sqr-sqrt_binary640.0
if -1.56042859089457964 < x < 0.0363653530648352999Initial program 63.1
Taylor expanded around 0 0.1
Simplified0.1
if 0.0363653530648352999 < x Initial program 0.0
rmApplied div-sub_binary640.0
Final simplification0.1
herbie shell --seed 2020295
(FPCore (x)
:name "sintan (problem 3.4.5)"
:precision binary64
(/ (- x (sin x)) (- x (tan x))))