function [eigValue,eigVector] = EigPowerMethod (A, epsilon)
[m,n]=size(A);
if (m != n)
  disp("Matrix Should Be Square");
else
  u = [];
  t = 0;
  r = ones(m);
  u = [u transpose(r(1,:))];
  i = 2;
  d = 1;
  while(d > epsilon)
    newU = A*u(:,i-1);
    t = newU(m);
    newU = newU/t;
    u =[u newU];
    d = norm(u(:,i)-u(:,i-1));
    i++;
  end
  fprintf("repeats loop count : %d \r\n",i)
  eigValue=t;
  eigVector=u(:,i-1);
endif
endfunction