plot re (log(zeta(0.5+ d+4i)/zeta (0.5- d+4i))) -(-d log(4/2/pi) )from -1 to 1