A novel method is reported for calculating the nonlinear resonant change in refractive index of a doped saturable medium at the short-pulse excitation. The derived system of equations for energy density and phase of a short pulse propagating in the doped medium allows one to address a rather general case of the resonant index non-linearity caused by the dope's population perturbation (bleaching) effect. As an example, the energy-dependent changes in the refractive index and transmission coefficient of a Cr 4+ :YAG crystal under the ns-range (wavelength 1.06 µm) excitation are modeled. It is known [1] that propagation of a short pulse through a resonantly amplifying (or absorbing) medium with a "slow" relaxation of the excited state is addressed by the Avizonis-Grotbeck's equation:where ε = E p s σ GSA hυ is the ratio of incident pulse energy density E p /s (where E p and s are the pulse energy and laser beam geometrical cross-section) and saturating energy density hυ/σ GSA (where hυ is the energy quanta at the resonant wavelength and σ GSA is the ground-state absorption (GSA) of resonantly-absorbing dope centers embedded in the material), α 0 = σ GSA N 0 is the weak-signal resonant absorption coefficient (where N 0 is the concentration of dope centers), γ is the residual (non-saturable) loss coefficient, and z is the coordinate inside the sample. Eq. (1) holds if the coefficient γ relates (i) to all the "passive" losses (owing to scattering or presence of non-resonant impurities) and (ii) to the loss stemming from the excited-state absorption (ESA) [2] in the system of dope centers (if such a source of loss exists, see Fig. 1). Notice that "slow" relaxation of the excited state means that decay time τ rel of the latter is essentially longer than pulse-width τ p .Solving Eq. (1), one is able to calculate nonlinear (energy-dependent) transmission coefficient of the medium as the ratio of normalized input (ε in ) and output (ε out ) pulse energy densities: T (ε in ) = ε out /ε in . Meanwhile, to the best of our knowledge there is not known an analogous equation concerning the resonant nonlinear change ∆n res (ε in ) in refractive index n of the medium that should accompany the transmission coefficient nonlinearity. Despite there have been published a few re-