% Métodos numéricos - Cap. 2.
% Resolução de sistemas pelo método iterativo 
%   de Gauss-Seidel.
% Resolução do exercicio 1 d) 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]';
x=zeros(3,1);
nb=norm(b,2);

tol=1e-8;
resn=1;
k=0;
while (resn>tol)
    k=k+1;
    y=x;
    x(1,1)=(10-0.4*x(2,1)-2.04*x(3,1))/2.44;
    x(2,1)=(-1-1.52*x(1,1)-0.97*x(3,1))/2.49;
    x(3,1)=(1.2-2.44*x(1,1)-3.96*x(2,1))/6.33;
    resn=norm(x-y,2)/norm(x,2);
end
k     % numero de iteraçoes
resn  % erro final