R=sum(m,2)  %calcolo la somma delle righe della tabella di contingenza 
C=sum(m)    %calcolo la somma delle colonne della tabella di contingenza
n_matrix=1; %inizializzo indice matrice delle tabelle di contingenza equivalenti

%ciclo per calcolare tutte le possibili tabelle di contingenza equivalenti
for a=1:R(1)    
    b=R(1)-a;   %individuo elementi prima riga in modo che la loro somma non cambi
    for c=1:R(2);
        d=R(2)-c;   %individuo elementi prima riga in modo che la loro somma non cambi
        if a+c==C(1) && b+d==C(2) && a~=m(1,1) && d~=m(2,2)    %verifico che la somma delle colonne non cambi\non ripetere la tabella iniziale
           M(:,:,n_matrix)=[a b; c d];
           n_matrix=n_matrix+1;
         end
    end
end
n_matrix