@@ -58,18 +58,18 @@ the `~astropy.modeling.Model` to be used in the fitting routine.
5858.. code-block :: python
5959
6060 >> > # Create the initial guesses
61- >> > fg1 = models.Gaussian1D(amplitude = 1900 * u.mJy, mean = 0.5 * u.um, stddev = 0.02 * u.um, \ \
62- bounds = {' amplitude' : (1700 * u.mJy, 2400 * u.mJy), \ \
63- ' mean' : (0.3 * u.um, 0.6 * u.um), \ \
64- ' stddev' : (0.001 * u.um, 0.03 * u.um)})
65- >> > fg2 = models.Gaussian1D(amplitude = 400 * u.mJy, mean = 0.63 * u.um, stddev = 0.03 * u.um, \ \
66- bounds = {' amplitude' : (200 * u.mJy, 600 * u.mJy), \ \
67- ' mean' : (0.6 * u.um, 0.7 * u.um), \ \
68- ' stddev' : (0.001 * u.um, 0.03 * u.um)})
69- >> > fg3 = models.Gaussian1D(amplitude = - 500 * u.mJy, mean = 0.79 * u.um, stddev = 0.02 * u.um, \ \
70- bounds = {' amplitude' : (- 600 * u.mJy, - 300 * u.mJy), \ \
71- ' mean' : (0.7 * u.um, 0.82 * u.um), \ \
72- ' stddev' : (0.001 * u.um, 0.03 * u.um)})
61+ >> > fg1 = models.Gaussian1D(amplitude = 1900 * u.mJy, mean = 0.5 * u.um, stddev = 0.02 * u.um,
62+ ... bounds = {' amplitude' : (1700 * u.mJy, 2400 * u.mJy),
63+ ... ' mean' : (0.3 * u.um, 0.6 * u.um),
64+ ... ' stddev' : (0.001 * u.um, 0.03 * u.um)})
65+ >> > fg2 = models.Gaussian1D(amplitude = 400 * u.mJy, mean = 0.63 * u.um, stddev = 0.03 * u.um,
66+ ... bounds = {' amplitude' : (200 * u.mJy, 600 * u.mJy),
67+ ... ' mean' : (0.6 * u.um, 0.7 * u.um),
68+ ... ' stddev' : (0.001 * u.um, 0.03 * u.um)})
69+ >> > fg3 = models.Gaussian1D(amplitude = - 500 * u.mJy, mean = 0.79 * u.um, stddev = 0.02 * u.um,
70+ ... bounds = {' amplitude' : (- 600 * u.mJy, - 300 * u.mJy),
71+ ... ' mean' : (0.7 * u.um, 0.82 * u.um),
72+ ... ' stddev' : (0.001 * u.um, 0.03 * u.um)})
7373
7474 Now comes the actual fitting:
7575
0 commit comments