%%% Simulated Card Experiment
clc
format compact
N=100000;
M=52;
subs=13;
results=zeros(1,M);
for k=[1:N]
    n=randperm(M);
    m=ceil(n*subs/M);
    finds=zeros(1,subs);
    for j=[1:subs]
        finds(j)=min(find(m==j));
    end
    maxfind=max(finds);
    [k, M-maxfind+1];
    K = M-maxfind+1;
    results(K)=results(K)+1;
end
[1:M;results]'
bar(results)