plot abs(zeta(1/2+is)) + arg(zeta(1/2+is))i from 0 to 30