function [lambda1_koz, x1_koz, k] = hatvanyModszer(A, maxk, tol);
  %% Ellenõrzések:
  [nA, mA] = size(A);
  if nA != mA||!isnumeric(A)||!isnumeric(maxk)||!isnumeric(tol)||...
    maxk <1||tol < 0 
    error('Hibás input paraméterek!');
  endif
  
  %% Inicializálás:
  lambda1_koz = -inf;
  hiba = inf;
  x1_koz = rand(nA,1);
  x1_koz = x1_koz/norm(x1_koz);
  k = 0;
  
  while k <= maxk && hiba > tol
    k = k+1;
    x1_new = A*x1_koz;
    lambda1_new  = x1_new(1)/x1_koz(1);
    hiba = abs(lambda1_new - lambda1_koz);
    lambda1_koz = lambda1_new;
    x1_koz = x1_new/norm(x1_new);
  end
endfunction

