
n=3;
A=hilb(n)
I=eye(3);
AI=kron(A,I)
IA=kron(I,A)

e=ones(n,1);
I=eye(n,n);
A=spdiags([e -2*e e],[-1 0 1],n,n);
A=full(A)
