n=4;b=1:n*n;
B=reshape(b,m,n)
for i=1:2:n, B(i,i)=-1;end
B
BB=reshape(b,m,n);
I=1:2:n;
K=sub2ind(size(BB),I,I)
BB(K)=-1
