function A = poisson(n)
I = speye(n,n);
 E = sparse(2:n,1:n-1,1,n,n);
 D = E+E'-2*I;
 A = -(kron(D,I)+kron(I,D));
