n=6;
for i=1:n
    for j=1:n
        A(i,j)=1/(i+j-1);
    end
end

x=ones(n,1);

power_method(A,x,10^(-4),100)