Propagation of an asymmetric Gaussian beam in a medium with saturated nonlinear refractive index is analyzed using a "collective variable approach" to solve the general nonlinear Schrödinger equation and compared with that of a symmetric beam in both lossless and lossy media. For a lossless medium, we construct a diagram which defines regions of oscillatory and diffractive propagation of an asymmetric beam and compare it with that of a symmetric beam. We detect breathing dynamics of the widths and amplitude of the asymmetric beam in the oscillatory regime of propagation and identify two different types of width and amplitude beating, Type 1 and 2, depending on the initial beam energy and saturation constant of the medium. This is in contrast to a cubic-quintic medium where only one type of beating is obtained. Spatial solitons have become the focus of many theoretical and experimental investigations [1-4] due to their possible applications in optical communications technology and alloptical switching. It is well known that stable self-guided beams can be created in two transverse dimensions if nonlinearities of higher orders (compared to the Kerr nonlinearity) are taken into account. This is true for both lossless cubic-quintic media [5] and saturable media [6]. It has also been shown that the propagation of Gaussian beams can be accompanied by periodic oscillations of the amplitude (breathing) and diffractive beam propagation in cubic-quintic media [5] and in saturable media [7,8].The theoretical description of wave propagation in both types of nonlinear media is based on an analysis of the general nonlinear Schrödinger equation (GNLSE). Due to the absence of exact analytical solutions of the GNLSE, approximate techniques are used to find spatially localized solutions that preserve their shape during propagation. Although numerical solutions are preferable for accuracy, they are very time consuming [9]. However, considerable physical insight can be gained if the GNLSE is solved approximately in a semianalytical fashion [1,2] using the "collective variable approach" (CVA) [10,11] for conservative systems. The CVA technique is based on a trial function, such as a Gaussian function, with a finite number of variables, which is usually a function of the propagation coordinate that evolves subject to the constraints of the system.The CVA technique was recently used to investigate propagation of an asymmetric Gaussian beam in a lossy cubicquintic medium [12], which is characterized by competition between a self-focusing χ (3) Kerr nonlinearity and a selfdefocusing χ (5) nonlinearity. It was shown that there are two main regions of behavior, periodic and non-periodic. The periodic region is further divided into two parts: an oscillatory self-focusing part and an oscillatory diffractive part. However, the sharp boundary between them that exists for a symmetric beam does not exist in the case of an asymmetric beam. For such a beam, an interesting phenomenon was found in the oscillatory region: the amplitude and widths of ...