/*幂法程序，函数名：pow.m*/functionm,u,index=pow(A,ep,N)%A 为矩阵； ep 为精度要求；N 为最大迭代次数；m 为绝对值最大的特征值；u 为对应最大特征值的特征向量。N=100;ep=1e-6;n=length(A);u=ones(n,1);index=0;k=0;m1=0;while k=Nv=A*u;vmax,i=max(abs(v);m=v(i);u=v/m;if abs(m-m1)epindex=1;break;end m1=m;k=k+1;end输入 A=2 -1 0;-1 2 -1;0 -1 2;m,u,index=pow(A,1e-6)
