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

xv(:,1)=x/norm(x,2);
k=1;
flagct=0;
while flagct==0,
      y=A*xv(:,k);
      rv(k)=y'*xv(:,k);
      xv(:,k+1)=y/norm(y,2);
      xn(k)=norm(xv(:,k+1)-xv(:,k),2);
      if xn(k)>tol,
          k=k+1;
          if k>kmax,
              flagct=1;
              flag=1;
              disp('Symmetric 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('Symmetric Power Method converges')
      end
end
xv
rv