scipy.special.zeta

scipy.special.zeta(x, q) returns the Riemann zeta function of two arguments: sum((k+q)**(-x), k=0..inf) = <ufunc 'zeta'>