plot log(dirichlet eta(0.5+ s i) ) from 0 to 40