x - \frac{2.30753 + x \cdot 0.27061000000000002}{1 + \left(0.992290000000000005 + x \cdot 0.044810000000000003\right) \cdot x}x - \frac{1}{\sqrt[3]{{\left(\frac{1 + \left(0.992290000000000005 + x \cdot 0.044810000000000003\right) \cdot x}{2.30753 + x \cdot 0.27061000000000002}\right)}^{3}}}double f(double x) {
double r83633 = x;
double r83634 = 2.30753;
double r83635 = 0.27061;
double r83636 = r83633 * r83635;
double r83637 = r83634 + r83636;
double r83638 = 1.0;
double r83639 = 0.99229;
double r83640 = 0.04481;
double r83641 = r83633 * r83640;
double r83642 = r83639 + r83641;
double r83643 = r83642 * r83633;
double r83644 = r83638 + r83643;
double r83645 = r83637 / r83644;
double r83646 = r83633 - r83645;
return r83646;
}
double f(double x) {
double r83647 = x;
double r83648 = 1.0;
double r83649 = 1.0;
double r83650 = 0.99229;
double r83651 = 0.04481;
double r83652 = r83647 * r83651;
double r83653 = r83650 + r83652;
double r83654 = r83653 * r83647;
double r83655 = r83649 + r83654;
double r83656 = 2.30753;
double r83657 = 0.27061;
double r83658 = r83647 * r83657;
double r83659 = r83656 + r83658;
double r83660 = r83655 / r83659;
double r83661 = 3.0;
double r83662 = pow(r83660, r83661);
double r83663 = cbrt(r83662);
double r83664 = r83648 / r83663;
double r83665 = r83647 - r83664;
return r83665;
}



Bits error versus x
Results
Initial program 0.0
rmApplied clear-num0.0
rmApplied add-cbrt-cube21.8
Applied add-cbrt-cube21.8
Applied cbrt-undiv21.8
Simplified0.0
Final simplification0.0
herbie shell --seed 2020060
(FPCore (x)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, D"
:precision binary64
(- x (/ (+ 2.30753 (* x 0.27061)) (+ 1 (* (+ 0.99229 (* x 0.04481)) x)))))