>>> bcpd = BCPD(max_tcp=2, max_scp=1, max_harmonic_order =10, random_seed=1, max_iter=10000)
>>> plot_change_points(data=df, cp_object=bcpd, cp_style="axvline")