%Octave:
%hold off
A=[.8 .2 .1;.1 .7 .3; .1 .1 .6]
axis([0 1 0 1]);axis('square')
plot(0,0);
%hold on
%title('Markov');xlabel('a');ylabel('b');grid
%nappula=1
%while nappula==1
 %  [a,b] = ginput
    [a,b]=[0.5,0.5]
    u=[a;b;1-a-b];
    x=u; k=0:7;
    while length(x)<=7
      u=A*u; x=[x u];
    endwhile
    % end
      plot(x(1,:),x(2,:),x(1,:),x(2,:),'o')
%endwhile
%end
hold off

      ////////////////

M=[0.95 0.03;0.05 0.97]
u0=[0.6;0.4]
u1=M*u0
u2=M*u1

