# to see contents of create_plots function, look at the top-most algorithm's code
baseline, params = baseline_fitter.pspline_psalsa(y, lam=(1e3, 5e0), k=0.5)
create_plots(y, baseline)