Shell

R code to plot this figure