function matrix=remat2(array,ilen,jlen);
matrix=( ones(ilen,1)*(0:jlen-1)*ilen+(1:ilen)'*ones(1,jlen) );
matrix(matrix)=array(:);
