We investigate stability of optical solitons in graded-index (GRIN) fibers by solving an effective nonlinear Schrödinger equation that includes spatial self-imaging effects through a length-dependent nonlinear parameter. We show that this equation can be reduced to the standard NLS equation for optical pulses whose dispersion length is much longer than the self-imaging period of the GRIN fiber. Numerical simulations are used to reveal that fundamental GRIN solitons as short as 100 fs can form and remain stable over distances exceeding 1 km. Higher-order solitons can also form, but they propagate stably over shorter distances. We also discuss the impact of third-order dispersion on a GRIN soliton.