function inversepowermethod(alpha,A,x,numtimes)

n= length(x);
for (r=1:numtimes)
  reduced = rref( [(A-alpha*eye(n,n)) x]);
  y = reduced(:,n+1);

  maxnr = 1;
  for (k=2:n)  
    if ( abs(y(k)) > abs(y(maxnr)) )
      maxnr = k;
    end
  end
  mu = y(maxnr);
  nu = alpha + (1/mu)
  x = (1/mu)*y
end