A=[0  -4  -8  -2
  -4  -7  -7  -8
  -9  -5  -4  -5
   0  -5  -9  -6];
n=length(A); x=(1:n)'; b=A*x;
x0=x-5;                                   % Starting vector
[s,W]=TEA1(A,b,x0,4); te=W(:,2:2:8)
q=b-A*x0; G=eye(size(A))-A;
[Y,X]=TEA(G,b,x0,q,n); Y