plot re(log((zeta(x+2.002117 pi i)) / (Zeta(1-x+2.002117 pi i)))) - i/236.31 (x-0.5)^3 from 0 to 1