contourplot[arg(zeta(x+y*I)^2)*8/Pi,{x=-5.0,5.0},{y=-22.0,22.0}]