%Ogrodnik eksperymentuje z nowym gatunkiem drzew. 
%Posadzi³ 20 sztuk i po dwóch latach zmierzy³ nastêpuj±ce
%¶rednice pni (w cm):
d=[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];

% Proszê obejrzeæ boxplot.
boxplot(d,'notch','on')
% Proszê znalesc srednia srednice i 
% 90\% przedzia³ ufnosci dla sredniej.
N=length(d);
N_rep=1e4;
sr=zeros(1,N_rep);
for i=1:N_rep
    x=randsample(d,N,true);
    sr(i)=mean(x);
    me(i)=median(x);
end
mi=mean(d);
gr=prctile(sr,[5 95]);
disp(sprintf('srednia i jej przedzial ufnosci[%.2f %.2f %.2f]',gr(1),mi,gr(2)));

% Proszê znale¼æ medianê i 90\% przedzia³ ufno¶ci dla mediany.
med=median(d);
gr=prctile(me,[5 95]);
disp(sprintf('mediana i jej przedzial ufnosci[%.2f %.2f %.2f]',gr(1),med,gr(2)));

% Wynik zilustrowaæ przy pomocy histogramu.
subplot(211)
hist(sr,30)

title('srednia')
subplot(211)
hist(me,sort(d))
title('mediana')