\sqrt[3]{x + 1} - \sqrt[3]{x}\frac{1}{\mathsf{fma}\left(\sqrt[3]{x + 1}, \sqrt[3]{x + 1}, \sqrt[3]{\sqrt[3]{x} \cdot \sqrt[3]{x}} \cdot \left(\sqrt[3]{\sqrt[3]{x}} \cdot \left(\sqrt[3]{x + 1} + \sqrt[3]{x}\right)\right)\right)}double f(double x) {
double r80689 = x;
double r80690 = 1.0;
double r80691 = r80689 + r80690;
double r80692 = cbrt(r80691);
double r80693 = cbrt(r80689);
double r80694 = r80692 - r80693;
return r80694;
}
double f(double x) {
double r80695 = 1.0;
double r80696 = x;
double r80697 = r80696 + r80695;
double r80698 = cbrt(r80697);
double r80699 = cbrt(r80696);
double r80700 = r80699 * r80699;
double r80701 = cbrt(r80700);
double r80702 = cbrt(r80699);
double r80703 = r80698 + r80699;
double r80704 = r80702 * r80703;
double r80705 = r80701 * r80704;
double r80706 = fma(r80698, r80698, r80705);
double r80707 = r80695 / r80706;
return r80707;
}



Bits error versus x
Initial program 29.9
rmApplied flip3--29.9
Simplified29.3
Simplified29.3
Taylor expanded around 0 0.5
rmApplied add-cube-cbrt0.5
Applied cbrt-prod0.6
Applied associate-*l*0.6
Final simplification0.6
herbie shell --seed 2020047 +o rules:numerics
(FPCore (x)
:name "2cbrt (problem 3.3.4)"
:precision binary64
(- (cbrt (+ x 1)) (cbrt x)))