function [lambda V cnt]= potencia(A,x,tol,max1)
% Método de la potencia para la obtención de un valor
% y un vector propios de A
%
 if nargin<4 max1=100; end, if nargin<3 tol=eps; end
 if nargin<2 n=length(A); x=rand(n,1); end
 lambda=0;cnt=0;err=1;
 while cnt<=max1 && err>tol
    y=A*x;
    c1=norm(y,inf);
    y=y/c1;
    dc=abs(lambda-c1); dv=norm(x-y); err=max(dc,dv);
    x=y; lambda=c1;
    cnt=cnt+1;
 end
 V=x/norm(x);
end