% concentration
clear
close all

p = .24;
n = 9;
T = 1000;
e = .01;

figure(1)
for i=1:n
    PR(i) = 0;
    for t = 1:T
        x = (rand(2^i,1) < p);
        phat(i,t) = mean(x);
    end
    subplot(3,3,i)
    hist(phat(i,:),(0:.01:1)')
    axis([0 1 0 1000])
end

