<br>
limite.central<-function(r=runif, distpar=c(0,1), m=.5, s=1/sqrt(12), n=c(1,3,10,30), N=1000) {<br>
  for(i in n){<br>
    if(length(distpar)==2){<br>
      x <- matrix(r(i*N, distpar[1], distpar[2]), nc=i)<br>
    }else{<br>
      x <- matrix(r(i*N, distpar), nc=i)<br>
    }<br>
    <br>
    x<-(apply(x, 1, sum) - i*m )/(sqrt(i)*s)<br>
    hist(x, col='light blue', probability=T,main=paste("n=",i),<br>
    ylim=c(0,max(.4, density(x)$y)))<br>
    lines(density(x), col='red', lwd=3)<br>
    curve(dnorm(x), col='blue', lwd=3, lty=3, add=T)<br>
  <br>
    if( N>100 ){<br>
      rug(sample(x,100))<br>
    }else{<br>
      rug(x)<br>
    }<br>
  }<br>
}