% wybory
x=ones(1,1500);
x(841:1500)=0;
Nboot=100000;
A=zeros(1,Nboot);

for i=1:Nboot
%ankieta=x(ceil(100*rand(1,1500)));
ankieta=randsample(x,1500,true);
A(i)=sum(ankieta)/1500;
end

lo=prctile(A,2.5);
hi=prctile(A,97.5);
disp(sprintf('przedzial ufnosci: %.3f - %.3f\n',lo,hi))
[n,x]=hist(A,30);
bar(x,n)
p=sum(A>=840/1500)/Nboot
line([lo lo] , [0 max(n)] ,'Color',[1 0 0] )
line([hi hi],[0 max(n)],'Color',[1 0 0] )