diff --git a/demos/sample_means/generate_plots.R b/demos/sample_means/generate_plots.R
new file mode 100644
index 0000000..eb3d59d
--- /dev/null
+++ b/demos/sample_means/generate_plots.R
@@ -0,0 +1,11 @@
+arg = commandArgs(TRUE)
+
+samples = rep(NA, 100000)
+for ( i in 1:100000 ){ samples[i] = mean(rexp(40, 0.2)) }
+
+jpeg(paste('plots/', arg, '.jpg', sep=""))
+
+hist(samples, main="", prob=T)
+lines(density(samples), col="darkblue", lwd=3)
+
+dev.off()
diff --git a/demos/sample_means/slurm_submit_generate_plots.sh b/demos/sample_means/slurm_submit_generate_plots.sh
new file mode 100644
index 0000000..feb5e40
--- /dev/null
+++ b/demos/sample_means/slurm_submit_generate_plots.sh
@@ -0,0 +1,9 @@
+#!/bin/bash
+#SBATCH --nodes=1
+#SBATCH --ntasks-per-node=1
+
+cd $SLURM_SUBMIT_DIR
+
+mkdir plots/
+
+R --vanilla -f generate.R --args "plot$SLURM_ARRAY_TASK_ID"