%
%  Power Method
%
A=[6 5 -5;2 6 -2;2 5 -1];
x=[-1;1;1];
X(1,:)=x';
format long
m=30;
%
for k=1:m
 y=A*x;
 r=phi(y)/phi(x);
 R(k)=r;
 x=y/norm(y);
 X(k+1,:)=x';
end;
