General (independent of models selected for surfactant molecular aggregates) analytical relations are derived to describe the initial stage of slow relaxation in micellar solutions with spherical micelles. This stage precedes the final stage of the relaxation occurring via an exponential decay of disturbances with time. The relations obtained are applicable throughout the interval of micellar solution concentrations from the first to the second critical micellization concentration. It is shown that the initial stage is characterized by power laws of variations in the concentrations of monomers and micelles with time, these laws being different for the relaxation processes proceeding from above and below toward equilibrium values of micellar solution parameters. Relations are derived for the duration of this stage, and the effect of initial conditions is studied. Characteristic times of the power-law stage are determined and compared with the characteristic time of the final exponent-law relaxation stage. The behavior of these times is investigated at surfactant solution concentrations in the vicinity of, and noticeably above, the first critical micellization concentration. On the basis of the droplet and quasi-droplet thermodynamic models of surfactant molecular aggregates, numerical solutions are found for nonlinearized equations of slow relaxation for the time dependence of surfactant monomer concentrations at all stages of the slow relaxation. Numerical results obtained from the models are compared with the results of a general analytical study.