powermethod.asv

来自「Numerical analysis （Gauss、Bisection、jaco」· ASV 代码 · 共 18 行

function s=powermethod(a,u,N)
if nargin<3
    return
end
i=0;
while i<N
    c=u(j);
    for j=1:length(u)-1
        if fabs(u(j))<fabs(u(j+1))
            c=u(j+1);
        end
    end
%    b=norm(a*u,inf);
    u=(a*u)/c;
    s=a*u;
    i=i+1;
end
return