Add plot1.R and output