Matlab实现乘幂法

function [u,r,res]=myPM(A,u)
res=1;
while (res>1e-6)
    y=A*u;  
    [a,b]=max(abs(y));
    r=y(b,1);
    u=y/a;
    res=norm(A*u-r*u);
end