fixed thresholding; added benchmark functions for simulations
argle
added plotting function for penalty, plot_scad