RFoptions(seed=0)
model <- RMexp()
x <- seq(0, 10, if (interactive()) 0.02 else 1)
plot(model, ylim=c(0,1))
plot(RFsimulate(model, x=x, seed=0))
plot(RFsimulate(RPgauss(model), x=x, seed=0), col=2) ## the same
RFoptions(seed=NA)
Run the code above in your browser using DataLab