symulacje

##########################################################################################
###                                                              Procedura realizująca metodę środków kwadratów                                                               ###
##########################################################################################


### Ustalenie wartości początkowej
x=7645372856 ### warto zmienić wartość początkową i ponownie uruchomić kod
y=c()
for (i in 1:10) ### zamiast 10 liczb możemy zadać wększą liczbę generowanych wartości
{
 x=x^2
  # otrzymujemy ciąg znaków  
 x=substr(x,6,15)
  #zmiana formatu na liczbowy
 x=as.numeric(x)
 y=c(y,x/10^10)
}
y ### wypisanie wygenerowanych wartości losowych

### histogram wygenerowanych wartości
hist(y)

### Generowanie 1000 wartości losowych
x=22483456290 ### warto zmienić wartość początkową i ponownie uruchomić kod
y=c()
for (i in 1:1000) ### zamiast 10 liczb możemy zadać wększą liczbę generowanych wartości
{
    x=x^2
    x=substr(x,6,15)
     x=as.numeric(x)
    y=c(y,x/10^10)
}
hist(y)