\frac{2.30753 + x \cdot 0.27061000000000002}{1 + x \cdot \left(0.992290000000000005 + x \cdot 0.044810000000000003\right)} - x\frac{2.30753 + x \cdot 0.27061000000000002}{1 + \log \left(e^{x \cdot \left(0.992290000000000005 + x \cdot 0.044810000000000003\right)}\right)} - xdouble code(double x) {
return (((2.30753 + (x * 0.27061)) / (1.0 + (x * (0.99229 + (x * 0.04481))))) - x);
}
double code(double x) {
return (((2.30753 + (x * 0.27061)) / (1.0 + log(exp((x * (0.99229 + (x * 0.04481))))))) - x);
}



Bits error versus x
Results
Initial program 0.0
rmApplied add-log-exp0.2
Final simplification0.2
herbie shell --seed 2020058
(FPCore (x)
:name "Numeric.SpecFunctions:invIncompleteGamma from math-functions-0.1.5.2, C"
:precision binary64
(- (/ (+ 2.30753 (* x 0.27061)) (+ 1 (* x (+ 0.99229 (* x 0.04481))))) x))