B = rand(4,4); A = B*B';
[V,D]  = eig(A); eigs = diag(D)

eigs =
    0.0041
    0.1935
    0.3852
    3.7365

VEKTORITERATION: der groesste Eigenwert
setzt sich durch: 

x = rand(4,1); x = A*x; x = x/norm(x); x.'
ans =
    0.6357    0.6206    0.3791    0.2590
>> x = A*x; x = x/norm(x); x.'
ans =
    0.6616    0.5906    0.3659    0.2820
>> x = A*x; x = x/norm(x); x.'
ans =
    0.6643    0.5873    0.3649    0.2841
>> x = A*x; x = x/norm(x); x.'
ans =
    0.6646    0.5869    0.3648    0.2843
>> x = A*x; x = x/norm(x); x.'
ans =
    0.6646    0.5869    0.3648    0.2843

ZUM VERGLEICH, der zum groessten
Eigenwert gehoerige Eigenvektor (normiert) ist,
das ist die vierte Spalte von V, als Zeilenvektor
dargestellt ergibt das:  

>> V(:,4)'
ans =
    0.6646    0.5869    0.3648    0.2843