# to see contents of create_data function, look at the top-most algorithm's code
figure, axes, handles = create_plots(data, baselines)
for i, (ax, y) in enumerate(zip(axes, data)):
    if i == 3:
        lam = 1e2
    else:
        lam = 1e3
    baseline, params = baseline_fitter.pspline_drpls(y, lam=lam)
    ax.plot(baseline, 'g--')