adding new solvers

these solvers take the initial composition on entrance to calculate the
chemical potentials of electrons and neutrons, then the baryon number
density is the only variable
debugging 'phi', which is a recent addition