plot re (log((zeta(-15.5+ s i)) / (Zeta(1+15.5+ s i)))) from -9 to 9