interesting.m

fit programme for a serials of data

%An interesting phenomena according to the text book
A=zeros(6,6);
n=1;
for i=1:6
   for j=1:i
       A(i,j)=n;
       A(j,i)=n;
       n=n+1;
   end
end
[L,U]=ltriangsec(A)
b=rand(6,10);
X=zeros(6,10);
for k=1:10
    disp(sprintf('b(%d)=',k))
    b(:,k)
    disp(sprintf('X(%d)=',k))
    X(:,k)=solute(A,b(:,k));
    X(:,k)
    if k==10
        break;
    end
    b(:,k+1)=X(:,k)/norm(X(:,k),inf);
end