Напишите код для повтора в функции:
run_function <- function(){
n=1
success = FALSE
while(!success){
n=n+1
if(rbinom(1, n, (n*n+1)^-1) == 1)
success = TRUE
}
n
}
Используйте replicate
, чтобы повторить его 1000 раз, и сохраните вектор результатов в x
.
x <- replicate(1000, run_function())