det_int = @(A) round(det(A)) ;

n = 3 ;
format rat ;
A = round( 10*rand(n,n) ) ;



% Scambio righe
fprintf( 'det_A = %d\n' , det_int(A) ) ;
A_mod = scambia_righe( A , 1 , 2 ) ;
fprintf( 'det_A dopo scambio righe = %d\n\n' , det_int(A_mod) ) ;

% Moltiplicazione riga
fprintf( 'det_A = %d\n' , det_int(A) ) ;
A_mod = moltiplica_riga( A , 1 , 10 ) ;
fprintf( 'det_A dopo moltiplicazione riga = %d\n\n' , det_int(A_mod) ) ;

% Somma di riga con multiplo di altra riga
fprintf( 'det_A = %d\n' , det_int(A) ) ;
A_mod = somma_righe( A , 1 , 10 , 2 ) ;
fprintf( 'det_A dopo somma multiplo riga = %d\n\n' , det_int(A_mod) ) ;

% Eliminazione Gauss
fprintf( 'det_A = %d\n' , det_int(A) ) ;
A_mod = eliminazione_Gauss( A ) ;
fprintf( 'det_A dopo Gauss = %d\n' , det_int(A_mod) ) ;
prodotto_diagonale = prod( diag( A_mod ) ) ;
fprintf( 'Prodotto termini diagonali = %d\n' , round(prodotto_diagonale) ) ;
format short ;