Nonlinear traveling waves of the Klein-Gordon equation with cubic nonlinearity are considered. These waves are described by the nonlinear ordinary differential equation of the second order having the energy integral. Linearized equation for variation obtained for such waves is transformed to the ordinary one using separation of variables. Then so-called algebraization by Ince is used. Namely, a new independent variable associated with the solution under consideration is introduced to the equation in variations. Integral of energy for the stationary waves is used in this transformation. An advantage of this approach is that an analysis of the stability problem does no need to use the specific form of the solution under consideration. As a result of the algebraization, the equation in variations with variable in time coefficients is transformed to equation with singular points. Indices of the singularities are found. Necessary conditions of the waves stability are obtained. Solutions of the variational equation, corresponding to boundaries of the stability/instability regions in the system parameter space, are constructed in power series by the new independent variable. Infinite recurrent systems of linear homogeneous algebraic equations to determine coefficients of the series can be written. Non-trivial solutions of these systems can be obtained if their determinants are equal to zero. These determinants are calculated up to the fifth order inclusively, then relations connecting the system parameters and corresponding to boundaries of the stability/ instability regions in the system parameter place are obtained. Namely, the relation between parameters of anharmonicity and energy of the waves are constructed. Analytical results are illustrated by numerical simulation by using the Runge-Kutta procedure for some chosen parameters of the system. A correspondence of the numerical and analytical results is observed.