% Métodos numéricos - Cap. 2.
% Resolução de sistemas pelo método iterativo 
%   de Gauss-Seidel.
% Resolução do exercicio 1 f) da ficha prática nº4.
%  Carlos Balsa 2007/2008.
%
clear all;

A=[2.44  0.40  2.04;
   1.52  2.49  0.97;
   2.44  3.96  6.33];
b=[10 -1 1.2]';

% Decomposiçao
D=diag(diag(A));
U=triu(A,1);
L=tril(A,-1);
M=L+D;
N=-U;
c=M\b;

x=zeros(3,1);
nb=norm(b,2);

tol=1e-8;
resn=1;
k=0;
while (resn>tol)
    k=k+1;
    y=x;
    x=(M\(N*x))+c
    resn=norm(x-y,2)/norm(x,2);
end


k
resn