commented out information to use source("Zhang13.R") for a quick simulation