x = exp(2565855.5315) ; psi(x) = Chebyshev function https://en.wikipedia.org/wiki/Chebyshev_function |psi(x)-x| ~= 1.14 sqrt(x) . Using: s1000 = (U)->-sum(X=1,1000,exp((I*w[X])*U)/(1/2+I*w[X])+exp((-I*w[X])*U)/(1/2-I*w[X])) w is a vector of the imaginary parts of the non-trivial zeta zeros with Im(rho)>0, i.e. rho_n = 1/2 + i*w[n] , n = 1… 100,000. ? rec = 0.1 ; for(X=0,0, z=2565855.5315 ; aa=real(s1000(z));… Continue reading Large values of psi(x), Chebychev function