T<-1; r=0.00; X<-2500; p<-0.004; seed<-250; v<-0.25; mean<-0; simulate<-function(seed) { samples<-rnorm(seed,sd=v/sqrt(seed),mean=mean/(seed)); return (S*cumprod(exp(samples))); } plot(1:seed,simulate(seed),ylim=c(1000,5000),xlab="Dzień",ylab="Cena"); for(i in 1:100) { lines(1:seed,simulate(seed),col=rgb(runif(1),runif(1),runif(1))); }