I am trying to make the following code work.
t = var('t')
f = zeta(1/2+i*t).abs()
ff = fast_callable(f, vars=[t], domain=CDF)
print find_root(ff, 0, 40)
There are actually 6 roots between 0 and 40. But find_root could not find any of them. Is there any walkaround?