function [l,x]=met_tragova(A,epsilon)

l=sum(diag(A));
k=2;
A0=A;
while 1
   A=A^2;
   prev=l;
   l=sum(diag(A))^(1/k);
   if (abs(l-prev)<epsilon) 
       break;
   end
   k=2*k;
end
x=A*rand(max(size(A)),1);
x=x/norm(x);


% Odredjujemo znak sopstvene vrednosti tako da Ax=lx daje manju gresku
if (abs(A0*x-l*x)>abs(A0*x+l*x))
    l=-l;
end