tcl.exp=function(n, N=200, titulo=" ", yl=c(0, .4)) { ## início da função – tcl.exp
  medias=numeric(N)
  for (i in 1:N) medias[i]= mean(rexp(n,1/3))
  hist(medias, xlim=c(-1,10), ylim=yl, freq=F, main=titulo)
  x=seq(-1,10, .02)
  points(x, dnorm(x, 3, 3*sqrt(1/n) ), type="l", lwd=3)
} ## fim da função
##graphics.off()
par(mfrow=c(3,3), mai=c(.3,.4,.1,.1))
tcl.exp(1,titulo="n=1")
tcl.exp(2,titulo="n=2")
tcl.exp(3,titulo="n=3")
tcl.exp(4,titulo="n=4")
tcl.exp(5,titulo="n=5")
tcl.exp(6,titulo="n=6")
tcl.exp(10,titulo="n=10",yl=c(0,.6))
tcl.exp(15,titulo="n=15",yl=c(0,.6))
tcl.exp(20,titulo="n=20",yl=c(0,.6))