function B=minor(A,i,j) 
%restituisce la matrice minore di A rispetto a i,j
%poi determinante(A)=sum ...
%(-1)^(j+1)*A(1,j)*determinante(minor(A,1,j))
X=A(1:i-1,1:j-1);
Y=A(1:i-1,j+1:end);
Z=A(i+1:end,1:j-1);
W=A(i+1:end,j+1:end);
B=[X Y; Z W];

