function A=assempoi(n)
yksin=ones(n,1);
%yksinm1=ones(n-1,1);
DB=spdiags([-yksin  4*yksin -yksin],[-1 0 1], n,n);
%%%A=full(DB);
SDB=-speye(n,n);
multD=-SDB;
multS=spdiags([yksin  yksin],[-1 1], n,n);
A=kron(multD,DB)+kron(multS,SDB);
