intervalli<-function(media=4,sqm=1,alfa=0.05,nN=100,nn=20)
{
	#Description
	#	Visualizza l'ampiezza degli intervalli di confidenza derivanti da campioni estratti
	#	da una distribuzione normale con media e varianza variabili
	#Usage
	#	intervalli(media=4,sqm=1,alfa=0.05,nN=100,nn=20)
	#Arguments
	#	media Media della variabile aleatoria normale da cui sono estratti i dati
	#	sqm Scarto quadratico medio della variabile aleatoria normale da cui sono estratti i dati
	#	alfa Errore di prima specie
	#	nN Numero dei campioni da estrarre
	#	nn Numerosità dei campioni da estrarre
	#Value
	#	Visualizza l'ampiezza dei campioni e restituisce il numero di quelli che non comprendono la media
	m.sim<-0
	for(i in 1:nN)
		{
			sim<-rnorm(nn,media,sqm)
			m.sim[i]<-mean(sim)
		}
	zeta<-qnorm(1-alfa/2,0,1)
	lim<-matrix(0,2,nN)
	for(i in 1:nN)
		{
			lim[1,i]<-m.sim[i]-zeta*sqm/sqrt(nn)
			lim[2,i]<-m.sim[i]+zeta*sqm/sqrt(nn)
		}
	y<-seq(1:nN)
	y<-matrix(y,nN,2)
	lim2<-t(lim)
	plot(y,lim2)
	for(i in 1:nN)
		lines(y[i,],lim2[i,])
	lines(c(0,100),c(media,media))
	vero<-matrix(0,nN,2)
	vero<-c(lim2[,1]>media,lim2[,2]<media)
	sum(vero)
}