Лаба 3 (Цепи Маркова и СМО) / st

function z = st(N,P_cum)
z=zeros(N,1);
z(1,1)=1;
r=rand(N,1);
n=2;
k=1;
while n~=N+1
if (k==1)&&(r(n-1,1)>=P_cum(1,1))
k=2;
z(n,1)=k;
n=n+1;
elseif (k==1)&&(r(n-1,1)<P_cum(1,1))
k=1;
z(n,1)=k;
n=n+1;
end
if (k==2)&&(r(n-1,1)>=P_cum(2,2))
k=3;
z(n,1)=k;
n=n+1;
elseif (k==2)&&(r(n-1,1)<=P_cum(2,2))
k=2;
z(n,1)=k;
n=n+1;
end
if (k==3)&&(r(n-1,1)>=P_cum(3,2))
k=3;
z(n,1)=k;
n=n+1;
elseif (k==3)&&(r(n-1,1)<P_cum(3,2))
k=2;
z(n,1)=k;
n=n+1;
end
if (k==4)&&(r(n-1,1)>=P_cum(4,3))
k=4;
z(n,1)=k;
n=n+1;
elseif (k==4)&&(r(n-1,1)<P_cum(4,3))
k=3;
z(n,1)=k;
n=n+1;
end
if (n==N)
z(n,1)=z(n-1,1);
break
end
end
x =1:N;
stairs(x,z);
end