symulacje

##########################################################################################
###                                                                                      Generator liniowy                                                                                     ###
##########################################################################################


### Generowanie 35 losowych wartości i prezentacja wykresu widmowego
x=7; a=11; c=0; m=31 ### parametry początkowe
y=c(x/m,(a*x+c) %% m) ### znak %% oznacza modulo (reszta z dzielenia)
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)
Sys.sleep(.5) ### spowolnienie wykonania programu, możemy obserwować kolejne punkty na wykresie
points(y[i],y[i+1])

}

### Generowanie 10000 losowych wartości i prezentacja wykresu widmowego
x=111; a=71; c=7; m=1371 ### parametry początkowe
y=c(x/m,(a*x+c) %% m) ### znak %% oznacza modulo (reszta z dzielenia)
plot(x[1],x[2],xlim=c(0,1),ylim=c(0,1))
for (i in 1:35)
{
    x=a*x+c
    x=x %% m
    y=c(y,x/m)
    points(y[i],y[i+1])
}