\frac{2.30753 + x \cdot 0.27061}{1 + x \cdot \left(0.99229 + x \cdot 0.04481\right)} - x\frac{1}{\sqrt{1 + x \cdot \left(0.99229 + x \cdot 0.04481\right)}} \cdot \frac{2.30753 + x \cdot 0.27061}{\sqrt{1 + x \cdot \left(0.99229 + x \cdot 0.04481\right)}} - x(FPCore (x) :precision binary64 (- (/ (+ 2.30753 (* x 0.27061)) (+ 1.0 (* x (+ 0.99229 (* x 0.04481))))) x))
(FPCore (x)
:precision binary64
(-
(*
(/ 1.0 (sqrt (+ 1.0 (* x (+ 0.99229 (* x 0.04481))))))
(/
(+ 2.30753 (* x 0.27061))
(sqrt (+ 1.0 (* x (+ 0.99229 (* x 0.04481)))))))
x))double code(double x) {
return ((double) ((((double) (2.30753 + ((double) (x * 0.27061)))) / ((double) (1.0 + ((double) (x * ((double) (0.99229 + ((double) (x * 0.04481))))))))) - x));
}
double code(double x) {
return ((double) (((double) ((1.0 / ((double) sqrt(((double) (1.0 + ((double) (x * ((double) (0.99229 + ((double) (x * 0.04481))))))))))) * (((double) (2.30753 + ((double) (x * 0.27061)))) / ((double) sqrt(((double) (1.0 + ((double) (x * ((double) (0.99229 + ((double) (x * 0.04481))))))))))))) - x));
}



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