clear all
d=0.2
A=[1.5*(1-d) 0.5*d 0.1*d; 0.75*d 1-d 0.1*d; 0.75*d 0.5*d 0.2*(1-d)]
[v lambda]=eig(A)

w=v(:,1)/sum(v(:,1))