function L=nma_lap1d(n)
%helper function to make sparse matrix

e=ones(n,1);
L=spdiags([e -3*e e],[-1 0 1],n,n);
end