Added plotting routines to show boundary condition etc
Added mu on two BCs plot