/* PROGRAM TO DISPLAY THE POWER OF CENTRAL LIMIT THEOREM */

loops = 500;
nobs = 50;
normmean = zeros(loops,1);
counter = 1;

do while counter <= loops;
	x = (rndn(nobs,1))^2;
	meanx = meanc(x);
	normmean[counter,1] = sqrt(nobs)*(meanx - 1)/sqrt(2);
	counter = counter +1;
endo;

library pgraph;
graphset;
title("Simulated Sampling Distribution of Transformed Mean from Chi-Square(1) pdf"); 
{b,m,f} = hist(normmean,20);