%% Metoda Rayleghových podílů

A=[4 3 3 -2; -2 5 -1 -2; 4 -1 3 2; 0 -1 -3 4];
eig(A)
x=[1;1;1;1]; % pocatecni iterace
for k=1:10
    x1=A*x;
    [x x1]
    lambda=x'*x1/(x'*x)
    x=x1;
    pause(1)
end

%
%% Potlačení divergence

x=[1;1;1;1]; % pocatecni iterace
for k=1:10
    x=x/norm(x);
    x1=A*x;
    [x x1]
    lambda=x'*x1
    x=x1;
    pause(1)
end