double code(double x) {
return ((double) (x - ((double) (((double) (2.30753 + ((double) (x * 0.27061)))) / ((double) (1.0 + ((double) (((double) (0.99229 + ((double) (x * 0.04481)))) * x))))))));
}
double code(double x) {
return ((double) (x - ((double) (((double) (1.0 / ((double) sqrt(((double) (1.0 + ((double) (x * ((double) (0.99229 + ((double) (x * 0.04481)))))))))))) * ((double) (((double) (2.30753 + ((double) (x * 0.27061)))) / ((double) sqrt(((double) (1.0 + ((double) (x * ((double) (0.99229 + ((double) (x * 0.04481))))))))))))))));
}



Bits error versus x
Results
Initial program 0.0
rmApplied add-sqr-sqrt0.1
Applied *-un-lft-identity0.1
Applied times-frac0.1
Simplified0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2020180
(FPCore (x)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, D"
:precision binary64
(- x (/ (+ 2.30753 (* x 0.27061)) (+ 1.0 (* (+ 0.99229 (* x 0.04481)) x)))))