## Wrap the main code into a function
lemniscate.run <- function(N = 1000, x = runif(N), y = runif(N)) {
  logic <- (x^2+y^2)^2 <= x^2-y^2
  K <- length(which(logic))
  return(K/N*4);
}

## Do 5 runs
results = rep(0,5)
for(i in 1:5) {
  results[i] = lemniscate.run();
}

## Summarize results
print(sprintf("The mean of 5 runs is %f with standard deviation of %f",
              mean(results),
              sd(results)));