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

A=[7 1 1; 1 5 1; 1 2 9];
b=[20 7 54]';
x=zeros(3,1);
nb=norm(b,2);

tol=1e-3;
resn=1;
k=0;
while (resn>tol)
    k=k+1;
    x(1,1)=(20-x(2,1)-x(3,1))/7;
    x(2,1)=(7-x(1,1)-x(3,1))/5;
    x(3,1)=(54-x(1,1)-2*x(2,1))/9;
    resn=norm(b-A*x,2)/nb;
end
k
resn