clear all a=1.5 A=[0 0 0 0 127 4 80; 0.6747*a 0.737 0 0 0 0 0; 0 0.0486 0.6610 0 0 0 0; 0 0 0.0147 0.6907 0 0 0; 0 0 0 0.0518 0 0 0; 0 0 0 0 0.8091 0 0; 0 0 0 0 0 0.8091 0.8089] [v lambda]=eig(A) %w=v(:,7)/sum(v(:,7))