function [xsol,rv,flag,k]=PowerMethod(A,x,n,tol,kmax)

xv(:,1)=x/norm(x,inf);
k=1;
flagct=0;
while flagct==0,
      y=A*xv(:,k);
      absy=abs(y);
      [r,I]=max(absy);
      rv(k)=y(I);
      xv(:,k+1)=y/y(I);
      xn(k)=norm(xv(:,k+1)-xv(:,k),inf);
      if xn(k)>tol,
          k=k+1;
          if k>kmax,
              flagct=1;
              flag=1;
              disp('Power Method does not converge in kmax iterations')
          end
      else
          xsol=xv(:,k+1);
          lambda=rv(k);
          flagct=1;
          flag=0;
          k=k-1;
          disp('Power Method converges')
      end
end