limite.centrale<-function(x,n=100,min=0,max=1)
{
	#Description
	#	Simula graficamente il teorema del limite centrale di dati provenienti da una distribuzione uniforme
	#Usage
	#	limite.centrale(x,n=100,min=0,max=1)
	#Arguments
	#	x Numero dei campioni da estrarre
	#	n Numerosità di ogni campione estratto
	#	min Valore minimo della distribuzione uniforme
	#	max Valore max della distribuzione uniforme
	#Value
	#	Vengono restituiti una serie di grafici che simulano il teorema del limite centrale
	#Note
	#	E' necessario aver caricato il pacchetto stepfun
	dd<-0
	for(i in 1:x)
		dd[i]<-mean(runif(n,min,max))
	dd<-(dd-(min+max)/2)/((max-min)/sqrt(12*n))
	aaa<-hist(dd,plot=F)$breaks
	par(mfrow=c(2,2))
	hist(dd,main="Istogramma delle frequenze assolute")
	hist(dd,prob=T,main="Istogramma delle frequenza relative")
	irq<-summary(dd)[5]-summary(dd)[1]
	lines(density(dd,width=irq))
	hist(dd,prob=T,main="Distribuzione media campionaria")
	bbb<-seq(min(aaa),max(aaa),0.01)
	lines(bbb,dnorm(bbb,0,1))
	plot(ecdf(dd),verticals=T,do.p=F,main="Fz distr. emp. e rip. N(0,1)")
	lines(sort(dd),pnorm(sort(dd),mean=0,sd=1))
	par(mfrow=c(1,1))
}