trying solver with pressure condition

n_b as variable, mu_e and mu_n determined
converging solver, solutions wrong when P<10^-4