Ex 5.24

Simulering för sannolikheten att födelsedagar sammanfaller

clear;
n = 25;
N = 10000;

antal_dubletter = 0;
for i = 1:N
    bday = ceil(rand(n,1)*356);
    sbday = sort(bday);
    for j = 2:n
        if sbday(j) == sbday(j-1)
            antal_dubletter = antal_dubletter + 1;
            break;
        end
    end
end

sannolikhet = antal_dubletter / N;
disp(['Sannolikheten är: ',num2str(sannolikhet*100),'%']);

Sannolikheten är: 58.37%