clear; %clear variables

xmean=100;
xsigma=10;

Ns=5:5:100; %number of samples for calculating mean

kk=1;
for N=____
 jj=1;
 for samples=1:25 %number of means to calculate stdev from
  for ii=1:N
   a(ii)=mygaussian(_______, ______);
  end
  m(______)=mymean(a);
  jj=jj+1;
 end
 s(kk)=mystdev(m);
 kk=kk+1;
end
plot(Ns,______,'b;samples;');
hold on; 
plot(Ns,______,'r;sigma/sqrt(N);');
hold off;
