N = 100000;
n = 1:N;
xbar = zeros(N,1);
p = 0.3;
for i=1:N
	xbar(i) = binornd(n(i),p)/n(i);
end
subplot(1,1,1)
plot(n,xbar);
