#test de rnorm
rm(list=ls())
Nt<-10000 #Nombre de pas en temps
M<- 8 #On va afficher le résultat de 8 expériences
nu<- 0.2 #Parametre de Drift
D<-0.1 #Coefficient de diffusion
Dt=min( c( 1e-2/abs(nu) , (1e-4/(2*D)) ) ) #Pas en temps

colors = c("red", "yellow", "green", "violet", "orange", "blue", "pink", "cyan") 
VM=nu*Dt*Nt #Valeur moyenne
ET=sqrt(2*D*Dt*Nt) #Ecart type
yborne<-c( -abs(VM)-2*ET , abs(VM)+2*ET )


for(iexp in 1:M)
{
  X=c(0) #Va contenir les positions successives de la particules
  for(i in 2:Nt)
  {
    X=c( X , X[i-1]+nu*Dt+sqrt(2*D)*Dt*rnorm(1,0,1/sqrt(Dt)) )
  }
  plot(X,type="l",ylim=yborne,col=colors[iexp]);par(new=TRUE);
}
par(new=FALSE)