%  chain.m simulates a Markov chain on {1, 2, ..., n} given an initial %  distribution and a transition matrix. % Program to simulate a Markov chain% You might want to set up another .m file that defines mu and P% The program assumes that the states are labeled 1, 2, ...% Below is a sample, which you can use by uncommenting these lines% mu=[1 0 0];     % initial distribution% P=[[.6 .3 .1]; [.3 .3 .4]; [.4 .1 .5]];  % transition matrixn=80;           % number of time steps to takex=zeros(1,n+1); % clear out any old valuest=0:n;          % time indicesx(1)=rando(mu); % generate first x value (time 0, not time 1)for i=1:n,  x(i+1) = rando(P(x(i),:));endplot(t, x, '*');axis([0 n 0 (length(mu)+1)]);