python import numpy as np from astropy import units as u from astropy.timeseries import TimeSeries from astropy.visualization import quantity_support times = np.linspace(0, 10, 100) * u.day magnitudes = np.sin(times.value) * u.mag fluxes = 10**(-0.4 * magnitudes) ts = TimeSeries(time=times, flux=fluxes) quantity_support() ts.plot()