%% BINARIS SZAMSOROZATOK HOSSZA
nRun = 10000000;
nN   = [10 100 1000];
%% futasi valtozok
nBin  = floor(5+ 4.5*log2(nN(end)));
vBins = zeros(length(nN),nBin);
bins  = 1:nBin;
bb    = 1;
%% kiserletek
for iSize=nN;
   for ii=1:nRun;
      ossSor = 2*round(rand(1,iSize))-1;
      ossSor(:,iSize+1) = - ossSor(:,iSize);
      valt   = ossSor(:,2:end)-ossSor(:,1:(end-1));
      valt   = abs(valt)/2;
      indF   = find(valt==1);
      hossz  = indF - [0 indF(1:end-1)];
      [bBin,bins] = hist(hossz,bins);
      vBins(bb,:) = vBins(bb,:) + bBin;
      if ~mod(ii,10000);
          fprintf('%03d',ii/10000);
          if ~mod(ii,100000);
              fprintf('\n');
          end;
      end;
    end;
    bb=bb+1;
    fprintf('\n');
end;
%% kirajzolas
clf; hold on; box on;
plot(bins,vBins'/nRun+eps);
xlabel('Hossz');
ylabel('Gyakorisag');
title('Alsorozatainak gyakorisága - 10m futtatás');
legend('10 hosszúságú','100 hosszúságú','1000 hosszúságú')
