  We address stability of multi-solitons in the cubic NLS (nonlinear
Schr\"{o}dinger) equation on the line. By using the dressing transformation and
the inverse scattering transform methods, we obtain the orbital stability of
multi-solitons in the $L^2(\mathbb{R})$ space when the initial data is in a
weighted $L^2(\mathbb{R})$ space.
