# num_dices - number of dices, N - number of times we are throwing them
clt <- function (N, num_dices)
{
	r <- rep (0, N);
	for (i in 1:num_dices)
	{
		r <- r + sample (1:6, N, replace=T);
	}
	r <- r/num_dices;
	hist (r, 100, freq=F, main=num_dices);
	x <- seq(min(r), max(r), length=100);
	lines (x, Norm(x, sd(r), mean(r)), type="l");
}

for (n in c(1,10,100,1000,10000))
{
	clt (1000, n);
	Sys.sleep (3); 
}