N=10000;
h=zeros(1,N);
p=zeros(1,N);
for i=1:N
    x=randn(1,100);
    y=randn(1,100);
    [h(i),p(i)]=ttest2(x,y);
end
hist(p,30);
frakcja_odrzucona=sum(h)/N