% i due esempi di sistemi lineari della lezione
A = [2 1 -1 -3; 4 2 1 -3; 2 -1 3 5; 3 -1 -1 1]
b = [-5 -7 11 8]'
b2 = [-5 -7 11 9]'
% risolvere il sistema lineare  Ax = b
x = A\b
% la nostra soluzione a mano
sol = [1 -6 1 0]' + 13/6 * [0 2 -1 1]'
rank(A)
[A, b]
rank([A, b])
% soluzione di  A * x = 0
null(A)
v = [0 2 -1 1]';
norm(v),  sqrt(6)
v/norm(v)
det(A)
cond(A)
% fattorizzazione  P * A = L * U
[L U P] = lu(A)
sol(A,b)
help sol
x2 = A\b2
rank([A, b2])
quit
# Octave 3.6.4, Wed Dec 11 17:02:06 2013 W. Europe Standard Time <unknown@achilles-TOSH>
