number of particles as default input, chemical potential can be given as initial guess
bugs fixed