code for plot 4