Add plot1

Code and figure for plot1