load('seba.mat');

S_mean = mean(S);
S_var = var(S, 1);
D_mean = mean(D);
D_var = var(D, 1);

plot(normpdf(-0.1:0.001:0.2, S_mean, sqrt(S_var)));
hold on;
histogram(S);
istogram(D);
legend('S', 'D');
hold off;
title('Different Hamming distances for 1000 iterations');
xlabel('Normalized Hamming distance');
ylabel('Occurrences');
set(gca, 'fontsize', 8, 'fontname', 'Times New Roman');