\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 ((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))))) * ((2.30753 + (x * 0.27061)) / sqrt(1.0 + (x * (0.99229 + (x * 0.04481)))))) - x;
}



Bits error versus x
Results
Initial program 0.0
rmApplied add-sqr-sqrt_binary64_7820.1
Applied *-un-lft-identity_binary64_7600.1
Applied times-frac_binary64_7660.1
Final simplification0.1
herbie shell --seed 2021093
(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))