Script for plot 1