pub fn psi_int(n: i32) -> f64
This routine computes the digamma function \psi(n) for positive integer n. The digamma function is also called the Psi function.