added script to plot SFS