# Reuse yearly and prior from the previous example
components = yearly.sample_curve(prior, sum=False)

fig, _ = yearly.plot_decomposition(components)
fig.suptitle(
    "Yearly Fourier seasonality decomposed into basis functions"
)
plt.show()