%genera 60 bombas distribuidas uniformemente en 
%la 9 zonas de Londres
n=60; %nº de bombas %
k=9; %nº de zonas (categorias) %
L=500 ; %nº de simulaciones % 

K=zeros(L,1);
for l=1:L
b=unidrnd(k,n,1); %generacion de las localizaciones de las n bombas bajo H0 % 

% cuenta el numero de bombas en cada zona %
N=zeros(k,1);
for j=1:k
   for i=1:n
      if b(i)==j N(j)=N(j)+1; 
      end
   end
end

% Calculo del estadístico %
for j=1:k
   K(l)=K(l)+N(j)*N(j);
end
K(l)=-n+(k/n)*K(l);
end 

%calculo del umbral (región critica)%
u=prctile(K,95)
hist(K)
   