In [1]:
# rm(list=ls())
options(OutDec = ",") 
#==============================================================================
# Exemplo para gerar de uma normal padrao pelo metodo de Box-Muller
#==============================================================================
set.seed(54345) # Semente
n  <- 100001    # Tamanho da amostra da N(0,1)
#==============================================================================
m  <- n-2*floor(n/2)
a1 <- sqrt(-2*log(runif(floor(n/2)+m)))
a2 <- 2*pi*runif(floor(n/2)+m)
x  <- c(a1*cos(a2),a1*sin(a2))[1:n]
In [2]:
#==============================================================================
# Histograma
#==============================================================================
par(mfrow=c(1,1),lwd=2,cex.lab=1.5,cex.axis=1.5,lab=c(10,5,0),
    mar=c(4.5,5,2,1),bty="n")
hist(x,nclass=50,prob=T,main="",ylim=c(0,0.40),xlim=c(-5,5),col="darkgreen",
     ylab=expression(f(x)),xlab=expression(x))
xseq <- seq(-5,5,length=1000)
yseq <- dnorm(xseq)
lines(xseq,yseq,col="red",lwd=3)
No description has been provided for this image
In [3]:
# rm(list=ls()) 
#==============================================================================
# graphics.off()
#==============================================================================
# Fim
#==============================================================================