Added script for plot 4