>> n=6;
>> A=zeros(2*n,2*n);
>> for ii=1:2:2*n for jj=2:2:2*n A(ii,jj)=ones(1); end, end
>> p=sum(A(1,:))
>> for jj=1:n x(1,jj)=p+jj-1; end
>> X=[x;ones(n-1,n)]
>> Y=cumsum(X);
>> A(2:2:end,1:2:end)=Y