% Métodos numéricos - Cap. 2
% resolução de sistemas pelo método de Cholesky
% Resolução do exercicio 1 b) e c) da ficha prática nº3
%  Carlos Balsa 2007/2008
%
A=[0.96 1.52 2.44; 1.52 2.49 3.96 ; 2.44 3.96 6.33];
% Factorizaçao de Cholesky
R=chol(A);
L=R';
% Resoluçao do sistema da alinea b)
disp('Soluçao sistema alinea b)');
b=[4.92 7.97 12.73]';
y=L\b;
x=L'\y

% Resoluçao do sistema da aliea c)
disp('Soluçao sistema alinea c)');
b=[1.00 1.57 2.51]';
y=L\b;
x=L'\y