s = RandStream('mt19937ar','Seed',1);
RandStream.setGlobalStream(s);
close all
sample=rand(50,1);
thetamle=max(sample)
%
maxes=[];
B=20000;
for i=1:B
bsample=bootsample(sample);
maxes=[maxes max(bsample)];
end
 
 
%parametric
maxes1=[];
B=20000;
for i=1:B
bsample=thetamle * rand(50,1);
maxes1=[maxes1 max(bsample)];
end
 
bihist(maxes', maxes1', 30)



