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