function ind=indize(i,j,nrow)
%function ind=indize(i,j,nrow)
ind=i'*ones(j/j(1))+ones(i/i(1))'*(j-1)*nrow;
ind=ind(:)';
