A=[0.4,0.4;0.2,0.1] B=[0.1;0.2] U=[0.5;0.3] for n=1:5, U=A*U+B, end