% drzewa
srednice=[8.5 7.6 9.3 5.5 11.4 6.9 6.5 12.9 8.7 4.8 4.2 8.1 6.5 5.8 6.7 2.4 11.1 7.1 8.8 7.2];
boxplot(srednice)
sigma=std(srednice);

me=mean(srednice);

Nboot=100000;
A=zeros(1,Nboot);
B=zeros(1,Nboot);
for i=1:Nboot
wartosci=randsample(srednice,20,true);
%A(i)=sum(wartosci)/20;
A(i)=mean(wartosci);
B(i)=median(wartosci);
end 
[n,x]=hist(A,30);
bar(x,n);
hi_me=prctile(A,95);
lo_me=prctile(A,5);

disp(sprintf('srednia i 90 procentowy przedzial %.3f %.3f %.3f\n',lo_me,me,hi_me))


hi_md=prctile(B,95);
lo_md=prctile(B,5);
md=median(srednice)

disp(sprintf('mediana i 90 procentowy przedzial %.3f %.3f %.3f\n',lo_md,md,hi_md))
hist([A' B'],sort(srednice))