In [1]:
# rm(list=ls())
options(OutDec = ",") 
#==============================================================================
# Exemplo: passeio aleatorio
#==============================================================================
set.seed(123456)

n   <- 120        # Tamanho da sequencia
x   <- rep(0,n+1) # Primeiro caminho   
y   <- rep(0,n+1) # Segundo caminho  
for(i in 1:n){
	x[i+1] <- x[i] + rnorm(1)
	y[i+1] <- y[i] + rnorm(1)
}
In [2]:
#==============================================================================
# Trajetorias de dois passeios aleatorios
#==============================================================================
par(mfrow=c(1,1),lwd=2.0,cex.lab=1.5,cex.axis=1.5,lab=c(8,8,5),
    mar=c(4.5,5,1,1),cex.main=2.0,bty="n")
plot(0:n,x,t="l",col="red",lwd=2,xlab=expression(n),ylab=expression(X[n]),
    ylim=c(-15,15),xlim=c(0,120))
points(0:n,x,pch=16,cex=1,col="red")
lines(0:n,y,col="blue",lwd=2)
points(0:n,y,pch=15,cex=1,col="blue")
No description has been provided for this image
In [3]:
# rm(list=ls()) 
#==============================================================================
# graphics.off()
#==============================================================================
# Fim
#==============================================================================