res[(2pi)^(1-s)*Gamma(s-1)*(zeta(s-1, 1/2 - i/2) - zeta(s-1, 1 - i/2))*zeta(s), {s,0}]