# The normalized spectrum is simply the old spectrum devided by the
# fitted continuum, which returns a new object:
#
spec_normalized = spectrum / y_continuum_fitted
#
f, ax = plt.subplots()  # doctest: +IGNORE_OUTPUT
ax.plot(spec_normalized.spectral_axis, spec_normalized.flux)  # doctest: +IGNORE_OUTPUT
ax.set_title("Continuum normalized spectrum")  # doctest: +IGNORE_OUTPUT
ax.grid(True)  # doctest: +IGNORE_OUTPUT