“…where E n ( r) is electric field envelope at output mirror, f is nonlinear gain transfer function, L r is cavity length, L a << L r is a thickness of amplifying medium relevant to experimental situation [12], R is reflectivity of output mirror, D( r) = R exp[−| r| 2 /D 2 0 ] is diaphragm function corresponding to transversely inhomogeneous linear losses, key parameter for comparison with experiments [3] was D 0 ∼ 50 − 2000µm , G n ( r) = σN n ( r)L a is transversely inhomogeneous linear gain, σ is simulated emission cross-section, N n ( r) is density of resonant ions per unit volume, δω is detuning from gain line center, T 2 and T 1 are transverse and longitudinal relaxation times correspondingly, δE( r) is random noise due to spontaneous emission term emulated as multimode random process [13], δN ( r) is fluctuating part of optical pump also defined as spatially smooth multimode random process, N f = D 2 0 /(λL r ) is Fresnel number, λ is lasing wavelength. The convergence rate to equilibrium solutions (stationary eigenmodes) for dicrete time step δt = 2L r n/c proved to be 20 -100 iterates [14].…”