Added N=%d to main
Added a=commandArgs(trailingOnly=T)
Revised two lines of comments
Added pch=3 for median plots