pub fn psi_1piy_e(x: f64) -> (Value, Result)
This routine computes the real part of the digamma function on the line 1+i y, \Re[\psi(1 + i y)].
