% Plots of CLT for various distributions

r = 1000000;
n = 10;
mu = 0;
sigma = 1;
s = zeros(r,1);

% z = randn(n,r);                      % standard normal distribution
% z = sqrt(12)*(rand(r,n) - .5);       % uniform distribution 
z = sqrt(18)*(sqrt(rand(r,n)) - 2/3);  % triangular distribution
x = mu + sigma*z;
s = sqrt(n)*(mean(x,2) - mu);

hist(s,100);

return