Hi all, this is an assignment for the birthday paradox simulation. The code below is working fine but I think we are meant to turn parts of the code into user-defined functions and call the functions instead of having it all written out. In the GUI simulation the user inputs both N for realisations and k1 for maximum group size. The functions we are meant to have separate are randi() and checking the random numbers to see if there is a match (ind = find(rdie == b);). The problem is I don't know how to create a separate function file for randi() and get the user inputted k1 value into that function file. I only know how to create the random number from within these loops. I tried creating a function called randomnumber but I need a k value that only comes from the GUI input. Can somebody point me in the right direction? Any help is greatly appreciated
tally = zeros(1, 365);
N = 1000;
k1 = 365;
for k = 2:k1
for m = 1:N
rdie = randi(365, 1, k);
for i = 1:length(rdie)
b = rdie(i);
ind = find(rdie == b);
if length(ind) > 1
tally(k) = tally(k) + 1;
break
end
end
end
end
probability = tally/N;
x = 2:k1;
y = probability;
L = min(length(x), length(y));
plot(x(1:L), y(1:L))
xlabel('Number of people')
ylabel('Probability')
