pub fn psi_e(x: f64) -> (Value, Result)
This routine computes the digamma function \psi(x) for general x, x \ne 0.