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