symulacje

Generator liniowy

x=7; a=11; c=0; m=31
y=c(x/m,(a*x+c) %% m)

#wykres rozrzutu – funkcja plot
plot(x[1],x[2],xlim=c(0,1),ylim=c(0,1))
for (i in 1:35)
{
  x=a*x+c
  # %% - operator modulo
  x=x %% m
  y=c(y,x/m)
  #zatrzymanie wykonania programu na 0,5 sek.
  Sys.sleep(.5)
  points(y[i],y[i+1])
}