PSIFN: a Fortran subprogram for calculating derivatives of the psi function. A double precision version is included. (See D.E. Amos, ACM TOMS 9 (1983) pp. 494-502.).