PSIFN

Algorithm 610: A portable Fortran subroutine for derivatives of the psi function.

This software is also peer reviewed by journal TOMS.