function [lambda,x,k]=Pow_Meth(A,x0,eps,N)
%幂法求成对比较阵A的主特征值lambda及其特征向量x
%x0为迭代初始向量,eps为预先给定的精度,N为最大迭代次数
%k为迭代次数
n=length(x0);
x0=x0/sum(x0); %初始向量x0归一化
for k=1:N
x1=A*x0;
x=x1/sum(x1);%归一化
err=max(abs(x-x0));
if err<=eps
break;
end
x0=x;
end
lambda=sum(x1./x0)/n;