Jonathan Leto writes: > Why is the error so large for gamma(x>18) ? It's computed with an approximation so the relative error is the appropriate measure. Note that 1/Gamma(19) is the first value which is less than DBL_EPSILON. > Would it be possible to just return factorial(x-1) if gamma is > given an integer argument? Sounds reasonable to me.