K=5000; %numero de atletas simulados %
m0=5;
m1=10;
sigma=5;
m=m0;
m1=[2:1:20];
for k=1:4
   sigma=k*2;
   X=ones(K,1)*m1+sigma*randn(K,length(m1)); %concentraciones%
   p=mean(X>5+1.64*sigma);%proporción de rechazo de H0%
   AXIS([2 20 0 1]),
   subplot(2,2,k), plot(m1,p),
end

   