function L2 = nma_lap2d(nx,ny) %helper function to make sparse 2D matrix Lx=nma_lap1d(nx); Ly=nma_lap1d(ny); Ix=speye(nx); Iy=speye(ny); L2=kron(Iy,Lx)+kron(Ly,Ix); end