source("psim.R")
set.seed(11333)
data <- rnorm(30,0,1)

kmax <- 10
vp <- numeric(kmax)
for (k in 1:kmax) {
    vp[k] <- psim(data,0)
    cat(sprintf("psim(data,0)=%g\n",vp[k]))
}
cat(sprintf("Estimated p=%g\n",mean(vp)))