A method has been developed for studying the secular evolution and stabilization of the shape of rings in small celestial bodies that do not have shepherd satellites. A model of a compound ring consisting of two close, generally non-coplanar elliptical Gaussian rings has been constructed. The self-gravitation of the ring is taken into account through the mutual gravitational energy of the boundary rings. The function is presented as a series with an accuracy of up to the 4th power of small eccentricities and mutual inclination of the rings. The secular evolution of a compound ring is described by differential equations in special (collective) variables. For rings without a central body (problem 1), a closed system of 8 differential equations is obtained using the mutual energy function. The evolution of rings in the azimuthally averaged potential of a rotating triaxial body is also studied (problem 2); a second system of eight differential equations is derived for it. In both problems, in addition to the general case, two particular ones are considered: i) the case of coplanar elliptical rings, and ii) the case of circular rings with a tilt. The theory is applied to study the recently discovered ring of the dwarf planet Haumea. It is shown that without taking into account self-gravity, the nodal precession time of the Haumea ring is equal to but taking into account the self-gravity of the ring can reduce this period. It is established that self-gravity does indeed contribute to the preservation of the ring shape without invoking the hypothesis of shepherd satellites. Criteria for the preservation of the ring shape are obtained, which made it possible to estimate the interval for the ratio of the ring mass to the mass of Haumea. Taking into account the optical thickness of the ring, it is shown that the Haumea ring with a mass can consist of ice particles of size d0≈0.7÷1 m.