R code for plot4
picture for plot4