plot log(dirichlet eta (1 + i s)/zeta (1 + i s)) from 0 to 30