solve mu = n * exp(theta)/(exp(theta)+1) for theta