fixed error in mu_n function call
trying p_ext as a variable in simple solver