%% Mocninná metoda

A=[1 2 3;4 5 6;7 8 9];
x=[1;1;1]; % poc. iterace
x1=A*x % prvni iterace, normujeme ji
M=max(abs(x1))
ii=find(M==abs(x1));
ii=ii(1); % index max. prvku
x=x1/M;
if x1(ii)<0
    x=-x;
end
x % prvni iterace

%
% 10 iterací

for k=1:10
    x1=A*x;
    M=max(abs(x1));
    ii=find(M==abs(x1));
    ii=ii(1); % index max. prvku
    x=x1/M;
    if x1(ii)<0
        x=-x;
    end
end
x
x1=A*x

%
% Porovnáme s výsledkem Matlabu

