clear;
clc;

effettuati=0;
successi=0;

while(effettuati < tot)
        %crea il sacchetto
        sacc = mod(floor(rand(1,2)*10000),2);   % 0 = bianco, 1 = nero
        sacc(1)=1;                              % impongo almeno una delle due sia nera
       
        %estraggo la prima pallina attraverso un ciclo random che "sceglie"
        % 1 o 2, valori che corrispondono alla posizione dell'array sacc
        primo_estratto=mod(floor(rand(1,1)*10000),2)+1;
        if(sacc(primo_estratto) == 1)           %se il primo estratto vale 1
            effettuati=effettuati+1;            %incremento il contatore delle prove effettuate
            if(primo_estratto==1)               %un if "brutto" per determinare qual'e' l'atra pallina da estrarre.
                secondo_estratto=2;
            else
                secondo_estratto=1;
            end
            if(sacc(secondo_estratto) == 1)     %se anche la seconda e' nera allora incrementa il contatore dei successi.
                successi=successi+1;
            end
        end
end
statistica=successi/effettuati                  %stampa brutalemnte il rapporto fra eventi favorevoli e numero di tentativi