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