\sqrt[3]{x + 1} - \sqrt[3]{x}\begin{array}{l}
\mathbf{if}\;x \leq -58140.775394494696 \lor \neg \left(x \leq 75523.64683228085\right):\\
\;\;\;\;\frac{\sqrt[3]{x}}{x} \cdot \left(0.3333333333333333 + \frac{-0.1111111111111111}{x}\right) + \left(\sqrt[3]{x} - \sqrt[3]{x}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt[3]{1 + {x}^{3}}}{\sqrt[3]{x \cdot x + \left(1 - x\right)}} - \sqrt[3]{x}\\
\end{array}(FPCore (x) :precision binary64 (- (cbrt (+ x 1.0)) (cbrt x)))
(FPCore (x)
:precision binary64
(if (or (<= x -58140.775394494696) (not (<= x 75523.64683228085)))
(+
(* (/ (cbrt x) x) (+ 0.3333333333333333 (/ -0.1111111111111111 x)))
(- (cbrt x) (cbrt x)))
(- (/ (cbrt (+ 1.0 (pow x 3.0))) (cbrt (+ (* x x) (- 1.0 x)))) (cbrt x))))double code(double x) {
return cbrt(x + 1.0) - cbrt(x);
}
double code(double x) {
double tmp;
if ((x <= -58140.775394494696) || !(x <= 75523.64683228085)) {
tmp = ((cbrt(x) / x) * (0.3333333333333333 + (-0.1111111111111111 / x))) + (cbrt(x) - cbrt(x));
} else {
tmp = (cbrt(1.0 + pow(x, 3.0)) / cbrt((x * x) + (1.0 - x))) - cbrt(x);
}
return tmp;
}



Bits error versus x
Results
if x < -58140.7753944946962 or 75523.646832280851 < x Initial program 60.4
Taylor expanded around -inf 64.0
Simplified0.7
if -58140.7753944946962 < x < 75523.646832280851Initial program 0.2
rmApplied flip3-+_binary64_4170.2
Applied cbrt-div_binary64_4430.2
Simplified0.2
Simplified0.2
Final simplification0.4
herbie shell --seed 2020299
(FPCore (x)
:name "2cbrt (problem 3.3.4)"
:precision binary64
(- (cbrt (+ x 1.0)) (cbrt x)))