%karma dla zwierzat
no_zwierzat=12;
time=4;
przyrost=[15.43 16.92 14.43 12.94 15.92 17.42 18.91 16.92 14.93 14.49 15.92 15.43];
me=mean(przyrost)
przyrost_mean=15.80;
%hist(przyrost)
sigma=std(przyrost);

Nboot=100000;
A=zeros(1,Nboot);
for i=1:Nboot
wartosci=randsample(przyrost,12,true);
A(i)=sum(wartosci)/12;
end 
[n,x]=hist(A,30);
bar(x,n)
p=sum(A<15)/Nboot
p2=tcdf((15-me)*sqrt(12)/sigma,11)
p3=normcdf((15-me)*sqrt(12)/sigma)