import matplotlib.pyplot as plt
import numpy as np
from skued import spectrum_colors

t = np.linspace(0, 10, num = 1000)
y = np.exp(-t/0.4) + 10 *(1 - np.exp(-t/2)) + np.sin(2*np.pi*t)

fig, ax = plt.subplots(1,1)
ax.scatter(t, y, c = list(spectrum_colors(y.size)))

ax.set_title('Some realistic dynamics')
ax.set_xlabel('Time (ps)')
ax.set_ylabel('Intensity (a.u.)')

plt.show()