/* ESTIMATED SAMPLING DISTRIBUTION FOR THE MEAN OF A N(0,1) R.V. */

nobs = 3;
z = zeros(500,1);
 
for ii(1,500,1);
  x = rndn(nobs,1);
  xbar = meanc(x);
  z[ii,1] = xbar;
endfor;

library pgraph;
graphset;
xtics(-2,2,0.2,1);
Title("Estimated Sampling Distribution (n=3) for the Mean of a N(0,1) R.V.");
{b,m,f} = hist(z,20);