%The program finds the steady state

%eigenvector for a 3X3 state transition matrix

%

P = (1/3)*[1 1 1

1 1 1

1 1 1]

m = size(P,1);

[X,D] = eig(P)

X_inv = inv(X)
 
 

for i=1:m,
 
 

sel = zeros(m);
 
 

sel(i,i) = 1; % sel is a matrix to deterine A_i
 
 

A(:,:,i) = X*sel*X_inv;
 
 

end

A