scipy.special.zetac

scipy.special.zetac(x[, out])

y=zetac(x) returns 1.0 - the Riemann zeta function: sum(k**(-x), k=2..inf)