scipy.special.zeta

scipy.special.zeta(x1, x2[, out])

y=zeta(x,q) returns the Riemann zeta function of two arguments: sum((k+q)**(-x),k=0..inf)