function sameday_prob = birthday(n)

% PEng 5.28 (5.24)
% Returnerar sannolikheten att 2+ personer av n har samma födelsedag

% Gör ntries tester där n födelsedagar slumpas fram
% Räkna hur många av dessa som har 2+ likadana födelsedagar och lagra i 
% nsame
% sameday_prob ? nsame/ntries 

% NUMMEASPEKT: BÄTTRE JU STÖRRE ntries ÄR. HUR STOR MÅSTE ntries VARA? 
% HUR FORT AVTAR FELET?

ntries = 5000;
nsame = 0;
for i = 1:ntries
		
	random_birthdays = ceil(365*rand(1, n));
	% Antal UNIKA födelsedagar i slumpvektorn
	nunique = length(unique(random_birthdays));
	
	if (n > nunique)
		nsame = nsame + 1;
	end
	
end

sameday_prob = nsame/ntries;
