clc;clear;
a=[0,0,0,1,0,1;1,0,0,0,0,0;0,1/2,0,0,0,0;0,1/2,1/3,0,0,0;0,0,1/3,0,0,0;0,0,1/3,0,1,0];
x=zeros(6,100);
x(:,1)=[1/6,1/6,1/6,1/6,1/6,1/6];
for i=1:99;
x(:,i+1)=a*x(:,i);
end
bar(x(:,100));
弄了一会，和老师做的答案貌似不一样，应该是因为模型的迭代矩阵有问题，看后面的马尔科夫链学习完以后有没有什么新的进展，这个暂且这样，大概迭代了30次左右达到稳定值，柱形图如下。
