clear;
p = [0.0001 0.0135 0.1800 0.5040 0.3024]';
N = 10;
k = 5;
alpha = 0.05;

c = [];
for i = 1:N
  x    = floor(rndLCG(k)*10);
  c(i) = length(unique(x));
end

h = hist(c, 1:1:k);
h = [h /sum(h)]';

V    = sum( (p-h).^2 ./ p )
pval = 1 - chi2cdf(V, N-1)

if (pval > alpha)
  printf("Δεν απορρίπτεται η Η_0\n\n");
else
  printf("Απορρίπτεται η Η_0\n\n");
end