A method is introduced to calc., for a model set of mol. levels, the spectral line shapes expected for a variety of conventional laser expts. including absorption, hole burning, fluorescence line narrowing, and Raman scattering. The method allows the incident laser field to have arbitrary intensity. Furthermore, the effects of model gaussian or lorenzian inhomogeneous distributions are readily incorporated. Earlier results for a 2-level system are easily obtained and new results are presented for inhomogeneously broadened 2- and 3-level systems, and for the effects of pure dephasing on the strong field spectra. The differences between fluorescence and Raman in strong fields, and the effect of strong fields on the spontaneous emission of inhomogeneously broadened transitions were described. Some predictions are made regarding line narrowing expts. in the strong-field limit.