solver with better initial guess

exponential goes to infinity as solver tries larger values for chemical
potentials
added table of Andrew's equilibrium distribution
added residual tolerance subroutine