%%reklamy

n=12;
p=2/3;
k=10;

%czy prawda
b=binopdf(k,n,p)
% czy wiecej niz 10
b=zeros(1);
for k=10:12
    b= b + binopdf(k,n,p);
end
disp(b)

%symulacja
p0=1-binocdf(9,n,p);
N=1000000;
y=random('binomial',n,p,1,N);
hist(y,[1:n]);
p1=sum(y>=k)/N;
disp(['ze wzoru ',num2str(b),', a z symulacji ',num2str(p0)])