N=1000; ilosc=0;
for proba=1:N
    urodziny=ceil(365*(rand(20,1))); % generujemy 20 dat urodzin
    koincydencja=0;osoba=1;          % i sprawdzamy, czy która¶ siê powtarza
    while osoba<20-1 & koincydencja==0 
        osoba=osoba+1;
        for druga_osoba = osoba+1:20
            if urodziny(druga_osoba)==urodziny(osoba)
                koincydencja=1;break;
            end
        end
    end
    ilosc=ilosc+koincydencja;
end
ilosc   %wypisujemy ilo¶æ przypadków koincydencji 
