P=[0,1/3,1/3,1/3;0,0,0,1;0.5,0,0,0.5;0,0,1,0]
P^1000
Q=[0,1/3,1/3;0,0,0;0.5,0,0]
(eye(3)-Q)^(-1)
R=[0,1/3,1/3,1/3;0,0,0,1;0.5,0,0,0.5;1,0,0,0]
[v1,v2]=eig(R')
v1(:,1)
((v1(:,1)')/sum(v1(:,1))) /0.4
1.2*(((v1(:,1)')/sum(v1(:,1))) /0.4)

S=[0,0.5;1,0]
((eye(2)-S)^(-1)) * [0.5,0;0,0]

T=[1,0;0,1]
T^10
T^11

