plot |zeta (1+si)| from 0 to 30