Power Iterations

A=[ 3 1,

1 3]

[v,d]=eig(A);

x0=[0,

1];

xi=x0;

for i=1:100,

x1=xi;

i

xi=A*xi

lambda=norm(xi,inf)/norm(x1,inf)

$\%$ Normalize

xi=xi/norm(xi,2)

end;