% Generovani realizaci Markovskeho retezce 
% vektor pocatecnich pravdepodobnosti p(0)=[1/2 1/3 1/6]
% matice prechodu: P=[1/3 2/3 0;1/2 1/4 1/4;0 3/4 1/4]
% generovani n realizaci nahodne veliciny s rozlozenim Rs(0,1):
n=8;
%u=unifrnd(0,1,8,1);
u=[0.5466 0.8948 0.1474 0.366 0.7574 0.2307 0.8395 0.5207]
% pocatecni krok:
if u(1)<1/2 h0=0;end
if u(1)>=1/2 & u(1)<=5/6 h0=1;end
if u(1)>=5/6 h0=2;end
% první krok
if h0==0 & u(2)<1/3 h1=0;end
if h0==0 & u(2)>=1/3 h1=1;end
if h0==1 & u(2)<1/2 h1=0;end
if h0==1 & u(2)>=1/2 & u(2)<3/4 h1=1;end
if h0==1 & u(2)>=3/4 h1=2;end
if h0==2 & u(2)<3/4 h1=1;end
if h0==2 & u(2)>=3/4 h1=2;end
% druhy krok
if h1==0 & u(3)<1/3 h2=0;end
if h1==0 & u(3)>=1/3 h2=1;end
if h1==1 & u(3)<1/2 h2=0;end
if h1==1 & u(3)>=1/2 & u(3)<3/4 h2=1;end
if h1==1 & u(3)>=3/4 h2=2;end
if h1==2 & u(3)<3/4 h2=1;end
if h1==2 & u(3)>=3/4 h2=2;end
% treti krok
if h2==0 & u(4)<1/3 h3=0;end
if h2==0 & u(4)>=1/3 h3=1;end
if h2==1 & u(4)<1/2 h3=0;end
if h2==1 & u(4)>=1/2 & u(4)<3/4 h3=1;end
if h2==1 & u(4)>=3/4 h3=2;end
if h2==2 & u(4)<3/4 h3=1;end
if h2==2 & u(4)>=3/4 h3=2;end
% ctvrty krok
if h3==0 & u(5)<1/3 h4=0;end
if h3==0 & u(5)>=1/3 h4=1;end
if h3==1 & u(5)<1/2 h4=0;end
if h3==1 & u(5)>=1/2 & u(5)<3/4 h4=1;end
if h3==1 & u(5)>=3/4 h4=2;end
if h3==2 & u(5)<3/4 h4=1;end
if h3==2 & u(5)>=3/4 h4=2;end
% paty krok
if h4==0 & u(6)<1/3 h5=0;end
if h4==0 & u(6)>=1/3 h5=1;end
if h4==1 & u(6)<1/2 h5=0;end
if h4==1 & u(6)>=1/2 & u(6)<3/4 h5=1;end
if h4==1 & u(6)>=3/4 h5=2;end
if h4==2 & u(6)<3/4 h5=1;end
if h4==2 & u(6)>=3/4 h5=2;end
% sesty krok
if h5==0 & u(7)<1/3 h6=0;end
if h5==0 & u(7)>=1/3 h6=1;end
if h5==1 & u(7)<1/2 h6=0;end
if h5==1 & u(7)>=1/2 & u(7)<3/4 h6=1;end
if h5==1 & u(7)>=3/4 h6=2;end
if h5==2 & u(7)<3/4 h6=1;end
if h5==2 & u(7)>=3/4 h6=2;end
% sedmy krok
if h6==0 & u(8)<1/3 h7=0;end
if h6==0 & u(8)>=1/3 h7=1;end
if h6==1 & u(8)<1/2 h7=0;end
if h6==1 & u(8)>=1/2 & u(8)<3/4 h7=1;end
if h6==1 & u(8)>=3/4 h7=2;end
if h6==2 & u(8)<3/4 h7=1;end
if h6==2 & u(8)>=3/4 h7=2;end
realizace=[h0 h1 h2 h3 h4 h5 h6 h7]
plot(realizace,'o')
axis([-1 9 -0.2 2.2])
