n=3;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);
spy(A)
Q=full(A);
save('C:\STIMATH\Output\KronDelta.dat','Q','-ascii')
