plot log((zeta(-77.66- si) - Zeta(-77.66+si)))from -3 to 3