To run an example problem using the optimised code run hole_DG_solver.m This will solve a hole in an infinite plate, acting in plane stress. The optimised global stiffness matrix formulation occurs in Linear2D.m.