plot log((zeta(0.5 - si)Zeta(0.5 +si)))from -23 to 23