Eratosthenes <- function(n) {
if(n >= 2){
sieve <- seq(2, n)
primes <- c()
for (i in seq(2, n)) {
if (any(sieve == i)) {
primes <- c(primes,i)
sieve <- c(sieve[(sieve %% i) != 0], i)
}
}
for (i in primes){
if (primes[i+1] == primes[i]+2 ) {
twinprime.list <- list(primes[i],primes[i+1])
}
}
return (twinprime.list)
} else {
stop("Input value of n should be at least 2.")
}
}
Eratosthenes(1000)