Add plot1 script and figure