We have shown that the polarization component method can be extended to the case of broadband emission (multiple-frequency approximation) if we make the approximation of slowly varying amplitudes and statistical independence of the polarization components. From the results of numerical modeling for the case of a liquid dye laser, it follows that the multiple-frequency approximation is meaningful to use for analysis of the spectral, polarization, and energy parameters in the region near and below the lasing threshold, and also for the statistical characteristics but taking into account saturation effects for the laser emission.Introduction. In a series of papers [1-5], a new generalized approach was proposed and developed to describing the formation of polarized radiation in laser systems. It was based on the idea of laser emission as the superposition of polarization components whose polarization vectors span the entire set of possible polarized states. Based on the polarization component method, it was possible to explain a number of characteristic features on the energy and polarization characteristics of a liquid dye laser near the lasing threshold [1,2,5], and also to derive analytical relationships to describe their statistical properties [4]. In particular, distribution functions of the degree of polarization of the output radiation were found [4], and their variation as a function of the laser system parameters was studied. Based on the polarization component method, it was possible to quite correctly describe the appearance of the complex phase polarization structure of liquid dye laser radiation observed in the experiment [6,7]. We should note that, although initially the polarization component method was developed with an aim toward liquid dye lasers, it has been demonstrated that it can be used to analyze dynamic states in the polarization instability region in VCSEL systems [8].However, the developed polarization component method involves two aspects which considerably limit its applicability as a general approach to describing formation of polarized radiation in laser systems. The first aspect is connected with the single-frequency approximation, the basis for all subsequent conclusions. Generally speaking, such an approximation is widely used in laser theory, including to describe polarization effects [9]. At the same time, within the polarization component method, the process of formation of laser emission is described as successive enhancement of spontaneous emission, so taking into account the spectral composition of the latter seems quite reasonable. Furthermore, it has been shown [10-12] that, for example, in semiconductor injection lasers, we need to take into account enhanced luminescence in principle for correct determination of threshold and dynamic characteristics.A second limitation of the developed approach was neglect of laser intensity saturation effects, which considerably simplified the calculations and made it possible in principle to construct analytical expressions for the distrib...