plot re(log((zeta(x+2 pi i)) / (Zeta(1-x+2 pi i)))) + i/1887.68 (2x-1)(1-(2x-1)^2) from 0 to 1