2D nodal DPG method in Matlab