%  set up a n by n laplace matrix T
%
T=zeros(n);
for i=1:n-1
  T(i,i+1)=1;
  T(i+1,i)=1;
  T(i,i)=-2;
end
T(n,n)=-2;