gam1d function

Compute 1/Gamma(x+1) - 1 Accurately