p1=p_cut1;   %inizializzo la probabilità (1)
p2=p_cut2;  %inizializzo la probabilità (1)
%ciclo per tutte le tabelle di contingenza equivalenti
for n=1:size(M,3)
    a=M(1,1,n);     %calcolo elemeti (solo per maggior chiarezza ei passaggi successivi)
    b=M(1,2,n);
    c=M(2,1,n);
    d=M(2,2,n);
    
    test1=nchoosek(a+b,a)*nchoosek(c+d,c)/nchoosek(a+b+c+d,a+c)    %calcolo la prob per la tabella selezionata
    test2=pdf('Hypergeometric',a,a+b+c+d,a+b,a+c)
    
    if test1<=p_cut1  %verifico se la prob della tabella è più estrema o ugualmente probabile
        p1=p1+test1;
    end
    
    if test2<=p_cut2
    p2=p2+test2;
    end
end
p1
p2