library("doParallel")
library("foreach")
cluster <- makeCluster(7)
registerDoParallel(cluster)

r <- 1
v <- c(1,1)
inA <- function(x) sum((x-v)^2)<r^2
N <- 100000; K <- 1000
C <- foreach(k=1:K,.combine='c') %dopar% {
	X <- matrix(rnorm(2*N),N)
	sum(as.integer(apply(X,1,inA)))
}
cat(sprintf("P(A)=%g\n",sum(C)/N/K))

stopCluster(cluster)