Add plot4

Code and figure for plot3