% Métodos numéricos - Carlos Balsa 2007/2008.
% Resolução do exercicio 2 da ficha prática nº11.
%
clear all;

A=[51.2 0.0702 -17.9; 0.0702 2.02 1.27; -17.9 1.27 30.3];
b=[1 1 1]';


% alinea a)
dA=det(A)


% alinea d)

x=zeros(3,1);
tol=1e-5;
erro=1;

k=0;
while (erro>tol)
    k=k+1;
    y=x;
    x(1,1)=(1-0.0702*x(2,1)+17.9*x(3,1))/51.1;
    x(2,1)=(1-0.0702*x(1,1)-1.27*x(3,1))/2.02;
    x(3,1)=(1+17.9*x(1,1)-1.27*x(2,1))/30.3;
    erro=norm(x-y,2);
    
end
k     % numero de iteraçoes
erro  % erro final
x     % solução