pub fn zeta(x: f64) -> f64
This routine computes the Riemann zeta function \zeta(s) for arbitrary s, s \ne 1.