%File: p1.m created Feb. 4, 2000 by Fayez Gebali
%The program finds the steady state
%eigenvector for a 2X2 state transition matrix
%
P = [0.5 0.75
0.5 0.25]
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 determine A_i
A(:,:,i) = X*sel*X_inv;
end
A