More generally, the accuracy and working range of psigamma(x,deriv) can be improved by having it return the leading term of the asymptotic expansion (-1)^(deriv-1)*factorial(deriv-1)/x^deriv whenever deriv>=1 and x>=1e15 [[alternative HTML version deleted]]