a = 1
c = 1
M = 120
array <- c()
for(i in 1:1000){
if(i == 1){
array[i] = 0.5 #numero entre 0 e 1
}
else{
x <- ((a * array[i - 1] +c) %% M)
array[i] = x
}
}
array = array/M
hist(array)
array2 <- c()
for (j in 1:1000){
x2 <- -(log(1 -array[j]))/2
array2[j] = x2
}
hist(array2)
q3 <- rbinom(1000,10,0.2)
hist(q3)