IPM.m

% 反幂法
% 注释的地方需要修改

function fanminfa
A = [3,1,2;1,3,4;2,4,6];           % 所求矩阵
[m,u]=powinv(A,1e-4,500);               % 第二个参数为精度，第三个为最大迭代次数
vpa(m,7)                           % 最大特征值有效数字位数
vpa(u,7)                           % 对应的特征向量有效数字位数

function [m,u]=powinv(A,ep,Nmax)
n=length(A);u=ones(n,1);k=0;m1=0;invA=inv(A);
while k<=Nmax
    v=invA*u;[vmax,i]=max(abs(v));
    m=v(i);u=v/m;
    if abs(m-m1)<ep
        break;
    end
    m1=m;k=k+1;
end
m=1/m;