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]]______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-devel
