% Matrice e termine noto particolari
A = [ 0 1 2 ; 4 5 6 ; 7 7 8 ] ;
b = [ 1 ; 1 ; 3 ] ;

% Eliminazione di Gauss su concatenazione di A e b
format rat ;
Tf = eliminazione_Gauss( [ A b ] ) ;
T  = Tf(:,1:3) ; % estrazione matrice triang. sup. T 
f  = Tf(:,  4) ; % estrazione termini noti ridotti f

% Sostituzione all'indietro
x = sostituzione( T , f ) ;

% Verifica
fprintf('Ax - b =\n') ;
disp( A*x-b ) ;

fprintf('x - A\\b =\n') ;
disp( x - A\b ) ;

% Determinazione inversa
b = eye(3,3) ;

% Eliminazione di Gauss su concatenazione di A e b
Tf = eliminazione_Gauss( [ A b ] ) ;
T  = Tf(:,1:3) ; % estrazione matrice triang. sup. T 
f  = Tf(:,4:6) ; % estrazione termini noti ridotti f

% Sostituzione all'indietro
inv_A = sostituzione( T , f ) ;

% Verifica
fprintf('A * inv_A =\n') ;
disp( A * inv_A ) ;

fprintf('inv_A * A =\n') ;
disp( inv_A * A ) ;

fprintf('inv_A - inv(A) =\n') ;
disp( inv_A - inv(A) ) ;
format short ;