% Métodos Numéricos - Carlos Balsa - 2008/09
% Resolução do ex. 1 da ficha 11

A=[0.0406 2.02 1.27; -51.2 0.0702 -17.9; 11.1 -15.7 30.3];

b=[-7.49 7.51 12.8]';

% alinea a)
dA=det(A);

% alinea b)
[L, U, P]=lu(A);

% alinea c)
y=L\(P*b);
x=U\y;

% alinea d)
c=[3.33 -69.33 25.7]';
y=L\(P*c);
x=U\y;

% alinea e)
r=c-A*x;
nr=norm(r,2)