plot (log((zeta(-3+ s i)) / (Zeta(1+3+ s i)))) from -30 to 30