source("sim4.R")
system("vi sim4.R")
source("sim4.R")
