To find the zeros of the Reimann \zeta with real part 1/2, and imaginary part less than 1000 in norm, the following is what I’m using at the moment,

using SpecialFunctions
using Roots
g(x) = abs(zeta(.5 + x*im))
Z =[]
for i in 1:5:1000
append!(Z,find_zeros(g,i,i+5))
end
E=[-reverse(Z)...,Z...] # by symmetric

I’m breaking the line [0,1000] into segments since find_zeros(g,0,1000) will miss many of the zeros.
But this looks a bit awkward, I wonder what will be a better way to find the zeros.