pub fn psi_n_e(n: i32, x: f64) -> (Value, Result)
This routine computes the polygamma function \psi^{(n)}(x) for n >= 0, x > 0.