n=10; R = zeros(3*n,3*n); R_temp = [1 2 3;4 5 6;7 8 9]; for i=1:3:3*n-2; R(i:i+2,i:i+2) = R_temp; end