% Métodos Numéricos - Engs. Civil e Electrotécnica 
% exercico 3 da ficha 2
% Carlos Balsa - 2009/2010

A=[1 -2 5 8; 7 1 -1 1; 1 3 10 -3; 1 -9 3 1];

% alínea a)
[L,U,P]=lu(A)

% alínea b)
b=[44 14 36 8]';
y=L\(P*b)
x=U\y

% alínea c)
b=[12 8 11 -4]';
y=L\(P*b)
x=U\y

