We consider a symmetric five-body problem with three unequal collinear masses on the axis of symmetry. The remaining two masses are symmetrically placed on both sides of the axis of symmetry. Regions of possible central configurations are derived for the four- and five-body problems. These regions are determined analytically and explored numerically. The equations of motion are regularized using Levi-Civita type transformations and then the phase space is investigated for chaotic and periodic orbits by means of Poincar\'e surface of sections.

The equilibrium points and the curves of zero-velocity (Roche varieties) are analysed in the frame of the regularized circular restricted three-body problem. The coordinate transformation is done with Levi-Civita generalized method, using polynomial functions of n degree. In the parametric plane, five families of equilibrium points are identified. These families of points correspond to the five equilibrium points in the physical plane L1, L2, ..., L5. The zero-velocity curves from the physical plane are transformed in Roche varieties in the parametric plane. The properties of these varieties are analysed and the Roche varieties for n = {1,2,...,6} are plotted. The equation of the asymptotic variety is obtained and its shape is analysed. The slope of the Roche variety in L11 point is obtained. For n = 1 the slope obtained by Plavec and Kratochvil (1964) in the physical plane was found.

A family of polynomial coupled function of $n$ degree is proposed, in order to generalize the Levi-Civita regularization method, in the restricted three-body problem. Analytical relationship between polar radii in the physical plane and in the regularized plane are established; similar for polar angles. As a numerical application, trajectories of the test particle using polynomial functions of $2, 3,..., 8$ degree are obtained. For the polynomial of second degree, the Levi-Civita regularization method is found.

The regularization of a new problem, namely the three-body problem, using 'similar' coordinate system is proposed. For this purpose we use the relation of 'similarity', which has been introduced as an equivalence relation in a previous paper (see \cite{rom11}). First we write the Hamiltonian function, the equations of motion in canonical form, and then using a generating function, we obtain the transformed equations of motion. After the coordinates transformations, we introduce the fictitious time, to regularize the equations of motion. Explicit formulas are given for the regularization in the coordinate systems centered in the more massive and the less massive star of the binary system. The 'similar' polar angle's definition is introduced, in order to analyze the regularization's geometrical transformation. The effect of Levi-Civita's transformation is described in a geometrical manner. Using the resulted regularized equations, we analyze and compare these canonical equations numerically, for the Earth-Moon binary system.