/* example of usage of DiscreteMC.ys 

   random walk on 1,2,...,5; for i = 2,3,4, go to i-1,i,i+1 with
   probability p,1-2p,p; from 1 and 5 go to 2 and 4 */

Load("DiscreteMC.ys");

m:=ZeroMatrix(5,5);
ForEach(i,2 .. 4)
   [
   m[i][i-1]:=p;
   m[i][i+1]:=p;
   m[i][i]:=1-2*p;
   ];
m[1][2]:=1;
m[5][4]:=1;

pi:=solvedmc(m);

ForEach(i,1 .. 5) Echo(Simplify(pi[i]));