Fixed some numerical difficulties in the chemical potential
Just trying to speed up the code a little
Documented the entropy minimiser function