\frac{2.30753 + x \cdot 0.27061000000000002}{1 + x \cdot \left(0.992290000000000005 + x \cdot 0.044810000000000003\right)} - x\frac{1}{\sqrt{1 + x \cdot \left(0.992290000000000005 + x \cdot 0.044810000000000003\right)}} \cdot \frac{\mathsf{fma}\left(0.27061000000000002, x, 2.30753\right)}{\sqrt{1 + x \cdot \left(0.992290000000000005 + x \cdot 0.044810000000000003\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 (((1.0 / sqrt((1.0 + (x * (0.99229 + (x * 0.04481)))))) * (fma(0.27061, x, 2.30753) / sqrt((1.0 + (x * (0.99229 + (x * 0.04481))))))) - x);
}



Bits error versus x
Results
Initial program 0.0
rmApplied add-sqr-sqrt0.0
Applied *-un-lft-identity0.0
Applied times-frac0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2020056 +o rules:numerics
(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))