# 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] + (2*rbinom(1,1,0.5)-1)
	y[i+1] <- y[i] + (2*rbinom(1,1,0.5)-1)
}

#==============================================================================
# Trajetorias de dois passeios aleatorios
#==============================================================================
liminf = -15  # Declarar como inteiro
limsup =  15  # Declarar como inteiro
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(liminf,limsup),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")
#for(i in liminf:limsup)
#    lines(1:n,rep(i,n),lty=3,lwd=1)
#for(i in 1:n)
#    lines(rep(i,100),seq(liminf,limsup,length=100),lty=3,lwd=1)

# rm(list=ls()) 
#==============================================================================
# graphics.off()
#==============================================================================
# Fim
#==============================================================================