clc;
clear all;
a=[5 4 1;
4 3 2;
1 2 5];
A=a;

x0=[1;1;1];

nit=10;

eigval=0;
siz=size(A);
siz=siz(1);
for m=1:siz
x = x0;t=x0;
for n = 1:nit
xnew = A*x;
lambda = norm(xnew,inf);
 if m==3
     lambda=-lambda;
end
tnew=transpose(A)*t;
x1=x;
x=xnew/lambda;
x2=x;
t1=t;
t=tnew/lambda;
t2=t;
end
x = x/(norm(x));
t=t/(norm(t));
eigvect(:,m)=x;
eigval(m)=lambda;
teigvect(:,m)=t;
Anew=A-x*transpose(x)*A*x*transpose(x);
A=Anew;
fprintf('\n lambda= %f',lambda);
end
teigvect(:,[2 3])=teigvect(:,[3 2])