trying p_ext constraint without a variable
added root find for mu_n from pressure_n
added print statement, analytical jacobian still needs pressure constraint
added pressure_n declaration