for i:10
     R1 = [1,2;3,4]; Matrix is always 2x2 but different values each iteration
     cov = blkdiag(R1);
end

[ R1,0,0,0...,
   0,R1,0,0...]