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

task_id <- as.numeric(Sys.getenv("SLURM_ARRAY_TASK_ID","1"))
cat(sprintf("task_id=%d\n",task_id))
set.seed(33111+task_id*12345)

p <- psim(data,0)
cat(sprintf("psim(data,0)=%g\n",p))
saveRDS(p,sprintf("output/p%04d.rds",task_id))