function [y,nu,iter]=hatvanymodszer(A,kmax,toll);
% A program az A mátrix abszolút értékben domináns sajátértékét
% és a hozzá tartozó sajátvektort adja meg. kmax a maximális 
% iterációszám, toll pedig az a toleranciaszint, aminél, ha kevesebbel tér
% el két egymás utáni sajátértékközelítés, akkor leáll a program futása.
% Eredményképpen a sajátvektorközelítést (y), a sajátértékközelítést (nu) és az
% iterációszámot (iter) kapjuk meg.

n=max(size(A));   
   y=(ones(n,1));
   y=y/norm(y);
   nu=y'*A*y;
   err=1;
   iter=0;
while err>toll & iter<kmax
   iter=iter+1;
   y=A*y;
   y=y/norm(y);
   nuold=nu;
   nu=y'*A*y;
   err=abs(nu-nuold);
end;