Updated 1st version working qualitatively for maxwell 2D with pressure as Lagrange multiplier.
Updated for newer version of matlab.