added solver without mu_n for the outer crust only
trying mu variable, but n_b constrained without equ