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