Parametric Expansions of an Algebraic Variety Near Its Singularities II

Abstract

The paper is a continuation and completion of the paper Bruno, A.D.; Azimov, A.A. Parametric Expansions of an Algebraic Variety Near Its Singularities. Axioms 2023, 5, 469, where we calculated parametric expansions of the three-dimensional algebraic manifold $\Omega$, which appeared in theoretical physics, near its 3 singular points and near its one line of singular points. For that we used algorithms of Nonlinear Analysis: extraction of truncated polynomials, using the Newton polyhedron, their power transformations and Formal Generalized Implicit Function Theorem. Here we calculate parametric expansions of the manifold $\Omega$ near its one more singular point, near two curves of singular points and near infinity. Here we use 3 new things: (1) computation in algebraic extension of the field of rational numbers, (2) expansions near a curve of singular points and (3) calculation of branches near infinity.

1. Introduction

Here we continue and conclude the paper [1]. There, in Sections 1–5, we proposed a new method for solving the polynomial equation

$$f\left(x_1,\dots ,x_n\right)=0$$

near a singular point or curve of singular points of the polynomial f. In Sections 6–10, this method was applied to compute the solutions to such a 12th degree equation with $n=3$ that originated in theoretical physics. This new method is based on:

• Newton’s polyhedron to isolate the truncated equations,
• Power transformations to simplify those equations, and
• Formal Generalized Implicit Function Theorem to obtain solutions in the form of power expansions whose coefficients are rational functions of the parameters. Computer algebra is used in these calculations.

Newton’s polyhedron is a multidimensional generalization of Newton’s polygon (see [2,3,4,5,6,7]). Power transformations are a generalization of the sigma process used previously to resolve singularities of algebraic manifolds (see [8,9,10]). Algorithms for computing power transformations were proposed in [11]. The resolution of the singularity is done step-by-step until we come to the situation with a truncated equation containing a polynomial multiplier of degree one. If the roots of this multiplier are parameterized, the roots of the whole polynomial are obtained as a power expansion using a Generalized Implicit Function Theorem (Theorem 1 of [1]). All these are an application to algebraic equations of the general theory of Nonlinear Analysis [12], which is also suitable for differential equations. For its applications to systems of partial derivative equations, see [13].

According to [1] and [12] (Section 2) computational steps are the following:

Step 1. 

Introduction of local coordinates. For coordinates $X=(x_1,\dots ,x_n)$ and singular point $X^0=(x_1^0,\dots ,x_n^0)$, they are $Y=X-X^0$, i.e., $y_i=x_i-x_i^0$, $i=1,\dots ,n$.

Step 2. 

Writing the initial polynomial $f\left(X\right)$ in local coordinates

$$g\left(Y\right)=f(X^0+Y)=\sum g_Q{Y}^Q\mathrm{over}Q\in \mathbf{S}.$$

(1)

Here $Q=(q_1,\dots ,q_n)$, $Y^Q=y_1^{q_1}\dots y_n^{q_n}$, $g_Q$ are real or complex coefficients, the sum has not similar terms, the set $\mathbf{S}\left(g\right)=\{Q:a_Q\not\equiv 0\}$, is called as support of the sum $g\left(Y\right)$. Here $0⩽Q\in {\mathbb{Z}}^n$. Let the support $\mathbf{S}\left(g\right)$ consists of vectors $Q_1,\dots ,Q_k$.

Step 3. 

The Newton polyhedron $\Gamma \left(g\right)$ is computating as the convex hull of the support $\mathbf{S}$:

$$\Gamma \left(g\right)=\left\{Q=\sum _{i=1}^k{\lambda }_i{Q}_i,{\lambda }_i⩾0,i=1,\dots ,k,\sum _{i=1}^k{\lambda }_i=1\right\}.$$

The boundary $\partial \Gamma$ of the polyhedron $\Gamma \left(g\right)$ consists from its generalized faces ${\Gamma }_j^{\left(d\right)}$, where d is dimension, $0⩽d⩽n-1$, and j is a number of the face ${\Gamma }_j^{\left(d\right)}$. Each face ${\Gamma }_j^{\left(d\right)}$ corresponds to its truncated polynomial

$${\widehat{g}}_j^{\left(d\right)}\left(Y\right)=\sum g_Q{Y}^Q\mathrm{over}Q\in \mathbf{S}\cap {\Gamma }_j^{\left(d\right)}$$

and the normal cone ${\mathbf{U}}_j^{\left(d\right)}$, consisting of all normals to the face ${\Gamma }_j^{\left(d\right)}$, which are external to the polyhedron $\Gamma$. For their computation we use the PolyhedralSets package of the computer algebra system (CAS) Maple. In the steps below $n=3$. Then ${\Gamma }_j^{\left(2\right)}$ is two dimensional face and normal cone ${\mathbf{U}}_j^{\left(2\right)}$ is a ray, spanned by external normal $N_j$ to the face ${\Gamma }_j^{\left(2\right)}$.

Step 4. 

We select faces ${\Gamma }_j^{\left(2\right)}$ with normal $N_j⩽0$ and corresponding truncated polynomials ${\widehat{g}}_j^{\left(2\right)}\left(Y\right)$.

Step 5. 

For each selected truncated polynomial ${\widehat{g}}_j^{\left(2\right)}\left(Y\right)$, we compute corresponding power transformation

$$\left(ln{y}_1,ln{y}_2,ln{y}_3\right)=\left(ln{z}_1,ln{z}_2,ln{z}_3\right)\alpha ,$$

(2)

where $\alpha$ is an unimodular $3\times 3$ matrix, such that

$${\widehat{g}}_j^{\left(2\right)}\left(Y\right)=h(z_1,z_2)z_3^l$$

(3)

with integral l.

Step 6. 

If the curve $h(z_1,z_2)=0$ has parametrization

$$z_1=b_1\left(t\right),z_2=b_2\left(t\right),$$

then it is obtained with the algcurves package from the CAS Maple. In that case we make the power transformation (2) in the full polynomial (1) and write it as

$$g\left(Y\right)=T(z_1,z_2,z_3)z_3^l=z_3^l\sum _{k=0}^m{T}_k(z_1,z_2)z_3^k,$$

with some natural m. Here polynomials $T_k(z_1,z_2)$ are computed by the command coeff(T,z[k],m) in CAS Maple. Here $T_0(z_1,z_2)=h(z_1,z_2)$ from (3).

Step 7. 

If $T_1(b_1\left(t\right),b_2\left(t\right))\not\equiv 0$, we make the substitution

$$z_1=b_1\left(t\right)+\epsilon ,z_2=b_2\left(t\right)+\epsilon$$

(4)

into the polynomial $T(z_1,z_2,z_3)$, obtain function $u(\epsilon ,t,z_3)=T(z_1,z_2,z_3)$, apply to the equation $u(\epsilon ,t,z_3)=0$ the Formal Generalized Implicit Function Theorem 1 [1] and get the parametric expansion

$$\epsilon =\sum _{k=1}^{\infty }{c}_k\left(t\right)z_3^k.$$

(5)

Step 8. 

We compute several terms of expansion (5), substitute them into (4). The result is substituted in power transformation (2), and we obtain parametric expansion of Y in power series of $z_3$ with coefficients which are rational functions of t.

If $T_1(b_1\left(t\right),b_2\left(t\right))\equiv 0$, we continue computation with new Newton polyhedron etc.

The method is new, with parts: the Newton polyhedron $\Gamma \left(g\right)$, polyhedron’s faces ${\Gamma }_j^{\left(d\right)}$, polyhedron graph, normal cones ${\mathbf{U}}_j^{\left(d\right)}$ and power transformations (2) were proposed by the first author beginning 1962. Early such objects he calculated manually, but now there are programs for that.

In [1], this theory was applied to a problem arises in the study of Ricci flows (see [14,15,16,17,18,19,20,21,22]). The Ricci flows describe the evolution of Einstein’s metrics on a variety. The equations of the normalized Ricci flow are reduced to a system of two differential equations with three parameters: $a_1$, $a_2$ and $a_3$:

$$\begin{array}{cc}\hfill \frac{d{x}_1}{dt}& ={\tilde{f}}_1(x_1,x_2,a_1,a_2,a_3),\hfill \\ \hfill \frac{d{x}_2}{dt}& ={\tilde{f}}_2(x_1,x_2,a_1,a_2,a_3),\hfill \end{array}$$

(6)

here, ${\tilde{f}}_1$ and ${\tilde{f}}_2$ are certain given functions.

The singular points of this system are associated with the invariant Einstein’s metrics. At the singular (stationary) point $x_1^0$, $x_2^0$, system (6) has two eigenvalues, ${\lambda }_1$ and ${\lambda }_2$. If at least one of them is equal to zero, then the singular (stationary) point $x_1^0$, $x_2^0$ is said to be degenerate. It was proved in [14,15,16,17,18,19,20,21,22] that the set $\Omega$ of the values of the parameters $a_1$, $a_2$, $a_3$, in which system (6) has at least one degenerate singular point, is described by all solutions of the equation

$$\begin{array}{cc}\hfill Q(s_1,s_2,s_3)\stackrel{\mathrm{def}}{\equiv }& (2s_1+4s_3-1)\left(64s_1^5-64s_1^4+8s_1^3+240s_1^2{s}_3-1536s_1{s}_3^2-\right.\hfill \\ \hfill & \left.-4096s_3^3+12s_1^2-240s_1{s}_3+768s_3^2-6s_1+60s_3+1\right)-\hfill \\ \hfill & -8s_1{s}_2(2s_1+4s_3-1)(2s_1-32s_3-1)(10s_1+32s_3-5)-\hfill \\ \hfill & -16s_1^2{s}_2^2\left(52s_1^2+640s_1{s}_3+1024s_3^2-52s_1-320s_3+13\right)+\hfill \\ \hfill & +64(2s_1-1)s_2^3(2s_1-32s_3-1)+2048s_1(2s_1-1)s_2^4=0,\hfill \end{array}$$

where $s_1$, $s_2$, $s_3$ are elementary symmetric polynomials, equal, respectively, to

$$s_1=a_1+a_2+a_3,s_2=a_1{a}_2+a_1{a}_3+a_2{a}_3,s_3=a_1{a}_2{a}_3.$$

Here Q is different from Q’s in Steps 2 and 3, but the sign $\stackrel{\mathrm{def}}{\equiv }$ means only a new notation.

Hence the polynomial $P(a_1,a_2,a_3)=Q(s_1,s_2,s_3)$ has degree 12. In [23], for symmetry reasons, the coordinates $\mathbf{a}=(a_1,a_2,a_3)$ were changed to the coordinates $\mathbf{A}=(A_1,A_2,A_3)$ by the linear transformation

$$\left(\begin{array}{c}{a}_1\\ a_2\\ a_3\end{array}\right)=M·\left(\begin{array}{c}{A}_1\\ A_2\\ A_3\end{array}\right),M=\left(\begin{array}{ccc}(1+\sqrt{3})/6& (1-\sqrt{3})/6& 1/3\\ (1-\sqrt{3})/6& (1+\sqrt{3})/6& 1/3\\ -1/3& -1/3& 1/3\end{array}\right)$$

The resulting polynomial is

$$\begin{array}{c}\hfill R\left(\mathbf{A}\right)=P\left(\mathbf{a}\right)\end{array}$$

(7)

and has degree 12 again.

Definition 1.

Let $\phi \left(X\right)$ be some polynomial, where $X=(x_1,\dots ,x_n)$. A point $X=X^0$ of the set $\phi \left(X\right)=0$ is called the singular point of the k-order, if all partial derivatives of the polynomial $\phi \left(X\right)$ with respect to $x_1,\dots ,x_n$ turn into zero at this point, up to and including the k-th order derivatives, and at least one partial derivative of order $k+1$ is nonzero.

In [23], all singular points of the variety $\Omega$ in coordinates $\mathbf{A}=\left(A_1,A_2,A_3\right)$ were found. The five points of the third order are:

Name	Coordinates  $\mathbf{A}$
$P_1^{\left(3\right)}$	$(0,0,3/4)$
$P_2^{\left(3\right)}$	$(0,0,-3/2)$
$P_3^{\left(3\right)}$	$\left(-\frac{1+\sqrt{3}}{2},\frac{\sqrt{3}-1}{2},\frac{1}{2}\right)$
$P_4^{\left(3\right)}$	$\left(\frac{\sqrt{3}-1}{2},-\frac{1+\sqrt{3}}{2},\frac{1}{2}\right)$
$P_5^{\left(3\right)}$	$(1,1,1/2)$

three points of the second order

Name	Coordinates  $\mathbf{A}$
$P_1^{\left(2\right)}$	$\left(\frac{1+\sqrt{3}}{4},\frac{1-\sqrt{3}}{4},\frac{1}{2}\right)$
$P_2^{\left(2\right)}$	$\left(\frac{1-\sqrt{3}}{4},\frac{1+\sqrt{3}}{4},\frac{1}{2}\right)$
$P_3^{\left(2\right)}$	$(-1/2,-1/2,1/2)$

and three more algebraic curves of singular points of the first order:

$$\begin{array}{c}\mathcal{F}=\left\{a_1=a_2,16a_1^3+16a_1^2{a}_3-4a_1-2a_3+1=0\right\},\\ \mathcal{I}=\left\{A_1+A_2+1=0,A_3=\frac{1}{2}\right\},\\ \mathcal{K}=\left\{A_1=-\frac{9}{4}th\left(t\right),A_2=-\frac{9}{4}h\left(t\right),A_3=\frac{3}{4},h\left(t\right)=\frac{t^2+1}{(t+1)(t^2-4t+1)}\right\}.\end{array}$$

The points $P_3^{\left(3\right)}$, $P_4^{\left(3\right)}$ and $P_5^{\left(3\right)}$ are of the same type; they pass into each other when rotated in the plane $A_1,A_2$ by an angle $2\pi /3$, just as all points $P_1^{\left(2\right)}$, $P_2^{\left(2\right)}$, $P_3^{\left(2\right)}$. The curves $\mathcal{F}$, $\mathcal{I}$, $\mathcal{K}$ correspond to two more curves of the same type. Therefore, it is sufficient to study the variety $\Omega$ in the neighborhood of points $P_1^{\left(3\right)}$, $P_2^{\left(3\right)}$, $P_5^{\left(3\right)}$, $P_3^{\left(2\right)}$ and curves $\mathcal{F}$, $\mathcal{I}$ and $\mathcal{K}$. Moreover, in [23] there were computed sections of the variety $\Omega$ by planes $A_3=\mathrm{const}$, and was shown that in finite part of the space ${\mathbb{R}}^3=\left\{A_1,A_2,A_3\right\}$ the variety $\Omega$ consists of two dimensional branches $F_1$, $F_2$, $F_3$, $G_1$, $G_2$, $G_3$ divided into parts $F_i^{\pm }$, $G_i^{\pm }$ with boundaries at the plane $A_3=1/2$.

In the paper [24], three variants of the global parametrization of the variety $\Omega$ were proposed. These parametrizations were computed using the parametric description of the discriminant set of a monic cubic polynomial [25] and can be written in radical form [26]. Such a global description of the variety $\Omega$ cannot provide an adequate picture of the $\Omega$ structure in the vicinity of its singular points.

In [1], parametric expansions of the variety $\Omega$ near the singular points $P_1^{\left(3\right)}$ (Section 7), $P_1^{\left(3\right)}$ (Section 8), $P_2^{\left(2\right)}$ (Section 8), $P_3^{\left(2\right)}$ (Section 10) and near the line of singular points $\mathfrak{I}$ (Section 9) were computed. Here these expansions are computed near the singular point $P_5^{\left(3\right)}$ (Section 2), near the curves of singular points $\mathcal{H}$ (Section 3) and $\mathcal{F}$ (Section 4), and near infinity (Section 5). Together they cover a wide range of cases. The following tactic has developed: if the truncated equation contains linear multipliers, they are used to do a linear transformation of the coordinates followed by the computation of Newton’s polyhedron; and if they are nonlinear, a power transformation of the coordinates is done. To understand the present article it is necessary a knowledge with papers [1] (open access), [23] and the book [27].

2. The Structure of the Manifold $\Omega$ near the Singular Point ${\mathit{P}}_{\mathbf{5}}^{\left(\mathbf{3}\right)}$

2.1. Preliminary Computations

Near the point $P_5^{\left(3\right)}:\left(A_1,A_2,A_3\right)=(1,1,1/2)$ we introduce local coordinates $x_1,x_2,x_3$:

$$A_1=x_1+1,A_2=x_2+1,A_3=\frac{1}{2}+x_3.$$

(8)

Then, from the polynomial $R\left(\mathbf{A}\right)$ in (7) we get a polynomial of degree 12.

$$S_4\left(x_1,x_2,x_3\right)=R\left(\mathbf{A}\right)=Q\left(s_1,s_2,s_3\right).$$

We compute the support $\mathbf{S}$ of the polynomial $S_4$, the Newton polyhedron ${\Gamma }_4\left(S_4\right)$, its faces ${\Gamma }_j^{\left(2\right)}$ and their external normals using the PolyhedralSets package of the computer algebra system (CAS) Maple 2021 [27]. We obtain 5 faces ${\Gamma }_j^{\left(2\right)}$. The graph of the polyhedron ${\Gamma }_4$ is shown in Figure 1.

In the top row—the whole polyhedron, in the next—all two-dimensional faces. Further down are the edges, then the vertices, and at the bottom is the empty set. The external normals to its two-dimensional faces are

$$N_{71}=\left(-1,-1,-1\right),N_{143}=\left(1,1,1\right),N_{215}=\left(-1,0,0\right),N_{239}=\left(0,-1,0\right),N_{241}=\left(0,0,-1\right).$$

Since $x_1,x_2,x_3⟶0,$ we take the only normal that has all coordinates negative. This is $N_{71}=\left(-1,-1,-1\right)$. It corresponds to a truncated polynomial

$${\widehat{f}}_{71}=-16{\left(2x_1^2-8x_1{x}_2+2x_1{x}_3+2x_2^2+2x_2{x}_3-x_3^2\right)}^2/81.$$

(9)

The quadratic polynomial bracketed in (9) does not factorize in the field of rational numbers, but it does factorize in the extension of this field with $\sqrt{3}$. We get

$$f=2x_1^2-8x_1{x}_2+2x_1{x}_3+2x_2^2+2x_2{x}_3-x_3^2.$$

We proceed according to [27]:

• >alpha := RootOf(y^2-3);
• >factor (f,alpha);

$$\begin{array}{c}-\frac{1}{2}(\left(2x_{2}RootOf\left(Z^2-3\right)-x_{3}RootOf\left(Z^2-3\right)+2x_1-4x_2+x_3\right)·\hfill \\ ·\left(2x_{2}RootOf\left(Z^2-3\right)-x_{3}RootOf\left(Z^2-3\right)-2x_1+4x_2-x_3\right))\hfill \end{array}$$

i.e.,

$$f=-\frac{1}{2}\left(2x_1-\left(4-2\sqrt{3}\right)x_2+\left(1-\sqrt{3}\right)x_3\right)·\left(2x_1-\left(4+2\sqrt{3}\right)x_2+\left(1+\sqrt{3}\right)x_3\right).$$

Now we do a linear substitution of the coordinates

$$\begin{array}{cc}\hfill y_1& =2x_1-\left(4-2\sqrt{3}\right)x_2+\left(1-\sqrt{3}\right)x_3,\hfill \\ \hfill y_2& =2x_1-\left(4+2\sqrt{3}\right)x_2+\left(1+\sqrt{3}\right)x_3,\hfill \\ \hfill y_3& =x_3.\hfill \end{array}$$

Its inverse substitution is

$$\begin{array}{cc}\hfill x_1& =\frac{\left(2+\sqrt{3}\right)\sqrt{3}}{12}{y}_1+\frac{\left(-2+\sqrt{3}\right)\sqrt{3}}{12}{y}_2+\frac{1}{2}{y}_3,\hfill \\ \hfill x_2& =\frac{\sqrt{3}}{12}{y}_1-\frac{\sqrt{3}}{12}{y}_2+\frac{1}{2}{y}_3,\hfill \\ \hfill x_3& =y_3.\hfill \end{array}$$

(10)

We substitute it into the polynomial $S_4\left(\mathbf{x}\right)$ and get the polynomial $S_5\left(\mathbf{y}\right)=S_4\left(\mathbf{x}\right)$. For the polynomial $S_5\left(\mathbf{y}\right)$, we compute Newton’s polyhedron ${\Gamma }_5$. Its graph is shown in Figure 2. It has 11 two-dimensional faces with external normals

$$\begin{array}{c}{N}_{14397}=\left(-1,-1,0\right),N_{15959}=\left(-1,0,-1\right),N_{19269}=\left(-2,-2,-1\right),\\ N_{39917}=\left(0,-1,-1\right),N_{111761}=\left(1,1,1\right),N_{131145}=\left(-1,0,0\right),N_{132735}=\left(0,0,1\right),\\ N_{135677}=\left(0,1,0\right),N_{137855}=\left(1,0,0\right),N_{159459}=\left(0,-1,0\right),N_{162019}=\left(0,0,-1\right).\end{array}$$

We parse the first 4 of them that have 2 or 3 coordinates negative, dedicating a subsection to each of them. Below we use notation from Maple [27].

2.2. The Normal $(-2,-2,-1)$

The corresponding truncated polynomial is

$$\begin{array}{cc}ftr19269=\hfill & -104976(5+3\alpha )(864\alpha y_1{y}_3^6-3456\alpha y_2{y}_3^6+117\alpha y_1^2{y}_3^4-594\alpha y_1{y}_2{y}_3^4\hfill \\ & \hspace{1em}\hspace{1em}+1287\alpha y_2^2{y}_3^4-1728y_1{y}_3^6+6048y_2{y}_3^6-12\alpha y_1^2{y}_2{y}_3^2+48\alpha y_1{y}_2^2{y}_3^2-360\alpha y_2^3{y}_3^2-117y_1^2{y}_3^4\hfill \\ & \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}+990y_1{y}_2{y}_3^4-2223y_2^2{y}_3^4+6\alpha y_1^2{y}_2^2-24y_1^3{y}_3^2+24y_1^2{y}_2{y}_3^2-84y_1{y}_2^2{y}_3^2+624y_2^3{y}_3^2-10y_1^2{y}_2^2),\hfill \end{array}$$

where $\alpha =\sqrt{3}$. Here and below “$ftr$” number means truncated polynomial, corresponding to normal $N_j$ with written number j. According to [11], we compute the matrix $\gamma =\left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ -2& -2& 1\end{array}\right)$ such that $\left(-2,-2,-1\right)\gamma =\left(0,0,-1\right)$. Since ${\gamma }^{-1}=\left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 2& 2& 1\end{array}\right)$, then we do a power transformation

$$y_1=z_1{z}_3^2,y_2=z_2{z}_3^2,y_3=z_3.$$

(11)

We get after factorization

$$ftr19269=-104976(5+3\alpha )(6\alpha z_1^2{z}_2^2-12\alpha z_1^2{z}_2+48\alpha z_1{z}_2^2-360\alpha z_2^3-10z_1^2{z}_1^2{z}_2^2$$

$$+117\alpha z_1^2-594\alpha z_1{z}_2+1287z_2^2\alpha -24z_1^3+24z_1^2{z}_2-84z_1{z}_2^2+624z_2^3+864\alpha z_1$$

$$-3456\alpha z_2-117z_1^2+990z_1{z}_2-2223z_2^2-1728z_1+6048z_2)z_3^8,$$

where $\alpha =\sqrt{3}$. The large polynomial in the parentheses is denoted by

$$f_1(z_1,z_2)=6\alpha z_1^2{z}_2^2-12\alpha z_1^2{z}_2+48\alpha z_1{z}_2^2-360\alpha z_2^3-10z_1^2{z}_1^2{z}_2^2$$

$$+117\alpha z_1^2-594\alpha z_1{z}_2+1287z_2^2\alpha -24z_1^3+24z_1^2{z}_2-84z_1{z}_2^2+624z_2^3+864\alpha z_1$$

$$-3456\alpha z_2-117z_1^2+990z_1{z}_2-2223z_2^2-1728z_1+6048z_2.$$

Consider the curve $f_1\left(z_1,z_2\right)=0.$ It has intersections with the axes: $z_1=0,z_2=0.$ It is a curve of genus 0, with one singular point

$$\left(z_1,z_2\right)=\left(-6+6\sqrt{3};-6-6\sqrt{3}\right)\approx \left(4.39230485;-16.39230485\right)$$

and parameterization

$$\begin{array}{cc}\hfill z_1& =b_1\left(t\right)=\frac{\left(5+3\sqrt{3}\right)\beta \left(2042820\sqrt{3}t-521639194050t^2+64\sqrt{3}-7660575t-119\right)}{2{\left(510705t+2\right)}^2},\hfill \\ \hfill z_2& =b_2\left(t\right)=\frac{\left(153\sqrt{3}+265\right)\beta \left(2042820\sqrt{3}t-260819597025t^2+120\sqrt{3}-5617755t-212\right)}{2\left(510705t+2\right)},\hfill \\ \hfill \beta & =3\left(4\sqrt{3}+510705t-5\right).\hfill \end{array}$$

We simplify these expressions by transforming $t=t_1/510705$ and obtaining

$$\begin{array}{cc}\hfill z_1& =b_1\left(t_1\right)=\frac{3\left(5+3\sqrt{3}\right)\left(4\sqrt{3}{t}_1-2t_1^2+64\sqrt{3}-15t_1-119\right)\left(4\sqrt{3}+t_1-5\right)}{2{\left(t+2\right)}^2},\hfill \\ \hfill z_2& =b_2\left(t_1\right)=\frac{3\left(265+153\sqrt{3}\right)\left(4\sqrt{3}{t}_1-t_1^2+120\sqrt{3}-11t_1-212\right)\left(4\sqrt{3}+t_1-5\right)}{2\left(t_1+2\right)}.\hfill \end{array}$$

(12)

The curve is shown in Figure 3.

With $z_3$ fixed according to (8), (10) and (11) at the original coordinates $\mathbf{A}$, we obtain a curve

$$\begin{array}{cc}\hfill A_1& =1+\frac{(2+\sqrt{3})\sqrt{3}}{12}{z}_1{z}_3^2+\frac{(-2+\sqrt{3})\sqrt{3}}{12}{z}_2{z}_3^2+\frac{1}{2}{z}_3,\hfill \\ \hfill A_2& =1+\frac{\sqrt{3}}{12}{z}_1{z}_3^2-\frac{\sqrt{3}}{12}{z}_2{z}_3^2+\frac{1}{2}{z}_3,\hfill \end{array}$$

(13)

where $z_1=b_1\left(t_1\right)$, $z_2=b_2\left(t_1\right)$ according to (12).

If $z_3=-0.05$, $z_1=b_1\left(t_1\right)$, $z_2=b_2\left(t_1\right)$ then according to (8), (10) and (13) $A_3=9/20=0.45$,

$$\begin{array}{cc}\hfill A_1& =\frac{(2+\sqrt{3})\sqrt{3}}{4800}{z}_1+\frac{(-2+\sqrt{3})\sqrt{3}}{4800}{z}_2+\frac{39}{40},\hfill \\ \hfill A_2& =\frac{\sqrt{3}}{4800}{z}_1-\frac{\sqrt{3}}{4800}{z}_2+\frac{39}{40}.\hfill \end{array}$$

(14)

This curve is shown in Figure 4. This curve is similar to the curve of [23] (Figure 12) showing the cross-section of the variety $\Omega$ at $A_3=0.45$, with branches $F_1^{-}$ and $G_2^{-}$.

If $z_3=0.005$, then $A_3=101/200=0.505$, hence according to (13)

$$\begin{array}{cc}\hfill A_1& =\frac{(2+\sqrt{3})\sqrt{3}}{480000}{z}_1+\frac{(-2+\sqrt{3})\sqrt{3}}{480000}{z}_2+\frac{401}{400},\hfill \\ \hfill A_2& =\frac{\sqrt{3}}{480000}{z}_1-\frac{\sqrt{3}}{480000}{z}_2+\frac{401}{400}.\hfill \end{array}$$

(15)

When $z_1=b_1\left(t_1\right)$, $z_2=b_2\left(t_1\right)$, the curve (15) is shown in Figure 5. It is similar to Figure 15 in [23], which shows the section of the variety $\Omega$ at $A_3=0.505$, with coinsiding branches $F_1^{+}$ and $G_2^{+}$.

In fact, the parametric expansion of the variety $\Omega$ can also be computed here. To do this, we substitute (11) into the polynomial $S_5\left(\mathbf{y}\right)$ and get the polynomial of Step 6

$$T\left(\mathbf{z}\right)=z_3^8\sum _{k=0}^m{T}_k(z_1,z_2)z_3^k,$$

(16)

where polynomials $T_k(z_1,z_2)$ are uniquelly determined and can be computed by the command coeff(T,z[k],m).

In it, according to (12), we make the substitution

$$z_1=b_1\left(t\right)+\epsilon ,z_2=b_2\left(t\right)+\epsilon .$$

(17)

We obtain that $T\left(\mathbf{z}\right)/z_3^8=u(\epsilon ,z_3)$ with coefficients depending on t via $b_1\left(t\right)$ and $b_2\left(t\right)$. We apply Theorem 1 in [1] to the equation $u\left(\epsilon ,z_3\right)=0$ and obtain the expansion

$$\epsilon =\sum _{k=1}^{\infty }{c}_k\left(t\right)z_3^k.$$

Returning to initial coordinates $\mathbf{A}$ via (8), (10)–(12) and (17), we obtain expansions

$$\begin{array}{cc}\hfill A_1& =1+\frac{1}{2}{z}_3+\frac{2\sqrt{3}+3}{12}\left(b_1\left(t_1\right)+\sum _{k=1}^{\infty }{c}_k\left(t_1\right)z_3^k\right)z_3^2+\frac{3-2\sqrt{3}}{12}\left(b_2\left(t_1\right)+\sum _{k=1}^{\infty }{c}_k\left(t_1\right)z_3^k\right)z_3^2,\hfill \\ \hfill A_2& =1+\frac{1}{2}{z}_3+\frac{\sqrt{3}}{12}\left(b_1\left(t_1\right)+\sum _{k=1}^{\infty }{c}_k\left(t_1\right)z_3^k\right)z_3^2-\frac{\sqrt{3}}{12}\left(b_2\left(t_1\right)+\sum _{k=1}^{\infty }{c}_k\left(t_1\right)z_3^k\right)z_3^2,\hfill \\ \hfill A_3& =\frac{1}{2}+z_3\hfill \end{array}$$

(18)

with parameters $t_1\in \mathbb{R}$ and $z_3\in \mathbb{R}$ for small $|z_3|$. Formula (13) contain first terms of expansions (18).

2.3. The Normal $(-1,-1,0)$

It corresponds to a truncated polynomial

$$ftr14397=-11337408\left(1+\alpha \right)\left(\alpha y_2+y_1-2y_2\right)\left(4y_3-1\right){\left(y_3+2\right)}^3{y}_3^6.$$

(19)

Since it is linear on $y_1,y_2$, its root is the $y_3$-axis, i.e., $y_1=y_2=0$, which we denote by N. This line N lies on the variety $\Omega$ and through N it passes one of its branches, which in the first approximation has the form $y_1=(2-\alpha )y_2$.

According to (8) and (9) in the original coordinates $\mathbf{A}$, this line N has the following form

$$A_1=\frac{1}{2}{y}_3+1,A_2=\frac{1}{2}{y}_3+1,A_3=\frac{1}{2}+y_3.$$

This is the straight line $g_2$ of [23] (Figure 3). In [23] (Figures 4–15), the points of the line N lie in the plane $M:A_1=A_2$. Moreover, in Figure 6 $A_1\in \left(-1,0\right)$, in Figures 8–11 $A_1\in \left(0,1\right)$, in Figures 4, 13 and 14 $A_1\in \left(1,2\right)$. According to (19), there are 3 singular points on the line N

$$y_3=0,y_3=-2,y_3=1/4.$$

(20)

In $\mathbf{A}$ coordinates, they look like this:

$$\left(1,1,1/2\right),\left(0,0,-3/2\right),\left(9/8,9/8,3/4\right),$$

i.e., they are points $P_5^{\left(3\right)},P_2^{\left(3\right)}$, and a point with $t=1$ on the curve $\mathcal{H}$ of singular points. The structure of the variety $\Omega$ near the point $P_5^{\left(3\right)}$ is dealt with in this Section, near the point $P_2^{\left(3\right)}$ was dealt with in Section 7 of [1], and in the next Section we will deal with the structure of the variety $\Omega$ near the curve $\mathcal{H}$.

We can obtain an expansion of the variety $\Omega$ near the line N. To do this, we substitute (10) in the polynomial $S_4\left(\mathbf{x}\right)$ and obtain the polynomial $V\left(y_1,y_2,y_3\right)=S_4\left(\mathbf{x}\right),$ which we write as

$$V\left(y_1,y_2,y_3\right)=\sum V_{q_1{q}_2}\left(y_3\right)y_1^{q_1}{y}_2^{q_2},$$

where $0\le q_1,q_2\in \mathbb{Z}$, $V_{q_1{q}_2}\left(y_3\right)-$ polynomials. Thus according to (19)

$$V_{00}=0,$$

$$V_{10}=-11337408(1+\alpha )(4y_3-1){(y_3+2)}^3{y}_3^6,$$

$$V_{01}=-11337408\left(1+\alpha \right)\left(\alpha -2\right)\left(4y_3-1\right){\left(y_3+2\right)}^3{y}_3^6.$$

According to [1] (Theorem 1), the equation $V=0$ has a solution

$$y_1=\left(2-\alpha \right)y_2+\sum _{k=2}^{\infty }{c}_k\left(y_3\right)y_2^k.$$

Going to the original coordinates, we get the expansion for $\mathbf{A}$ with parameters $y_2$ and $y_3$. It is valid everywhere except the neighborhoods of the points (20).

2.4. The Normal $(-1,0,-1)$

It corresponds to a truncated polynomial

$$ftr15959=-18\left(-26+15\alpha \right){\left(-y_2+6+6\alpha \right)}^2{\left(-y_2+3+3\alpha \right)}^2·(-18y_2^2{y}_3^2\alpha +6\alpha y_1^2{y}_2$$

$$-108y_2{y}_3^2\alpha -y_1^2{y}_2^2+36y_2^2{y}_3^2-18\alpha y_1^2+6y_1^2{y}_2+108y_2{y}_3^2-36y_1^2)·y_2^2.$$

Let’s put

$$f_2=-18y_2^2{y}_3^2\alpha +6\alpha y_1^2{y}_2-108y_2{y}^2{y}_3^2\alpha -y_1^2{y}_2^2+36y_2^2{y}^2{y}_3^2-18\alpha y_1^2+6y_1^2{y}_2+108y_2{y}_3^2-36y_1^2.$$

According to [11], we compute the unimodular matrix $\gamma =\left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ -1& 0& 1\end{array}\right)$ such that $\left(-1,0,-1\right)\gamma =\left(0,0,-1\right)$. Since ${\gamma }^{-1}=\left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 1& 0& 1\end{array}\right)$, then we do a power transformation

$$y_1=z_1{z}_3,y_2=z_2,y_3=z_3.$$

(21)

We get

$$f_2\left(z_1,z_2,z_3\right)=\left(6\alpha z_1^2{z}_2-z_1^2{z}_2^2-18\alpha z_1^2-18\alpha z_2^2+6z_1^2{z}_2-108z_2\alpha -36z_1^2+36z_2^2+108z_2\right)z_3^2$$

Let’s denote

$$\begin{array}{c}\hfill g\left(z_1,z_2\right)=6\alpha z_1^2{z}_2-z_1^2{z}_2^2-18\alpha z_1^2-18\alpha z_1^2-18\alpha z_2^2+6z_1^2{z}_2-108z_2\alpha -36z_1^2\\ \hfill +36z_2^2+108z_2=-z_1^2{z}_2^2+6\left(\alpha +1\right)z_1^2{z}_2-18\left(\alpha +2\right)z_1^2+18\left(-\alpha +2\right)z_2^2+108(-\alpha +1)z_2.\end{array}$$

(22)

The curve $g\left(z_1,z_2\right)=0$ has genus 0, intersections with axes

$${(z}_1,z_2)=\left(0,0\right),{(z}_1,z_2)=\left(0,-\frac{6\left(\sqrt{3}-1\right)}{\sqrt{3}-2}\right)=(0,16.39230485),$$

(23)

parameterization

$$\begin{array}{cc}\hfill z_1& =b_1\left(t\right)=\hfill \\ \hfill & \frac{3(173465063\sqrt{3}-1091281895){\beta }_2(-176651t+52563+26043\sqrt{3})}{70130447(2376102210\sqrt{3}t-7846989697t^2-814835898\sqrt{3}+4779544314t-1425469788)},\hfill \\ \hfill z_2& =b_2\left(t\right)=-\frac{\left(2662513+1729681\sqrt{3}\right){\beta }_2^2}{1891308\left(-45683t+17823+9267\sqrt{3}\right)\left(-131t+25+11\sqrt{3}\right)},\hfill \\ \hfill {\beta }_2& =(717+381\sqrt{3}-397t).\hfill \end{array}$$

(24)

The singular points are reached at

$$t_1=\frac{717+381\sqrt{3}}{397}\approx 3.468290574,t_2=\frac{52563+26043\sqrt{3}}{176651}\approx 0.5529026114,$$

The plot of the curve is shown in Figure 6:

According to (8), (10) and (11) in coordinates $A_1,A_2,A_3$, we obtain

$$\begin{array}{cc}\hfill A_1& =1+\frac{(2+\sqrt{3})\sqrt{3}}{12}{z}_1{z}_3+\frac{(-2+\sqrt{3})\sqrt{3}}{12}{z}_2+\frac{1}{2}{z}_3,\hfill \\ \hfill A_2& =1+\frac{\sqrt{3}}{12}{z}_1{z}_3-\frac{\sqrt{3}}{12}{z}_2+\frac{1}{2}{z}_3,A_3=\frac{1}{2}+z_3.\hfill \end{array}$$

(25)

If $z_3=0$ then,

$$A_1=1+\frac{\left(-2+\sqrt{3}\right)\sqrt{3}}{12}{z}_2,A_2=1-\frac{\sqrt{3}}{12}{z}_2,A_3=\frac{1}{2}.$$

(26)

This is the singular line of [23] (Figure 5), which is obtained from the singular line $\mathfrak{I}$ by rotating by an angle $2\pi /3.$ If $z_3=-0.05$, then according to (25) $A_3=9/20=0.45$ and

$$\begin{array}{cc}\hfill A_1& =-\frac{\left(2+\sqrt{3}\right)\sqrt{3}}{240}{z}_1+\frac{\left(-2+\sqrt{3}\right)\sqrt{3}}{12}{z}_2+\frac{39}{40},\hfill \\ \hfill A_2& =-\frac{\sqrt{3}}{240}{z}_1-\frac{\sqrt{3}}{12}{z}_2+\frac{39}{40}.\hfill \end{array}$$

(27)

When parameterized by (24), this curve is shown in Figure 7. It is similar to the part of [23] (Figure 12) corresponding to $A_3=0.45$, with parts of branches $F_1^{-}$ and $G_2^{-}$.

If $z_3=0.005,$ then $A_3=\frac{101}{200}=0.505$ and

$$\begin{array}{cc}\hfill A_1& =-\frac{\left(2+\sqrt{3}\right)\sqrt{3}}{2400}{z}_1+\frac{\left(-2+\sqrt{3}\right)\sqrt{3}}{12}{z}_2+\frac{401}{400},\hfill \\ \hfill A_2& =\frac{\sqrt{3}}{2400}{z}_1-\frac{\sqrt{3}}{12}{z}_2+\frac{401}{400}.\hfill \end{array}$$

(28)

When parameterized by (24), this curve is shown in Figure 8. According to (23) in this figure, when $z_2<0$ and $z_2>\frac{6(\sqrt{3}-1)}{2-\sqrt{3}}$, 2 branches each merge into one line. These results are similar to those of Section 9 of [1], differing from them by the angle $2\pi /3$ rotation.

Figure 8 is similar to the part of Figure 15 in [23] corresponding to $A_3=0.505$ with parts of branches $F_1^{+}$ and $G_2^{+}$.

Here the branches are very close and a different scale is needed. We can draw each branch separately, as in Figure 9.

In that case also there exists a parametric expansions. In the polynomial $V\left(\mathbf{y}\right)=S\left(\mathbf{x}\right)$ we make the power transformation (21) and obtain

$$V\left(\mathbf{y}\right)=W\left(\mathbf{z}\right)=z_3^2\sum _{k=0}^{10}{w}_k\left(z_1,z_2\right)z_3^k.$$

Here $w_0(z_1,z_2)=g(z_1,z_2)$ from (22). According to (24) we substitute (17) into polynomials $w_k(z_1,z_2)$. We obtain $W\left(\mathbf{z}\right)/z_3^2=u_1(\epsilon ,z_3)$ with coefficients depending on t via $b_1\left(t\right)$ and $b_2\left(t\right)$. We apply Theorem 1 in [1] to the equation $u_1(\epsilon ,z_3)=0$ and obtain the expansion

$$\epsilon =\sum _{k=1}^{\infty }{\tilde{c}}_k\left(t\right)z_3^k.$$

So according to (25)

$$\begin{array}{cc}\hfill A_1& =1+\frac{(2\sqrt{3}+3)}{12}(b_1\left(t\right)+\epsilon )z_3+\frac{(3-2\sqrt{3})}{12}(b_2\left(t\right)+\epsilon )+\frac{1}{2}{z}_3,\hfill \\ \hfill A_2& =1+\frac{\sqrt{3}}{12}(b_1\left(t\right)+\epsilon )z_3-\frac{\sqrt{3}}{12}(b_2\left(t\right)+\epsilon )+\frac{1}{2}{z}_3,A_3=\frac{1}{2}+z_3.\hfill \end{array}$$

(29)

Formulas (26) and (27) give the initial terms of them.

2.5. The Normal $(0,-1,-1)$

It corresponds to a truncated polynomial

$$ftr39917=18\left(26+15\alpha \right){\left(y_1-3+3\alpha \right)}^2{\left(y_1-6+6\alpha \right)}^2{y}_1^2·$$

$$\left(18y_1^2{y}_3^2\alpha -6\alpha y_1{y}_2^2+108y_1{y}_3^2\alpha -y_1^2{y}_2^2+36y_1^2{y}_3^2+18y_2^2\alpha +6y_1{y}_2^2+108y_1{y}_3^2-36y_2^2\right).$$

Let’s denote

$$f_3=18y_1^2{y}_3^2\alpha -6\alpha y_1{y}_2^2+108y_1{y}_3^2\alpha -y_1^2{y}_2^2+36y_1^2{y}_3^2+18y_2^2\alpha +6y_1{y}_2^2+108y_1{y}_3^2-36y_2^2.$$

According to [11], we compute the matrix $\gamma =\left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& -1& 1\end{array}\right)$ such that $\left(0,-1,-1\right)\gamma =(0,0,-1)$. Since ${\gamma }^{-1}=\left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 1& 1\end{array}\right)$, then we do a power transformation.

$$y_1=z_1,y_2=z_2{z}_3,y_3=z_3$$

and we get

$$f_3=\left(-6\alpha z_1{z}_2^2-z_1^2{z}_2^2+18\alpha z_1^2+18\alpha z_2^2+6z_1{z}_2^2+108\alpha z_1+36z_1^2-36z_2^2+108z_1\right)z_3^2.$$

Let us denote

$$g_1\left(z_1,z_2\right)=-z_1^2{z}_2^2+\left(18\alpha +36\right)z_1^2+\left(-6\alpha +6\right)z_1{z}_2^2+\left(108\alpha +108\right)z_1+\left(18\alpha -36\right)z_2^2=$$

$$=-z_1^2{z}_2^2+6\left(-\alpha +1\right)z_1{z}_2^2+18\left(\alpha +2\right)z_1^2+18\left(\alpha -2\right)z_2^2+108\left(\alpha +1\right)z_1.$$

According to (22) the polynomial $g_1\left(z_1,z_2\right)$ is obtained from the polynomial $g\left(z_1,z_2\right)$ by the transposition

$$\begin{array}{c}\hfill (z_1,z_2,\alpha )\to (z_2,z_1,-\alpha ).\end{array}$$

(30)

So here everything is similar to Section 2.4, but in the plane $A_1$, $A_2$ rotated by an angle of $2\pi /3$. Also here there are expansions symmetric to (29) by transpositions (30) and

$$\begin{array}{c}\hfill ((b_1\left(t\right),b_2\left(t\right))\to ((b_2\left(t\right),b_1\left(t\right))\end{array}$$

(31)

• Result of Section 2

Theorem 1.

Near the singular point $P_5^{\left(3\right)}$ the variety Ω has 3 local singular parametric expansions (18) and (29) and symmetric to (29) by transpositions (30) and (31). The expansions (18) describe parts of branches $F_1$ and $G_2$ very near the point $P_5^{\left(3\right)}$. Expansions (29) and its symmetric describes branches $F_1$ and $G_2$ near line (26) and symmetry line to it. For $\epsilon \to 0$ branches $F_1$ and $G_2$ coincide at parts of these lines.

3. The Structure of the Manifold $\Omega$ near the Curve $\mathcal{H}$ of Singular Points

Recall that the curve $\mathcal{H}$ is given by equations

$$\begin{array}{cc}\hfill A_1& =b_1\left(t\right)=-\frac{9t\left(t^2+1\right)}{4\left(t+1\right)\left(t^2-4t+1\right)},\hfill \\ \hfill A_2& =b_2\left(t\right)=-\frac{9\left(t^2+1\right)}{4\left(t+1\right)\left(t^2-4t+1\right)},\hfill \\ \hfill A_3& =\frac{3}{4}.\hfill \end{array}$$

(32)

In the polynomial $R\left(\mathbf{A}\right)=Q\left(\mathbf{s}\right)$, substitute

$$A_3=\frac{3}{4}+\mu$$

and write the result as

$$\tilde{R}\left(A_1,A_2,\mu \right)=\sum _{k=0}^m{R}_k(A_1,A_2){\mu }^k.$$

(33)

The polynomials $R_0$, $R_1$, $R_2$, and $R_3$ are computed using the command coeff(R, mu, k) [27]. After factorization, the polynomials $R_0$, $R_1$, $R_2$ and $R_3$ have the form:

$$\begin{array}{c}\hfill R_0\left(A_1,A_2\right)=-(256A_1^6-1536A_1^5{A}_2+768A_1^4{A}_2^2+5120A_1^3{A}_2^3+768A_1^2{A}_2^4-1536A_1{A}_2^5+256A_2^6-5184A_1^4\\ \hfill -10368A_1^2{A}_2^2-5184A_2^4+17496A_1^2+17496A_2^2-19683){(4A_1^3-12A_1^2{A}_2-12A_1{A}_2^2+4A_2^3+9A_1^2+9A_2^2)}^2/2125764,\end{array}$$

$$\begin{array}{c}\hfill R_1\left(A_1,A_2\right)=(-1/177147)\left(8(4A_1^3-12A_1^2{A}_2-12{A_{1}A}_2^2+4A_2^3+9A_1^2+9A_2^2)\right(64A_1^8-384A_1^7{A}_2\\ \hfill +256A_1^6{A}_2^2+896A_1^5{A}_2^3+384A_1^4{A}_2^4+896A_1^3{A}_2^5+256A_1^2{A}_2^6-384{A_{1}A}_2^7+64A_2^8-576A_1^7\\ \hfill +1728A_1^6{A}_2+576A_1^5{A}_2^2+2880A_1^4{A}_2^3+2880A_1^3{A}_2^4+576A_1^2{A}_2^5+1728A_1{A}_2^6-576A_2^7-2592A_1^6\\ \hfill -7776A_2^2{A}_1^4-7776A_2^4{A}_1^2-2592A_2^6+3888A_1^5-11664A_1^4{A}_2-7776A_1^3{A}_2^2-7776A_1^2{A}_2^3\\ \hfill -11664{A_{1}A}_2^4+3888A_2^5+13122A_1^4+26244A_1^2{A}_1^2+13122A_2^4-6561A_1^3+19683A_1^2{A}_2\\ \hfill +19683{A_{1}A}_2^2-6561A_2^3-19683A_1^2-19368A_2^2\left)\right),\end{array}$$

$$\begin{array}{c}\hfill R_2\left(A_1,A_2\right)=-\frac{544}{81}{A}_2^6-\frac{16}{3}{A}_1^3-\frac{544}{81}{A}_1^6-\frac{16}{3}{A}_2^3-\frac{512}{27}{A}_2^5{A}_1+16A_1^2{A}_2+16A_1{A}_2^2\\ \hfill -\frac{512}{27}{A}_2{A}_1^5-\frac{544}{27}{A}_2^2{A}_1^4+\frac{5120}{81}{A}_2^3{A}_1^3-\frac{544}{27}{A}_2^4{A}_1^2+\frac{5120}{19683}{A}_1^2{A}_2^8+\frac{10240}{19683}{A}_1^4{A}_2^6+\frac{8192}{19683}{A}_1^3{A}_2^7\\ \hfill +\frac{5120}{19683}{A}_1^8{A}_2^2+\frac{8192}{19683}{A}_1^7{A}_2^3+\frac{10240}{19683}{A}_1^6{A}_2^4+\frac{28672}{19683}{A}_1^5{A}_2^5-\frac{2048}{6561}{A}_1^9{A}_2-\frac{7168}{2187}{A}_1^2{A}_2^6+\frac{256}{27}{A}_1^2{A}_2^5\\ \hfill +\frac{3584}{729}{A}_1^7{A}_2-\frac{7168}{2187}{A}_1^6{A}_2^2+\frac{256}{9}{A}_1^6{A}_2-\frac{25088}{2187}{A}_1^5{A}_2^3+\frac{256}{27}{A}_1^5{A}_2^2-\frac{3584}{729}{A}_1^4{A}_2^4+\frac{1280}{27}{A}_1^4{A}_2^3\\ \hfill -\frac{688}{9}{A}_1^4{A}_2-\frac{25088}{2187}{A}_1^3{A}_2^5+\frac{1280}{27}{A}_1^3{A}_2^4-\frac{1376}{27}{A}_1^3{A}_2^2+\frac{160}{3}{A}_1^2{A}_2^2-\frac{2048}{6561}{A}_1{A}_2^9+\frac{3584}{729}{A}_1{A}_2^7\\ \hfill +\frac{256}{9}{A}_1{A}_2^6-\frac{688}{9}{A}_1{A}_2^4-\frac{1376}{27}{A}_1^2{A}_2^3-\frac{16384}{6561}{A}_1^8{A}_2+\frac{14336}{6561}{A}_1^7{A}_2^2-\frac{2048}{19683}{A}_1^6{A}_2^3-\frac{75776}{6561}{A}_1^5{A}_2^4\\ \hfill -\frac{75776}{6561}{A}_1^4{A}_2^5-\frac{2048}{19683}{A}_1^3{A}_2^6+\frac{14336}{6561}{A}_1^2{A}_2^7-\frac{16384}{6561}{A}_2^8{A}_1-\frac{256}{27}{A}_2^7+\frac{688}{27}{A}_2^5+\frac{80}{3}{A}_2^4\\ \hfill +\frac{10240}{19683}{A}_1^{10}-\frac{1792}{2187}{A}_2^8+\frac{80}{3}{A}_1^4+\frac{688}{27}{A}_1^5-\frac{256}{27}{A}_1^7-\frac{1792}{2187}{A}_1^8+\frac{1024}{19683}{A}_2^{10}\\ \hfill +\frac{34816}{59049}{A}_2^9+\frac{34816}{59049}{A}_1^9,\end{array}$$

$$\begin{array}{c}\hfill R_3\left(A_1,A_2\right)=\frac{5120}{729}{A}_2^6-\frac{1216}{27}{A}_1^3+\frac{5120}{729}{A}_1^6-\frac{1216}{27}{A}_2^3-\frac{5120}{243}{A}_2^5{A}_1+\frac{1216}{9}{A}_1^2{A}_2+\frac{1216}{9}{A}_2^2{A}_1\\ \hfill -\frac{5120}{243}{A}_1^5{A}_2+\frac{5120}{243}{A}_2^2{A}_1^4+\frac{51200}{729}{A}_1^3{A}_2^3+\frac{5120}{243}{A}_2^4{A}_1^2-\frac{163840}{19683}{A}_1^2{A}_2^6+\frac{23552}{2187}{A}_1^2{A}_2^5+\frac{8192}{6561}{A}_1^7{A}_2\\ \hfill -\frac{163840}{19683}{A}_1^6{A}_2^2+\frac{23552}{729}{A}_1^6{A}_2-\frac{57344}{19683}{A}_1^5{A}_2^3+\frac{23552}{2187}{A}_1^5{A}_2^2-\frac{81920}{6561}{A}_1^4{A}_2^4+\frac{117760}{2187}{A}_1^4{A}_2^3-\frac{14272}{81}{A}_1^4{A}_2\\ \hfill -\frac{57344}{19683}{A}_1^3{A}_2^5+\frac{117760}{2187}{A}_1^3{A}_2^4-\frac{28544}{243}{A}_1^3{A}_2^2+\frac{1280}{27}{A}_1^2{A}_2^2+\frac{8192}{6561}{A_{1}A}_2^7+\frac{23552}{729}{A_{1}A}_2^6-\frac{14272}{81}{A_{1}A}_2^4\\ \hfill -\frac{28544}{243}{A}_1^2{A}_2^3-\frac{65536}{59049}{A}_1^8{A}_2+\frac{57344}{59049}{A}_1^7{A}_2^2-\frac{8192}{177147}{A}_1^6{A}_2^3-\frac{303104}{59049}{A}_1^5{A}_2^4-\frac{303104}{59049}{A}_1^4{A}_2^5-\frac{8192}{177147}{A}_1^3{A}_2^6\\ \hfill +\frac{57344}{59049}{A}_1^2{A}_2^7-\frac{65536}{59049}{A}_1{A}_2^8-\frac{23552}{2187}{A}_2^7+\frac{14272}{243}{A}_2^5+\frac{640}{27}{A}_2^4-\frac{40960}{19683}{A}_2^8+\frac{640}{27}{A}_1^4\\ \hfill +\frac{14272}{243}{A}_1^5-\frac{23552}{2187}{A}_1^7-\frac{40960}{19683}{A}_1^8+\frac{139264}{531441}{A}_2^9+\frac{139264}{531441}{A}_1^9.\end{array}$$

Let’s denote

$$\Phi \left(A_1,A_2\right)=4\left(A_1^3+A_2^3\right)-12\left(A_1^2{A}_2+A_1{A}_2^2\right)+\left(A_1^2+A_2^2\right).$$

Then $R_0$ is divided by ${\Phi }^2$, $R_1$ is divided by $\Phi$, and $R_2$ and $R_3$ are not divided by $\Phi$. The curve $\Phi \left(A_1,A_2\right)=0$ has genus 0, its parameterization is

$$\left\{A_1,A_2\right\}=\{-\frac{9t\left(t^2+1\right)}{4\left(t+1\right)\left(t^2-4t+1\right)},-\frac{9\left(t^2+1\right)}{4\left(t+1\right)\left(t^2-4t+1\right)}\}$$

(34)

and is shown in Figure 10.

The parameterization (34) is the same as the formulas of (32). The curve $\mathcal{H}$ goes to infinity at

$$t_1=-1,t_2=2+\sqrt{3}\approx 3.732050808,t_3=2-\sqrt{3}\approx 0.267949192.$$

According to (32) we substitute

$$A_1=b_1\left(t\right)+\epsilon ,A_2=b_2\left(t\right).$$

(35)

into the polynomials $R_k(A_1,A_2)$. Then the polynomial (33) will become a polynomial

$$\tilde{R}\left(A_1,A_2,\mu \right)=u\left(\epsilon ,\mu \right)=\sum u_{pq}\left(t\right){\epsilon }^p{\mu }^q,$$

(36)

whereby

$$u_{pq}=\frac{1}{p!}·\frac{{\partial }^p{R}_q}{\partial A_1^p},$$

where $A_1=b_1\left(t\right)$, $A_2=b_2\left(t\right)$ according to (32). In particular, we obtain

$$u_{00}=u_{10}=u_{01}\equiv 0,$$

(37)

$$\begin{array}{c}\hfill u_{20}\left(t\right)=\frac{1}{2}·\frac{{\partial }^2{R}_0}{\partial A_1^2}=-\frac{28160}{59049}{A}_1^9-\frac{9728}{59049}{A}_2^9+\frac{320}{243}{A}_1^8+\frac{512}{81}{A}_1^7+\frac{3}{2}{A}_2^2-\frac{112}{9}{A}_1^5-\frac{14}{9}{A}_2^4-\frac{70}{9}{A}_1^4+\frac{16}{27}{A}_2^5+\frac{9}{2}{A}_1^2\\ \hfill +\frac{320}{2187}{A}_2^8-\frac{4}{3}{A}_2^3+\frac{64}{243}{A}_2^6+\frac{20}{3}{A}_1^3+\frac{448}{243}{A}_1^6+\frac{80}{3}{A}_1^4{A}_2-\frac{5120}{243}{A}_1^4{A}_2^3+\frac{3200}{729}{A}_1^4{A}_2^4+\frac{3584}{729}{A}_1^5{A}_2^3\\ \hfill -\frac{28}{3}{A}_1^2{A}_2^2+\frac{160}{9}{A}_1^2{A}_2^3-\frac{512}{81}{A}_1^2{A}_2^5+\frac{1280}{729}{A}_1^2{A}_2^6+\frac{160}{27}{A}_1^3{A}_2^2-\frac{2560}{243}{A}_1^3{A}_2^4+\frac{17920}{2187}{A}_1^3{A}_2^5-\frac{3584}{243}{A}_1^6{A}_2\\ \hfill +\frac{8960}{2187}{A}_1^6{A}_2^2+\frac{716800}{59049}{A}_1^6{A}_2^4-\frac{512}{81}{A}_1^7{A}_2-\frac{16384}{59049}{A}_1^7{A}_2^3-\frac{71680}{19683}{A}_1^8{A}_2^2+\frac{80}{9}{A}_1{A}_2^4-\frac{81920}{59049}{A}_1^3{A}_2^7\\ \hfill -\frac{839680}{59049}{A}_1^4{A}_2^6+\frac{51200}{19683}{A}_1^2{A}_2^8-\frac{57344}{19683}{A}_1^5{A}_2^5-4A_1{A}_2^2-\frac{448}{81}{A}_1{A}_2^5-12A_1^2{A}_2+\frac{64}{27}{A}_1^2{A}_2^4\\ \hfill -\frac{4480}{243}{A}_1^3{A}_2^3+\frac{448}{27}{A}_1^5{A}_2+\frac{320}{81}{A}_1^4{A}_2^2+\frac{71680}{19683}{A}_1^3{A}_2^6+\frac{2048}{729}{A}_1^2{A}_2^7-\frac{3584}{6561}{A}_1{A}_2^8+\frac{35840}{6561}{A}_1^4{A}_2^5-\frac{38912}{6561}{A}_1^7{A}_2^2\\ \hfill -\frac{100352}{19683}{A}_1^6{A}_2^3+\frac{7168}{729}{A}_1^5{A}_2^4+\frac{225280}{177147}{A}_1^9{A}_2+\frac{2560}{729}{A}_1^8{A}_2+\frac{512}{729}{A}_1{A}_2^7-\frac{1024}{243}{A}_1{A}_2^6-\frac{80}{9}{A}_1^3{A}_2\\ \hfill +\frac{80}{9}{A}_1{A}_2^3-\frac{4096}{177147}{A}_1{A}_2^9-\frac{22528}{177147}{A}_1^{10}-\frac{14336}{177147}{A}_2^{10}\\ \hfill =\frac{243{\left(t^2+1\right)}^2{(t^4+6t^2-8t-3)}^2{(5t^6+24t^5+15t^4-80t^3+15t^2+24t+5)}^2}{1024{(t+1)}^8{(t^2-4t+1)}^8},\end{array}$$

$$\begin{array}{c}\hfill u_{11}\left(t\right)=\frac{\partial R_1}{\partial A_1}=32A_1{A}_2^2-\frac{224}{9}{A}_1^5+\frac{71680}{6561}{A}_1^4{A}_2^5+\frac{7168}{243}{A}_1^6{A}_2+32A_1^3-\frac{512}{243}{A}_2^8\\ \hfill -\frac{22528}{177147}{A}_1^{10}-\frac{56}{3}{A}_2^4-\frac{64}{9}{A}_2^5+\frac{512}{81}{A}_1^8+\frac{1024}{243}{A}_2^7+\frac{2048}{729}{A}_1^7-\frac{320}{9}{A}_1^4{A}_2-\frac{35840}{729}{A}_1^4{A}_2^3\\ \hfill -\frac{5120}{243}{A}_1^4{A}_2^4+\frac{2048}{243}{A}_1^5{A}_2^2+\frac{2048}{19683}{A}_2^{10}+\frac{280}{9}{A}_1^4+\frac{640}{9}{A}_1^2{A}_2^3-\frac{8192}{243}{A}_1^5{A}_2^3-\frac{112}{3}{A}_1^2{A}_2^2\\ \hfill +\frac{2048}{243}{A}_1^3{A}_2^4-\frac{4096}{243}{A}_1^3{A}_2^5+\frac{28672}{6561}{A}_1^6{A}_2^4-\frac{38912}{19683}{A}_1^8{A}_2^2-\frac{4096}{243}{A}_1^7{A}_2-\frac{224}{9}{A}_1{A}_2^4\\ \hfill -\frac{114688}{59049}{A}_1^7{A}_2^3+\frac{143360}{59049}{A}_1^4{A}_2^6-\frac{14336}{19683}{A}_1^2{A}_2^8+\frac{57344}{19683}{A}_1^5{A}_2^5+\frac{16384}{6561}{A}_1^3{A}_2^7+\frac{6400}{81}{A}_1^3{A}_2^3\\ \hfill +\frac{1600}{27}{A}_1^2{A}_2^4+\frac{640}{81}{A}_1{A}_2^5+\frac{4096}{2187}{A}_1^2{A}_2^7+\frac{1600}{81}{A}_1^4{A}_2^2+\frac{10240}{2187}{A}_1^5{A}_2^4+\frac{28672}{6561}{A}_1^6{A}_2^3\\ \hfill +\frac{20480}{6561}{A}_1^7{A}_2^2+\frac{20480}{6561}{A}_1^3{A}_2^6-\frac{1024}{243}{A}_1^8{A}_2-\frac{77824}{177147}{A}_1{A}_2^9+\frac{20480}{19683}{A}_1^9{A}_2+\frac{2048}{729}{A}_1{A}_2^6\\ \hfill +\frac{5120}{6561}{A}_1{A}_2^8+\frac{640}{9}{A}_1^5{A}_2-\frac{224}{3}{A}_1^3{A}_2-\frac{224}{9}{A}_1{A}_2^3+\frac{5120}{6561}{A}_1^9-\frac{1024}{2187}{A}_2^9-\frac{2240}{81}{A}_1^6\\ \hfill +\frac{320}{27}{A}_2^6-\frac{4096}{243}{A}_1^2{A}_2^6-\frac{448}{9}{A}_1^3{A}_2^2-\frac{7168}{243}{A}_1^2{A}_2^5=\\ \hfill =\frac{729{\left(t^2+1\right)}^4\left(t^4+6t^2-8t-3\right){\left(5t^6+24t^5+15t^4-80t^3+15t^2+24t+5\right)}^2}{512{(t+1)}^8{(t^2-4t+1)}^8},\end{array}$$

(38)

$$\begin{array}{c}{u}_{02}\left(t\right)=R_2\left(b_1\left(t\right),b_2\left(t\right)\right)=\hfill \\ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}=-\frac{243{\left(t^2+1\right)}^3\left(t^2+8t+1\right)\left(11t^4+8t^3-42t^2+8t+11\right){\left(5t^6+24t^5+15t^4-80t^3+15t^2+24t+5\right)}^2}{1024{\left(t+1\right)}^8{\left(t^2-4t+1\right)}^8}.\hfill \end{array}$$

From the Formulas (37) and (38), we can see that the Newton’s polygon $\Gamma$ of the polynomial $u\left(\epsilon ,\mu \right)$ (36) in the plane $p,q$ has an edge ${\Gamma }_1^{\left(1\right)}$ containing the points $(2,0)$, $(1,1)$, $(0,2)$ (Figure 11) with external normal $N_1=\left(-1,-1\right)$.

So we’re doing a power transformation

$$\epsilon =\mu \delta .$$

(39)

Then the polynomial $u\left(\epsilon ,\mu \right)$ becomes a polynomial.

$${\mu }^{2}V\left(\delta ,\mu \right)={\mu }^2\sum V_{pr}\left(t\right){\delta }^p{\mu }^r=u\left(\epsilon ,\mu \right),$$

(40)

where

$$u_{p,q}\left(t\right)=V_{p,p+q-2}\left(t\right).$$

The support and Newton’s polygon for the polynomial $V\left(\delta ,\mu \right)$ are shown in Figure 12.

The truncated equation corresponding to the edge ${\tilde{\Gamma }}_1^{\left(1\right)}$ is

$$u_{20}\left(t\right){\delta }^2+u_{11}\left(t\right)\delta +u_{02}\left(t\right)=0.$$

(41)

It has two roots

$$\begin{array}{cc}\hfill {\delta }_1\left(t\right)& =\frac{-3{(t^2+1)}^2+2\sqrt{5t^8+24t^7+20t^6-56t^5+30t^4-56t^3+20t^2+24t+5}}{t^4+6t^2-8t-3},\hfill \\ \hfill {\delta }_2\left(t\right)& =\frac{-3{\left(t^2+1\right)}^2-2\sqrt{5t^8+24t^7+20t^6-56t^5+30t^4-56t^3+20t^2+24t+5}}{t^4+6t^2-8t-3}.\hfill \end{array}$$

(42)

Their denominator goes to zero at two real points

$$t_4\approx 1.324070992,t_5\approx -0.3044219351.$$

(43)

These points are indicated in Figure 10.

Next, we consider the expansions of the $\Omega$ manifold for the cases of two roots (42).

In the polynomial $V\left(\delta ,\mu \right)$ of (40), we make the substitutions

$$\delta ={\delta }_i+{\varkappa }_i,i=1,2.$$

(44)

where ${\delta }_i$ are given by the formula (42). We get

$$W_i\left({\varkappa }_i,\mu \right)=V\left(\delta ,\mu \right)=\sum W_{isr}\left(t\right){\varkappa }_i^s{\mu }^r,i=1,2,$$

where integers $s,r\ge 0,$ $r+s\ge 1$. In this case.

$$W_{isr}\left(t\right)=\sum _{p\ge s}{V}_{ps}\left(t\right)C_p^s{\delta }_i^{p-s},i=1,2,$$

where $C_p^s={\displaystyle \frac{p!}{s!(p-s)!}}$ are binomial coefficients. In particular, according to (38), (41) and (42), we have

$$W_{i00}\equiv 0,W_{i10}\equiv 2u_{20}\left(t\right)·{\delta }_i\left(t\right)+u_{11}\left(t\right),i=1,2.$$

More specifically,

$$\begin{array}{cc}\hfill W_{110}& =2u_{20}\left(t\right){\delta }_1\left(t\right)+u_{11}\left(t\right)\hfill \\ \hfill & =\frac{243\left(t^4+6t^2-8t-3\right){\left(t^2+1\right)}^2{\left(5t^6+24t^5+15t^4-80t^3+15t^2+24t+5\right)}^2}{256{\left(t+1\right)}^8}\hfill \\ \hfill & \times \frac{\sqrt{\left(t^2+1\right)\left(5t^6+24t^5+15t^4-80t^3+15t^2+24t+5\right)}}{{\left(t^2-4t+1\right)}^8},\hfill \end{array}$$

(45)

$$\begin{array}{cc}\hfill W_{210}& =2u_{20}\left(t\right){\delta }_2\left(t\right)+u_{11}\left(t\right)\hfill \\ \hfill & =-\frac{243\left(t^4+6t^2-8t-3\right){\left(t^2+1\right)}^2{\left(5t^6+24t^5+15t^4-80t^3+15t^2+24t+5\right)}^2}{256{\left(t+1\right)}^8}\hfill \\ \hfill & \times \frac{\sqrt{\left(t^2+1\right)(5t^6+24t^5+15t^4-80t^3+15t^2+24t+5)}}{{(t^2-4t+1)}^8},\hfill \end{array}$$

(46)

i.e., $W_{210}=-W_{110}\not\equiv 0$. Hence Theorem 1 in [1] is applicable, which for solutions of equations $W_i\left({\varkappa }_i,\mu \right)=0$ according to (45) and (46) gives the expansions

$${\varkappa }_i\left(t\right)=\sum _{k=1}^{\infty }{c}_{ik}\left(t\right){\mu }^k,i=1,2.$$

Let’s go from coordinates $\left({\varkappa }_i,\mu \right)$ to coordinates $\left(\delta ,\mu \right)$ by (44), with decompositions

$${\delta }^{\left(1\right)}={\delta }_1\left(t\right)+{\varkappa }_1\left(t\right)={\delta }_1+\sum _{k=1}^{\infty }{c}_{1k}{\mu }^k$$

and

$${\delta }^{\left(2\right)}={\delta }_2\left(t\right)+{\varkappa }_2\left(t\right)={\delta }_2+\sum _{k=1}^{\infty }{c}_{2k}{\mu }^k.$$

Then we go to the coordinates $\epsilon ,\mu :$ ${\epsilon }^{\left(i\right)}=\mu {\delta }^{\left(i\right)},i=1,2$ by (39) to the coordinates $A_1,\mu$ by (35):

$$A_1^{\left(1\right)}=b_1\left(t\right)+{\epsilon }^{\left(1\right)}=b_1\left(t\right)+\mu \left({\delta }_1+{\varkappa }_1\right)=b_1\left(t\right)+{\delta }_1\left(t\right)\mu +\mu \sum _{k=1}^{\infty }{c}_{1k}\left(t\right){\mu }^k,$$

(47)

$$A_1^{\left(2\right)}=b_1\left(t\right)+{\epsilon }^{\left(2\right)}=b_1\left(t\right)+\mu \left({\delta }_2+{\varkappa }_2\right)=b_1\left(t\right)+{\delta }_2\left(t\right)\mu +\mu \sum _{k=1}^{\infty }{c}_{2k}\left(t\right){\mu }^k,$$

(48)

$$A_2=b_2\left(t\right)=-\frac{9\left(t^2+1\right)}{4\left(t+1\right)\left(t^2-4t+1\right)}.$$

(49)

With $\mu$ fixed, the Formulas (47) and (49) are defined in the plane $A_1,A_2$ with $A_3=3/4+\mu$ the first curve $K^{\left(1\right)}$, and the Formulas (48) and (49) define there the second curve $K^{\left(2\right)}$. We obtain four curves. Restricting ourselves to the initial terms of the expansions, we draw them. When $\mu =\frac{1}{8}$, the curve $K^{\left(1\right)}$ is

$$\begin{array}{cc}{K}^{\left(1\right)}=\hfill & \left\{A_1^{\left(1\right)}=b_1\left(t\right)+{\delta }_1\left(t\right)\mu ,A_2=b_2\left(t\right)\right\}\hfill \\ & =\left\{-\frac{9t\left(t^2+1\right)}{4\left(t+1\right)\left(t^2-4t+1\right)}+\frac{-3{\left(t^2+1\right)}^2+2\sqrt{\left(t^2+1\right)(5t^6+24t^5+15t^4-80t^3+15t^2+24t+5)}}{8·\left(t^4+6t^2-8t-3\right)},b_2\left(t\right)\right\},\hfill \end{array}$$

and it is shown in Figure 13.

The curve $K^{\left(2\right)}$ is

$$\begin{array}{cc}{K}^{\left(2\right)}=\hfill & \left\{A_1^{\left(2\right)}=b_1\left(t\right)+{\delta }_2\left(t\right)\mu ,A_2=b_2\left(t\right)\right\}\hfill \\ & \hfill =\left\{-\frac{9t\left(t^2+1\right)}{4\left(t+1\right)\left(t^2-4t+1\right)}+\frac{-3{\left(t^2+1\right)}^2-2\sqrt{\left(t^2+1\right)(5t^6+24t^5+15t^4-80t^3+15t^2+24t+5)}}{8·\left(t^4+6t^2-8t-3\right)},b_2\left(t\right)\right\},\end{array}$$

It’s shown in Figure 14.

When $\mu =-1/8$, the curve

$$\begin{array}{cc}{K}^{\left(1\right)}=\hfill & \left\{A_1^{\left(1\right)}=b_1\left(t\right)+{\delta }_1\left(t\right)\mu ,A_2=b_2\left(t\right)\right\}\hfill \\ & =\left\{-\frac{9t\left(t^2+1\right)}{4\left(t+1\right)\left(t^2-4t+1\right)}-\frac{-3{\left(t^2+1\right)}^2+2\sqrt{\left(t^2+1\right)\left(5t^6+24t^5+15t^4-80t^3+15t^2+24t+5\right)}}{8·\left(t^4+6t^2-8t-3\right)},b_2\left(t\right)\right\}\hfill \end{array}$$

is shown in Figure 15.

Figure 15. The curve $K^{\left(1\right)}$ at $\mu =-1/8$.

Figure 15. The curve $K^{\left(1\right)}$ at $\mu =-1/8$.

and the curve

$$\begin{array}{cc}{K}^{\left(2\right)}=\hfill & \left\{A_1^{\left(2\right)}=b_1\left(t\right)+{\delta }_2\left(t\right)\mu ,A_2=b_2\left(t\right)\right\}\hfill \\ & =\left\{-\frac{9t\left(t^2+1\right)}{4\left(t+1\right)\left(t^2-4t+1\right)}+\frac{3{\left(t^2+1\right)}^2+2\sqrt{\left(t^2+1\right)\left(5t^6+24t^5+15t^4-80t^3+15t^2+24t+5\right)}}{8·\left(t^4+6t^2-8t-3\right)},b_2\left(t\right)\right\}\hfill \end{array}$$

is shown in Figure 16.

The curves of Figure 13 and Figure 14 correspond to $A_3=7/8$ and are similar to the curves of Figure 13 of [23], showing the section of the variety $\Omega$ by the plane $A_3=1$. The curves of Figure 15 and Figure 16 correspond to $A_3=5/8$ and they are similar to the curves of Figure 14 of [23], showing the cross section of $\Omega$ by the plane $A_3=5/8.$ This confirms the correctness of the calculated expansions.

In Figure 13, Figure 15 and Figure 16, there are discontinuities in the curves at the places of the roots (43) of the denominators ${\delta }_1\left(t\right)$ and ${\delta }_2\left(t\right)$ in (42). They can be eliminated by substituting

$$A_1=b_1\left(t\right)+\epsilon ,A_2=b_2\left(t\right)+\epsilon$$

instead of substitution (35) and calculate the corresponding expansion.

• Result of Section 3

Theorem 2.

Near the curve $\mathcal{H}$ of singular points the variety Ω has two singular parametric expansions (47), (49) and (48), (49). They represent parts of branches $G_2^{+}$ and $G_3^{-}$ correspondingly. At $A_3=3/4$ they coincide with curve $\mathcal{H}$.

4. The Structure of the Variety $\Omega$ near the Curve $\mathcal{F}$ of Singular Points

We take the polynomial $Q\left(\mathbf{s}\right)=Q\left(s_1,s_2,s_3\right)$, where $s_1=a_1+a_2+a_3$, $s_2=a_1·a_2+a_1·a_3+a_1·a_3$, $s_3=a_1·a_2·a_3$ are elementary symmetric polynomials, and we substitute $a_1=a_2$. Then the polynomial $Q\left(\mathbf{s}\right)$ takes the form

$$\tilde{Q}\left(a_1,a_3\right)=-\left(1+2a_3\right)\left(8a_1{a}_3+8a_3^2-4a_1-4a_3+1\right){\left(16a_1^3+16a_1^2{a}_3-4a_1-2a_3+1\right)}^3.$$

Let’s write the polynomial $16a_1^3+16a_1^2{a}_3-4a_1-2a_3+1$ in $\mathbf{A}$ coordinates, substituting

$$a_1=\frac{1+\sqrt{3}}{6}{A}_1+\frac{1-\sqrt{3}}{6}{A}_1+\frac{1}{3}{A}_3,a_3=-\frac{1}{3}{A}_1-\frac{1}{3}{A}_1+\frac{1}{3}{A}_3$$

with $A_1=A_2$.

Then we get a polynomial $-\frac{1}{27}\left(16A_1^3-48A_1{A}_3^2-32A_3^3+54A_3-27\right)$. We put

$$\mathcal{F}\left(A_1,A_3\right)=16A_1^3-48A_1{A}_3^2-32A_3^3+54A_3-27.$$

The curve $\mathcal{F}=0$ consists of singular points, has genus 0, and parameterization

$$\left[A_1,A_3\right]=\left[b_1\left(t\right)=-\frac{\left(5t+2\right){\left(t+4\right)}^2}{6t\left(t^2-16t-8\right)},b_2\left(t\right)=\frac{11t^3-48t^2-48t-16}{6t\left(t^2-16t-8\right)}\right].$$

(50)

The curve is shown in Figure 17.

In [23] (Figure 3), the components $f_1^{\pm }$, $f_2$, $f_3^{\pm }$ of this curve are shown in gray. The scales on the axes are different there. In the polynomial $R\left(\mathbf{A}\right)=Q\left(\mathbf{s}\right)$ we substitute

$$A_1=B_1,A_2=B_1+B_2,{A_3=B}_3,$$

(51)

and we get a polynomial depending on three variables,

$$K\left(B_1,B_2,B_3\right)=\sum _{l=0}^{12}{K}_l(B_1,B_3)B_2^l.$$

(52)

We factorize $K_l$ for $l=0,1,2,3$ because they are need for our computation and get

$$\begin{array}{cc}-531441K_0(B_1,B_3)=(-{2B}_3-3+\hfill & 4B_1)\left(16B_1^2-40B_1{B}_3+16B_3^2+12B_1-24B_3+9\right)\hfill \\ & \times {\left(16B_1^3-48B_1{B}_3^2-32B_3^3+54B_3-27\right)}^3=\hfill \\ & =\left(-{2B}_3-3+4B_1\right)\left(16B_1^2-40B_1{B}_3+16B_3^2+12B_1-24B_3+9\right){\mathcal{F}}^3(B_1,B_3).\hfill \end{array}$$

and

$$-177147K_1\left(B_1,B_3\right)=8B_1{\left(16B_1^3-48B_1{B}_3^2-32B_3^3+54B_3-27\right)}^2$$

$$\times \left(256B_1^4-704B_1^3{B}_3+96B_1^2{B}_3^2+736B_1{B}_3^3-320B_3^4+378B_1{B}_3-432B_3^2-189B_1+216B_3\right)$$

$$\begin{array}{c}\hfill =8B_1\left(256B_1^4-704B_1^3{B}_3+96B_1^2{B}_3^2+736B_1{B}_3^3-320B_3^4+378B_1{B}_3-432B_3^2-189B_1+216B_3\right){\mathcal{F}}^2\left(B_1,B_3\right).\end{array}$$

Then

$$\begin{array}{c}\hfill K_2\left(B_1,B_3\right)=-\frac{45056}{729}{B}_1^7{B}_3+\frac{1856}{27}{B}_1^4{B}_3-\frac{131072}{19683}{B}_1^{10}-\frac{467}{27}{B}_1^4+\frac{22528}{729}{B}_1^7+\frac{1310720}{177147}{B}_3^{10}\\ \hfill -\frac{32768}{2187}{B}_3^8+\frac{16384}{2187}{B}_3^7-\frac{128}{3}{B}_3^3+\frac{2048}{81}{B}_3^4+\frac{1024}{81}{B}_3^5-\frac{1024}{81}{B}_3^6+\frac{64}{3}{B}_3^2-\frac{32}{9}{B}_3-\frac{65536}{729}{B}_1^2{B}_3^5\\ \hfill +\frac{106496}{2187}{B}_1^3{B}_3^4+\frac{136192}{729}{B}_1^4{B}_3^3+\frac{512}{27}{B}_1^2{B}_3^2+\frac{2048}{81}{B}_1^3{B}_3-\frac{212992}{2187}{B}_1^3{B}_3^5-\frac{40960}{729}{B}_1^5{B}_3^2\\ \hfill -\frac{177152}{2187}{B}_1^6{B}_3+\frac{131072}{729}{B}_1^2{B}_3^6+\frac{81920}{729}{B}_1^5{B}_3^3-\frac{272384}{729}{B}_1^4{B}_3^4+\frac{354304}{2187}{B}_1^6{B}_3^2+\frac{3080192}{177147}{B}_1^9{B}_3\\ \hfill -\frac{1638400}{19683}{B}_1^7{B}_3^3-\frac{5472256}{59049}{B}_1^6{B}_3^4+\frac{2621440}{19683}{B}_1^5{B}_3^5+\frac{2768896}{19683}{B}_1^4{B}_3^6+\frac{671744}{19683}{B}_1^8{B}_3^2-\frac{1441792}{19683}{B}_1^2{B}_3^8\\ \hfill -\frac{3407872}{59049}{B}_1^3{B}_3^7+\frac{2048}{27}{B}_1^2{B}_3^4+\frac{8192}{81}{B}_1^3{B}_3^3-\frac{1856}{27}{B}_1^4{B}_3^2-\frac{8192}{81}{B}_1^3{B}_3^2-\frac{2048}{27}{B}_1^2{B}_3^3,\\ \hfill K_3\left(B_1,B_3\right)=\frac{1024}{27}{B}_1{B}_3^4+\frac{65536}{729}{B}_1{B}_3^6-\frac{1024}{27}{B}_1{B}_3^3-\frac{32768}{729}{B}_1{B}_3^5+\frac{256}{27}{B}_1{B}_3^2-\frac{1952}{243}{B}_1^3\\ \hfill +\frac{336896}{6561}{B}_1^6-\frac{3276800}{531441}{B}_1^9-\frac{917504}{531441}{B}_3^9-\frac{90112}{6561}{B}_3^7+\frac{3424}{243}{B}_3^3-\frac{5632}{243}{B}_3^4+\frac{5632}{243}{B}_3^5+\frac{45056}{6561}{B}_3^6\\ \hfill -\frac{112}{9}{B}_3^2+\frac{56}{3}{B}_3-\frac{720896}{19683}{B}_1{B}_3^8-\frac{28}{27}-\frac{851968}{59049}{B}_1^2{B}_3^7+\frac{23953408}{177147}{B}_1^3{B}_3^6+\frac{655360}{6561}{B}_1^4{B}_3^5\\ \hfill -\frac{6324224}{59049}{B}_1^5{B}_3^4-\frac{13991936}{177147}{B}_1^6{B}_3^3+\frac{1015808}{59049}{B}_1^8{B}_3-\frac{53248}{2187}{B}_1^2{B}_3^5-\frac{2772992}{6561}{B}_1^3{B}_3^4+\frac{20480}{243}{B}_1^4{B}_3^3\\ \hfill -\frac{2048}{81}{B}_1^2{B}_3^2+\frac{7808}{243}{B}_1^3{B}_3+\frac{604160}{2187}{B}_1^5{B}_3^2-\frac{673792}{6561}{B}_1^6{B}_3+\frac{720896}{19683}{B}_1^7{B}_3^2+\frac{512}{81}{B}_1^2{B}_3+\frac{26624}{2187}{B}_1^2{B}_3^4\\ \hfill +\frac{1386496}{6561}{B}_1^3{B}_3^3-\frac{10240}{243}{B}_1^4{B}_3^2-\frac{302080}{2187}{B}_1^5{B}_3-\frac{7808}{243}{B}_1^3{B}_3^2+\frac{2048}{81}{B}_1^2{B}_3^3.\end{array}$$

The multiplier $\mathcal{F}(B_1,B_3)$ enters in $K_0\left(B_1,B_3\right)$ in the third degree, in $K_1\left(B_1,B_3\right)$ in the second degree, in $K_2\left(B_1,B_3\right)$, and in $K_3\left(B_1,B_3\right)$ it does not enter. Then $K_0$ is divisible by ${\mathcal{F}}^3$, $K_1$ is divisible by ${\mathcal{F}}^2$, and $K_2$ and $K_3$ are not divisible by $\mathcal{F}$. The curve $\mathcal{F}\left(B_1,B_3\right)=0$ has genus 0, parameterization (50).

The curve $\mathcal{F}=0$ goes to infinity at

$$t_1=0,t_2=2\left(4+3\sqrt{2}\right)\approx 16.48528137,t_3=2\left(4-3\sqrt{2}\right)\approx -0.485281372.$$

Into the polynomials $K_l(B_1,B_3)$ we substitute

$$B_1=b_1\left(t\right)+\epsilon ,B_3=b_2\left(t\right)$$

(53)

according to (50). Then the polynomial (52) will become a polynomial

$$K\left(B_1,B_3,B_2\right)=u\left(\epsilon ,B_2\right)=\sum _{p,q\ge 0}{u}_{pq}\left(t\right){\epsilon }^p{B}_2^q,\mathrm{whereby}{u}_{pq}=\frac{1}{p!}·\frac{{\partial }^p{K}_q}{\partial B_1^p}\left(b_1\left(t\right),b_2\left(t\right)\right),$$

(54)

where $B_1=b_1\left(t\right)$, $B_3=b_2\left(t\right)$ according to (50). In particular, we obtain

$$\begin{array}{cccc}\hfill u_{00}=& u_{10}=u_{01}\equiv 0,\hfill & \hfill u_{20}\left(t\right)=& \frac{1}{2}·\left(\frac{{\partial }^2{K}_0}{\partial B_1^2}\right)\left(b_1\left(t\right),b_2\left(t\right)\right)=0,\hfill \\ \hfill u_{11}\left(t\right)=& \left(\frac{\partial K_1}{\partial B_1}\right)\left(b_1\left(t\right),b_2\left(t\right)\right)=0,\hfill & \hfill u_{02}\left(t\right)=& K_2\left(b_1\left(t\right),b_2\left(t\right)\right)=\frac{64{(5t+2)}^2{(t+4)}^4{(t-2)}^4{(t+1)}^4}{t^2{(t^2-16t-8)}^6}.\hfill \end{array}$$

(55)

$$\begin{array}{c}\hfill u_{30}=\frac{1}{6}·\frac{{\partial }^3{K}_0}{\partial B_1^3}==\frac{40960}{81}{B}_1^2{B}_3^3+\frac{14417920}{59049}{B}_1^8{B}_3+\frac{10240}{81}{B}_1^2{B}_3-\frac{286720}{729}{B}_1^4{B}_3^2-\frac{57671680}{531441}{B}_1^9\\ \hfill +\frac{360448}{6561}{B}_3^7+\frac{4259840}{6561}{B}_1^3{B}_3^3+\frac{3670016}{531441}{B}_3^9+\frac{22528}{243}{B}_3^4-\frac{22528}{243}{B}_3^5-\frac{180224}{6561}{B}_3^6\\ \hfill +\frac{1835008}{2187}{B}_1^5{B}_3^2+\frac{102400}{243}{B}_1^3{B}_3+\frac{2621440}{6561}{B}_1^7{B}_3^2-\frac{2981888}{6561}{B}_1^6{B}_3-\frac{150470656}{177147}{B}_1^6{B}_3^3\\ \hfill -\frac{2883584}{19683}{B}_1{B}_3^8+\frac{4096}{27}{B}_1{B}_3^4-\frac{4096}{27}{B}_1{B}_3^3-\frac{131072}{729}{B}_1{B}_3^5+\frac{262144}{729}{B}_1{B}_3^6+\frac{573440}{729}{B}_1^4{B}_3^3\\ \hfill +\frac{1490944}{6561}{B}_1^6-\frac{224}{9}{B}_3-\frac{13696}{243}{B}_3^3+\frac{448}{9}{B}_3^2-\frac{25600}{243}{B}_1^3+\frac{112}{27}-\frac{917504}{2187}{B}_1^5{B}_3\\ \hfill -\frac{102400}{243}{B}_1^3{B}_3^2-\frac{17039360}{59049}{B}_1^2{B}_3^7+\frac{103546880}{177147}{B}_1^3{B}_3^6+\frac{1024}{27}{B}_1{B}_3^2-\frac{40370176}{59049}{B}_1^5{B}_3^4\\ \hfill -\frac{1064960}{2187}{B}_1^2{B}_3^5+\frac{18350080}{19683}{B}_1^4{B}_3^5+\frac{532480}{2187}{B}_1^2{B}_3^4-\frac{40960}{81}{B}_1^2{B}_3^2-\frac{8519680}{6561}{B}_1^3{B}_3^4\\ \hfill =-\frac{8192\left(5t+2\right){\left(t+4\right)}^2{\left(t-2\right)}^2{\left(t+1\right)}^3{\left(8t^3-3t^2+24t+8\right)}^3}{6561t^5{\left(t^2-16t-8\right)}^6}.\end{array}$$

(56)

From the Formulas (55) and (57), we can see that the Newton’s polygon $\Gamma$ of the polynomial $u\left(\epsilon ,B_2\right)$ given by (54) in the plane $p,q$ has an edge ${\Gamma }_1^{\left(1\right)}$, containing the points $(3,0)$, $(0,2)$ (Figure 18) with external normal $N_1=\left(-2,-3\right)$.

A truncated polynomial corresponds to this edge

$${\epsilon }^3{u}_{30}\left(t\right)+B_2^2{u}_{02}\left(t\right)=0$$

(57)

According to [11] we find the unimodular matrix $\alpha =\left(\begin{array}{cc}3& -1\\ -2& 1\end{array}\right)$ for $N_1$ such that $N_1\alpha =\left(0,-1\right)$. Therefore, we need to do a power transformation

$$\left(ln\delta ,lnD\right)=\left(ln\epsilon ,ln{B}_2\right)·\alpha ,$$

where $\delta$ and D are new variables, i.e.,

$$\left(ln\epsilon ,ln{B}_2\right)=\left(ln\delta ,lnD\right)·{\alpha }^{-1}.$$

Since ${\alpha }^{-1}=\left(\begin{array}{cc}1& 1\\ 2& 3\end{array}\right)$, then $\epsilon =\delta D^2,B_2=\delta D^3$. Hence we can write

$$\begin{array}{c}\hfill K\left(B_1,B_3,B_2\right)=u\left(\epsilon ,B_2\right)=\sum u_{pq}\left(t\right){\epsilon }^p{B}_2^q=\sum u_{pq}\left(t\right){\delta }^{p+q}{D}^{2p+3q}={\delta }^2{D}^{6}V\left(\delta ,D\right).\end{array}$$

Then the polynomial $u\left(\epsilon ,B_2\right)$ becomes a polynomial

$${\delta }^2{D}^{6}V\left(\delta ,D\right)={\delta }^2{D}^6\sum V_{rs}\left(t\right){\delta }^r{D}^s=u\left(\epsilon ,B_2\right),$$

where $V_{r,s}\left(t\right)=V_{p+q,2p+3q}\left(t\right)=u_{p,q}\left(t\right)$. Thus the polygon $\Gamma$ of Figure 18 takes the form shown in Figure 19. For the polynomial $V\left(\delta ,D\right)$ the polygon is shown in Figure 20. The truncated Equation (57) takes the form of

$${\delta }^2{D}^6\left(u_{30}\left(t\right)\delta +u_{02}\left(t\right)\right)=0.$$

From where

$${\delta }_0\left(t\right)=c_0\left(t\right)=-\frac{u_{02}\left(t\right)}{u_{30}\left(t\right)}=\frac{6561\left(5t+2\right){·\left(t+4\right)}^2·{\left(t-2\right)}^2\left(t+1\right)t^3}{128{\left(8t^3-3t^2+24t+8\right)}^3}.$$

The only real root of the denominator is

$$t_4=-\frac{3{\left(13+16\sqrt{2}\right)}^{\frac{1}{3}}}{8}+\frac{21}{8{\left(13+16\sqrt{2}\right)}^{\frac{1}{3}}}+\frac{1}{8}\approx -0.3111842957.$$

(58)

In this case.

$$V_{g,h}\left(t\right)=V_{p+q-2,2p+3q-6}\left(t\right)=u_{p,q}\left(t\right).$$

After substitution $\delta ={\delta }_0\left(t\right)+\xi$, into the polynomial $V\left(\delta ,D\right)$, we obtain

$$W\left(\xi ,D\right)=V\left({\delta }_0\left(t\right)+\xi ,D\right).$$

When $\xi =0$, the polynomial $W\left(0,D\right)$ is calculated using the command $coeff\left(M\left(0,D\right),0\right)$ [27]. The quotient at D of degree zero is zero. The coefficient on the first degree of D is obtained by

$$a\left(t\right)=coeff\left(M\left(0,D\right),D,1\right)=\frac{1594323t^6{\left(5t+2\right)}^2{\left(t+4\right)}^4{\left(t-2\right)}^4{\left(t+1\right)}^4}{2}.$$

Therefore, Theorem 1 in [1] is applied to equation $W\left(\xi ,D\right)=0$, and according to it a solution is

$$\xi =\sum _{k=1}^{\infty }{c}_k\left(t\right)D^k.$$

(59)

When we get the truncated equation $u_{30}\left(t\right)·\xi +a\left(t\right)·D=0$, then it follows

$$\xi =-\frac{a\left(t\right)}{u_{30}\left(t\right)}·D=\frac{10460353203t^{11}\left(t+1\right)\left(5t+2\right){\left(t+4\right)}^2{\left(t-2\right)}^2{\left(t^2-16t-8\right)}^6}{16384{\left(8t^3-3t^2+24t+8\right)}^3}·D=c_1\left(t\right)D.$$

Now let’s go back and get an approximation

$$\epsilon =\delta D^2=({\delta }_0\left(t\right)+\xi )D^2\approx {\delta }_0\left(t\right)D^2+c_1\left(t\right)D^3,$$

(60)

$$B_2=({\delta }_0\left(t\right)+\xi )D^3\approx {\delta }_0\left(t\right)D^3+c_1\left(t\right)D^4.$$

(61)

Therefore, from the formula (51) we get

$$A_1=B_1=b_1\left(t\right)+{\delta }_0\left(t\right)D^2+c_1\left(t\right)D^3,$$

(62)

$$A_2=B_1+B_2=b_1\left(t\right)+{\delta }_0\left(t\right)D^2+(c_1\left(t\right)+{\delta }_0\left(t\right){)D}^3+c_1\left(t\right)D^4,$$

$$A_3=B_3=b_2\left(t\right)=\frac{11t^3-48t^2-48t-16}{6t\left(t^2-16t-8\right)}.$$

(63)

The curves (62) and (63) at $D=\pm 0.1$ are shown in Figure 21 and Figure 22, respectively. The gaps in these curves are the neighborhoods of the point $t_4$ of (58). They can be filled in if instead of substituting (51) we do

$$B_1=b_1\left(t\right)+\epsilon ,B_3=b_2\left(t\right)+\epsilon .$$

The closeness of these curves to the curve of Figure 17 confirms the correctness of the found parametric expansion of the (59) of the variety $\Omega$ near the curve of singular points.

According to (61) and (62) branches $F_i$ intersect curve $\mathcal{F}$ with singularity of type

$$\sqrt{\frac{A_1-b_1\left(t\right)}{{\delta }_0\left(t\right)}}\approx {\left(\frac{B_2}{{\delta }_0\left(t\right)}\right)}^{1/3}.$$

• Result of Section 4

Theorem 3.

Near the curve $\mathcal{F}$ of singular points the variety Ω has one singular parametric expansions (59)–(61) and (63). They represent parts of branches $F_1^{\pm }$, $F_2$, $F_3^{\pm }$. At $A_1=A$ they coincide with curve $\mathcal{F}$, having points of curve $\mathcal{F}$ as singular points.

5. The Variety $\Omega$ at Infinity

The number of branches of the variety $\Omega$ at infinity exceeds their number near its singularities. Their complete study would exceed 7 sections on branches in the finite domain (4 Section in [1] and 3 Section here). Therefore, we study here only those branches corresponding to the first nonlinear polynomial multiplier included in the truncated polynomial in degree one.

5.1. Reducing the Study at Infinity to the Study in the Finite Domain

In the polynomial $S\left(\mathbf{A}\right)=Q\left(\mathbf{s}\right)$, we do a power transformation

$$A_1=B_1{B}_3,A_2=B_2{B}_3,A_3=B_3.$$

(64)

The resulting polynomial is divided by $B_3^{12}$ and factorized, we get

$$S(B_1{B}_3,B_2{B}_3,B_3)/B_3^{12}=-\frac{1}{531441}\tilde{T}\left(B_1,B_2,B_3\right).$$

In the sum $\tilde{T}\left(B_1,B_2,B_3\right)$, which is not a polynomial, we substitute

$$B_1=C_1,B_2=C_2,B_3=C_3^{-1}.$$

(65)

The resulting polynomial is $J\left(C_1,C_2,C_3\right)=\tilde{T}\left(B_1,B_2,B_3\right)$. Let us explain the meaning of these transformations for the two-dimensional case, restricting ourselves to coordinates $A_2$ and $A_3$. The polyhedron of the original polynomial $S\left(\mathbf{A}\right)$ has the form shown in Figure 23.

After replacing (64), it takes the form shown in Figure 24.

The polyhedron of sum $T\left(\mathbf{B}\right)$ is shown in Figure 25.

After substituting (65), we get the polynomial $J\left(\mathbf{C}\right),$ whose polyhedron is shown in Figure 26.

In Figure 23, Figure 24, Figure 25 and Figure 26, the edge that corresponds to infinity in coordinates $\mathbf{A}$ is bolded.

Now we need to study the polynomial $J\left(\mathbf{C}\right)$ at $C_1,C_2⟶const,C_3\to 0$. For this purpose, we compute the Newton’s polyhedron ${\Gamma }_8$ of the polynomial $J\left(\mathbf{C}\right)$. Its graph is shown in Figure 27.

It has 4 two-dimensional faces with external normals

$$N_{53}=\left(1,1,1\right),N_{71}=\left(-1,0,0\right),N_{77}=\left(0,-1,0\right),N_{79}=\left(0,0,-1\right).$$

Since $C_1$ and $C_2\to \mathrm{const}$, and $C_3\to 0$, we select the only normal $N_{79}=\left(0,0,-1\right)$, which has only the third coordinate negative. After factorization, the corresponding truncated polynomial $ftr79$ has the form:

$$\begin{array}{cc}ftr79=\hfill & 1024{\left(C_1^2-4C_1{C}_2+C_2^2-2C_1-2C_2-2\right)}^2\hfill \\ & \times {\left(C_1^2-4C_1{C}_2+C_2^2+4C_1+4C_2-8\right)}^2{\left(C_1+C_2+2\right)}^2{\left(C_1+C_2-1\right)}^2.\hfill \end{array}$$

(66)

We will devote a separate subsection to each of its multipliers.

5.2. The First Multiplier in (66)

Multiplier

$$f_1=C_1^2-4C_1{C}_2+C_2^2-2C_1-2C_2-2$$

does not factorize in the field of rational numbers, but does factorize in the extension of that field with $\alpha =\sqrt{3}$. It is the product of two linear forms $f_1=D_1{D}_2$, and we will consider the whole thing as a coordinate substitution, where

$$D_1=C_1-C_2\left(\sqrt{3}+2\right)-\left(\sqrt{3}+1\right),$$

$$D_2=C_1+C_2\left(\sqrt{3}-2\right)+\left(\sqrt{3}-1\right),$$

$$D_3=C_3.$$

and put $\mathbf{D}=\left(D_1,D_2,D_3\right)$. Its inverse substitution is

$$C_1=\frac{(-2+\sqrt{3})\sqrt{3}}{6}{D}_1+\frac{(2+\sqrt{3})\sqrt{3}}{6}{D}_2-1,$$

$$C_2=-\frac{\sqrt{3}}{6}{D}_1+\frac{\sqrt{3}}{6}{D}_2-1,C_3=D_3.$$

(67)

We substitute it into the polynomial $J\left(\mathbf{C}\right)$ and get the polynomial $S_1\left(\mathbf{D}\right)=J\left(\mathbf{C}\right)$. For the polynomial $S_1\left(\mathbf{D}\right)$, we compute its Newton’s polyhedron ${\Gamma }_9$.

Its graph is shown in Figure 28. It has 11 two-dimensional faces with external normals

$$\begin{array}{c}\hfill N_{5645}=\left(-1,0,-1\right),N_{6101}=\left(-1,-1,-1\right),N_{13631}=\left(0,-1,-1\right),N_{116147}=\left(1,1,1\right),\\ \hfill N_{122463}=\left(-1,0,0\right),N_{122607}=\left(0,3,2\right),N_{124121}=\left(0,1,0\right),N_{133225}=\left(0,0,-1\right),\\ \hfill N_{150921}=\left(3,0,2\right),N_{164051}=\left(1,0,0\right),N_{175461}=\left(0,-1,0\right).\end{array}$$

Since $D_1,$$D_2$ and $D_3\to 0$, we select the only normal that has all coordinates negative. This is $N_{6101}=\left(-1,-1,-1\right)$. It corresponds to the truncated polynomial

$$\begin{array}{c}\hfill ftr6101=-23328\left(-2+\alpha \right)(288\alpha D_1^4{D}_3^2-128\alpha D_1^3{D}_2^3+256\alpha D_1^2{D}_2^4-288\alpha D_1^2{D}_2^2{D}_3^2\\ \hfill +1728\alpha D_1{D}_2^3{D}_3^2-1782\alpha D_1{D}_2{D}_3^4-4320\alpha {D_2^{4}D}_3^2+4212\alpha D_2^2{D}_3^4-2916\alpha D_3^6+64D_1^4{D}_2^2\\ \hfill -576D_1^4{D}_3^2-256D_1^3{D}_2^3+432D_1^3{D}_2{D}_3^2+448D_1^2{D}_2^4-576D_1^2{D}_2^2{D}_3^2+1053D_1^2{D}_3^4\\ \hfill +3024D_1{D}_2^3{D}_3^2-3564D_1{D}_2{D}_3^4-7488D_2^4{D}_3^2+7371D_2^2{D}_3^4-5832D_3^6).\end{array}$$

After the power transformation

$$D_1=M_1{M}_3,D_2=M_2{M}_3,D_3=M_3,$$

(68)

we get $ftr6101=F_{10}(M_1,M_2)·M_3^6$, where

$$\begin{array}{c}\hfill F_{10}\left(M_1,M_2\right)=(-128\alpha M_1^3{M}_2^3+256\alpha M_1^2{M}_2^4+64M_1^4{M}_2^2-256M_1^3{M}_2^3+448M_1^2{M}_2^4+288\alpha M_1^4\\ \hfill -288\alpha M_1^2{M}_2^2+1728\alpha M_1{M}_2^3-4320\alpha M_2^4-576M_1^4+432M_1^3{M}_2-576M_1^2{M}_2^2+3024M_1{M}_2^3\\ \hfill -7488M_2^4-1782\alpha M_1{M}_2+4212\alpha M_2^2+1053M_1^2-3564M_1{M}_2+7371M_2^2-2916\alpha -5832)M_3^6.\end{array}$$

The curve $F_{10}=0$ has genus 0, and parameterization

$$\begin{array}{cc}{M_1=b}_1\left(t\right)=\hfill & -\left(3(35677231+53951067\alpha )\right(226041519229686484434944000\alpha t^3\hfill \\ & \hspace{1em}\hspace{1em}+312865218289492809482240000t^4+29692450221454838768025600\alpha t^2+\hfill \\ & \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}2752150186688110972108800t^3+19470907467358707865979520\alpha t\hfill \\ & \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}+205018459636432060173312000t^2+431078468622082108802982\alpha \hfill \\ & \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}+552278881165356491537664t+2976481999975581125816265\left)\right)/\hfill \\ & \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}(81242985303506296261120t \times \hfill \\ & \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\times (2346560t+357378+757301\alpha )(56960\alpha t+204800t^2-26880t+25527)),\hfill \end{array}$$

(69)

$$\begin{array}{c}\hfill {M_2=b}_2\left(t\right)=-\left((-327921327+145854781\alpha )\right(84568711555214157742080000\alpha t^3\\ \hfill +15278292434425670008832320000t^4-582100343712291111321600\alpha t^2\\ \hfill -8379475635071933625131392000t^3+22597125012513902404916160\alpha t\\ \hfill +129684691757394325122969600t^2+862156937244164217605964\alpha \\ \hfill +7644517497617158949374080t+4294936642585162134439377\left)\right)/\\ \hfill (123980512928512598138\left(56960\alpha t+204800t^2-26880t+25527\right)\\ \hfill \left(56960t+8509\alpha \right)\left(-8509+8960t\right)),\end{array}$$

(70)

and is shown in Figure 29.

In (70), the denominator has 2 real roots

$$t_1=0,t_2\approx -0.7112802609.$$

(71)

In (70), the denominator has 2 real roots

$$t_3=\frac{8509}{8960}\approx 0.9496651786,t_4=-\frac{8509\sqrt{3}}{56960}\approx -0.2587433343$$

(72)

In fact, the parametric expansion of the variety $\Omega$ can also be calculated here. To do this, we perform a power transformation (68) to the polynomial $S_1\left(\mathbf{D}\right)$ and similarly to (16) get the polynomial

$$-\frac{T\left(\mathbf{M}\right)}{M_3^6}=\sum _{k=0}^6{T}_k\left(M_1,M_2\right)M_3^k,$$

where $\mathbf{M}=(M_1,M_2,M_3)$. We substitute into the polynomials $T_k\left(M_1,M_2\right)$ according to (69) and (70)

$$M_1=b_1\left(t\right)+\epsilon ,M_2=b_2\left(t\right)+\epsilon .$$

(73)

We obtain the polynomial $u\left(\epsilon ,M_3\right)=-T\left(M_1,M_2,M_3\right)/M_3^6$ with coefficients depending on t via $b_1\left(t\right)$ and $b_2\left(t\right)$. In this polynomial.

$$u\left(\epsilon ,M_3\right)=\sum _{k=0}^6{T}_k(b_1+\epsilon ,b_2+\epsilon )M_3^k=\sum _{p,q\ge 0}{u}_{pq}{\epsilon }^p{M}_3^q,$$

where

$$u_{pq}=\sum _{p_1+p_2=p\ge 1}\frac{1}{p_1!p_2!}·\frac{{\partial }^p{T}_q}{\partial M_1^{p_1}\partial M_2^{p_2}},$$

(74)

when $M_i=b_i\left(t\right)$, $i=1,2$, $p_1,p_2\ge 0$, $p\ge 1$,

Specifically

$$\begin{array}{cc}\hfill u_{00}\equiv & 0,\hfill \\ \hfill u_{10}=& \frac{\partial T_0\left(M_1,M_2\right)}{\partial M_1}+\frac{\partial T_0\left(M_1,M_2\right)}{\partial M_2}=\hfill \\ \hfill & 2985984\alpha M_1^4{M}_2+5971968\alpha M_1^3{M}_2^2-5971968\alpha M_1^2{M}_2^3-2985984\alpha M_1{M}_2^4\hfill \\ \hfill & -5971968M_1^4{M}_2-2985984M_1^3{M}_2^2-2985984M_1^2{M}_2^3-5971968M_1{M}_2^4\hfill \\ \hfill & -97417728\alpha M_1^3+30233088\alpha M_1^2{M}_2-30233088\alpha M_1{M}_2^2+97417728\alpha M_2^3\hfill \\ \hfill & +167961600M_1^3-47029248M_1^2{M}_2-47029248M_1{M}_2^2+167961600M_2^3\hfill \\ \hfill & +49128768\alpha M_1-49128768\alpha M_2-56687040M_1-56687040M_2\stackrel{\mathrm{def}}{=}H\left(b_1\left(t\right),b_2\left(t\right)\right),\hfill \\ \hfill u_{01}=& T_1\left(b_1\left(t\right),b_2\left(t\right)\right)\stackrel{\mathrm{def}}{=}G\left(b_1\left(t\right),b_2\left(t\right)\right).\hfill \end{array}$$

Here the sign $\stackrel{\mathrm{def}}{=}$ means new notation.

Indeed functions $u_{10}\left(t\right)$ and $u_{01}\left(t\right)$ have very complicated forms. So we omit them and give only some their properties.

The function $u_{10}\left(t\right)$ has two multiple roots

$$\begin{array}{cc}\hfill t_5=& -\frac{5553288233\sqrt{3}}{11245663040}-\frac{1415395569}{2811415760}+\frac{6443874209\sqrt{6}}{22491326080}+\frac{828223515\sqrt{2}}{2249132608}\approx -0.1361976710,\hfill \\ \hfill t_6=& -\frac{5553288233\sqrt{3}}{11245663040}-\frac{1415395569}{2811415760}-\frac{6443874209\sqrt{6}}{22491326080}-\frac{828223515\sqrt{2}}{2249132608}\approx -2.581322779\hfill \end{array}$$

(75)

of multiplicity 6, and the function $u_{01}\left(t\right)$ has the same roots of multiplicity 8. The denominators of the functions $u_{10}\left(t\right)$ and $u_{01}\left(t\right)$ each have four multiple roots of (71) and (72). By the implicit function theorem [1] (Theorem 1), the equation $u\left(\epsilon ,M_3\right)=0$ has a solution as a power series on $M_3$

$$\epsilon =\sum _{k=1}^{\infty }{c}_k\left(t\right)M_3^k,$$

(76)

where $c_k\left(t\right)$ are the rational functions of t, which are expressed through the coefficients $u_{pq}\left(t\right),$ which in turn are expressed through $b_1\left(t\right)$ and $b_2\left(t\right)$ according to (74). This expansion is valid for all values of $t,$ except maybe the neighborhood of the roots of (75). In particular,

$$\begin{array}{c}{c}_1\left(t\right)=-\frac{u_{01}}{u_{10}}=-\frac{G}{H}=\left(3\left(5000596138840425\sqrt{3}-6061042284824999\right)\left(7108208938240\sqrt{3}t+15394617958400t^2+\right.\right.\hfill \\ 1623856300668\sqrt{3}-1613861867520t+6311014555365)\left(7108208938240\sqrt{3}t+7197224345600t^2\right.\\ \left.+541285433556\sqrt{3}+7246825313280t+1592795378919)\right)/\left(42266280808032016 \times \right.\\ \times \left(69374221130894188740608000\sqrt{3}{t}^2+303911073479526952468480000t^4-9161742964858347934924800\sqrt{3}{t}^2\right.\\ -168426828577166939652096000t^3+21223516032095931917445120\sqrt{3}t\\ +80511690505906612625817600t^2+2536997119105720789608561\sqrt{3}\\ \hfill \left.+21428792681795614620161280t+5952963999951162251632530)\right)\end{array}$$

where the denominator has no real roots. According to (76) approximately.

$$\epsilon \approx c_1\left(t\right)M_3.$$

By the sequence of substitutions (64), (65), (67), (68) and (73), we return to the original coordinates, which at small $\left|M_3\right|$ by $\Omega$ are approximated with

$$\begin{array}{c}\hfill D_1=\left(b_1\left(t\right)+c_1\left(t\right)M_3\right)M_3,D_2=\left(b_2\left(t\right)+c_1\left(t\right)M_3\right)M_3,{D_3=M}_3;\end{array}$$

$$\begin{array}{cc}\hfill C_1=& -\frac{\left(2\sqrt{3}-3\right)}{6}·\left(-2c_1\left(t\right)\left(2\sqrt{3}+3\right)M_3^2+\left(b_1\left(t\right)-4\sqrt{3}{b}_2\left(t\right)-7b_2\left(t\right)\right)M_3+4\sqrt{3}+6\right),\hfill \\ \hfill C_2=& -\frac{\sqrt{3}}{6}·\left(b_1\left(t\right)M_3-b_2\left(t\right)M_3+2\sqrt{3}\right),\hfill \\ \hfill C_3=& M_3;\hfill \end{array}$$

According to (65)

$$\begin{array}{c}\hfill B_1=C_1,B_2=C_2,B_3=\frac{1}{C_3}=\frac{1}{M_3}.\end{array}$$

According to (64)

$$A_1=B_1{B}_3=b_1\left(t\right)\left(\frac{1}{2}-\frac{\sqrt{3}}{3}\right)+b_2\left(t\right)\left(\frac{1}{2}+\frac{\sqrt{3}}{3}\right)+c_1\left(t\right)M_3-\frac{1}{M_3},$$

(77)

$$A_2=B_2{B}_3=-\frac{\sqrt{3}}{6}·b_1\left(t\right)+\frac{\sqrt{3}}{6}·b_2\left(t\right)-\frac{1}{M_3},$$

(78)

$$A_3=\frac{1}{M_3}.$$

When $M_3=-0.1$, the curve (77) and (78) is shown in Figure 30.

At $M_3=0.1$, the curve (77) and (78) is shown in Figure 31.

5.3. Second Multiplier in (66)

Polynomial

$$f_2=C_1^2-4C_1{C}_2+C_2^2+4C_1+4C_2-8$$

does not factorize in the field of rational numbers, but does factorize in the extension of that field with $\alpha =\sqrt{3}$. It is the product of two linear forms $f_2=D_1·D_2$, which we treat as coordinate substitutions

$$D_1=C_1-C_2\left(\sqrt{3}+2\right)+\left(2\sqrt{3}+2\right),D_2=C_1+C_2\left(\sqrt{3}-2\right)-\left(2\sqrt{3}-2\right),D_3=C_3.$$

Its inverse substitution is

$$\begin{array}{cc}\hfill C_1=& \frac{(-2+\sqrt{3})\sqrt{3}}{6}{D}_1+\frac{(2+\sqrt{3})\sqrt{3}}{6}{D}_2+2,\hfill \\ \hfill C_2=& -\frac{\sqrt{3}}{6}{D}_1+\frac{\sqrt{3}}{6}{D}_2+2,\hfill \\ \hfill C_3=& D_3.\hfill \end{array}$$

(79)

We substitute it into the polynomial $J\left(\mathbf{C}\right)$ and get the polynomial $S_2\left(\mathbf{D}\right)=J\left(\mathbf{C}\right).$ For the polynomial $S_2\left(\mathbf{D}\right)$, we calculate Newton’s polyhedron ${\Gamma }_{10}$.

Its graph is shown in Figure 32. It has 11 two-dimensional faces with external normals

$$\begin{array}{c}{N}_{2049}=\left(-2,-2,-1\right),N_{6293}=\left(-2,0,-1\right),N_{52673}=\left(0,-2,-1\right),\\ N_{117443}=\left(1,1,1\right),N_{122283}=\left(-2,1,0\right),N_{122987}=\left(0,1,0\right),N_{123867}=\left(-1,0,0\right),\\ N_{132855}=\left(0,-1,0\right),N_{150911}=\left(1,0,0\right),N_{158011}=\left(0,0,-1\right),N_{164049}=\left(1,-2,0\right).\end{array}$$

Since $D_1,$$D_2$ and $D_3\to 0$, we select the only normal that has all coordinates negative. This is the normal $N_{2049}=\left(-2,-2,-1\right)$. It corresponds to the truncated polynomial

$$ftr2049=26244{(3{\alpha D}_3^2+3D_3^2-8D_1)}^2{(3{\alpha D}_3^2-3D_3^2+8D_2)}^2.$$

Assume

$$x_1=-8D_1+(3\alpha +3{)D}_3^2,x_2=8D_2+(3\alpha -3{)D}_3^2,x_3=D_3.$$

The inverse transformation is

$$D_1=-\frac{1}{8}{x}_1+\frac{\left(3+3\alpha \right)}{8}{x}_3^2,D_2=\frac{1}{8}{x}_2+\frac{(3-3\alpha )}{8}{x}_3^2,D_3=x_3.$$

(80)

We substitute it into the polynomial $S_2\left(\mathbf{D}\right)=J\left(\mathbf{C}\right)$ and get the polynomial $S_3\left(\mathbf{x}\right)=J\left(\mathbf{C}\right)$. For the polynomial $S_3\left(\mathbf{x}\right)$, we compute Newton’s polyhedron ${\Gamma }_{11}$.

Its graph is shown in Figure 33. It has 9 two-dimensional faces with external normals

$$\begin{array}{c}{N}_{1919}=\left(-2,0,-1\right),N_{2049}=\left(-2,-2,-1\right),N_{4559}=\left(0,-2,-1\right),\\ N_{12467}=\left(2,2,1\right),N_{14571}=\left(-1,0,0\right),N_{15095}=\left(0,1,0\right),\\ N_{15335}=\left(1,0,0\right),N_{18043}=\left(0,0,-1\right),N_{19131}=\left(0,-1,0\right).\end{array}$$

Since $x_1$, $x_2$, and $x_3\to 0$, we select the only normal that has all coordinates negative. This is $N_{2049}=\left(-2,-2,-1\right).$ According to results of our program, it corresponds to the truncated polynomial

$$\begin{array}{c}\hfill ftr2049=-13695130288521216{\alpha }^2{x}_3^8-126806761930752{\alpha }^2{x}_1^2{x}_3^4-507227047723008{\alpha }^2{x}_1{x}_2{x}_3^4\\ \hfill -126806761930752{\alpha }^2{x}_2^2{x}_3^4+41085390865563648x_3^8+380420285792256x_1^2{x}_3^4\\ \hfill +1521681143169024x_3^4{x}_1{x}_2+380420285792256x_2^2{x}_3^4+7044820107264x_1^2{x}_2^2.\end{array}$$

Doing the power transformation

$$x_1=y_1{y}_3^2,x_2=y_2{y}_3^2,x_3=y_3$$

(81)

and factorize we get

$$ftr2049=-7044820107264\left(18{\alpha }^2{y}_1^2+72{\alpha }^2{y}_1{y}_2+18{\alpha }^2{y}_2^2-y_1^2{y}_2^2+1944{\alpha }^2-54y_1^2-216y_1{y}_2-54y_2^2-5832\right)y_3^8.$$

If we substitute $\alpha =\sqrt{3}$, into the polynomial in parentheses, it is equal to $-y_1^2{y}_2^2$. Therefore, the power transformation (81) is substituted into the large polynomial $S_3\left(\mathbf{x}\right)$ and divided by $(-7+4\sqrt{3})y_3^8$ we get the polynomial $S_4\left(\mathbf{y}\right)$. We compute its Newton polyhedron ${\Gamma }_{12}$.

Its graph is shown in Figure 34. It has 11 two-dimensional faces with external normals

$$N_{4367}=\left(-1,0,0\right),N_{5643}=\left(-1,0,-1\right),N_{6099}=\left(-1,-1,-1\right),$$

$$N_{13629}=\left(0,-1,-1\right),N_{57365}=\left(0,-1,0\right),N_{111755}=\left(2,2,-1\right),N_{122609}=\left(0,1,0\right),$$

$$N_{124119}=\left(0,2,-1\right),N_{150923}=\left(1,0,0\right),N_{162019}=\left(0,0,1\right),N_{164049}=\left(2,0,-1\right).$$

Since $y_1$, $y_2$, and $y_3\to 0$, we select the only normal that has all coordinates negative. This is $N_{6099}=\left(-1,-1,-1\right)$. The corresponding shortening

$$ftr6099=-28179280429056y_1^2{y}_2^2\alpha +63403380965376y_1^2{y}_2{y}_3\alpha -31701690482688y_1^2{y}_3^2\alpha$$

$$+232479063539712y_1{y}_2^2{y}_3\alpha -507227047723008y_1{y}_2{y}_3^2\alpha +285315214344192y_1{y}_3^3\alpha$$

$$-475525357240320y_2^2{y}_3^2\alpha +1046155785928704y_2{y}_3^3\alpha -570630428688384y_3^4\alpha$$

$$-49313740750848y_1^2{y}_2^2+105672301608960y_1^2{y}_2{y}_3-63403380965376y_1^2{y}_3^2$$

$$+401554746114048y_1{y}_2^2{y}_3-887647333515264y_1{y}_2{y}_3^2+475525357240320y_1{y}_3^3$$

$$-824243952549888y_2^2{y}_3^2+1806996357513216y_2{y}_3^3-998603250204672y_3^4$$

Doing the power transformation.

$$y_1=z_1{z}_3,y_2=z_2{z}_3,y_3=z_3.$$

(82)

and we get

$$\begin{array}{c}\hfill ftr6099pow=-1761205026816z_3^4(16\alpha z_1^2{z}_2^2-36\alpha z_1^2{z}_2-132\alpha z_1{z}_2^2+28z_1^2{z}_2^2+18\alpha z_1^2\\ \hfill +288\alpha z_1{z}_2+270\alpha z_2^2-60z_1^2{z}_2-288z_1{z}_2^2-162\alpha z_1-594\alpha z_2+36z_1^2+504z_1{z}_2+468z_2^2+324\alpha \\ \hfill -270z_1-1026z_2+567)=-1761205026816z_3^4·F_{20}(z_1,z_2),\end{array}$$

where addition “$pow$” indicate that it is after power transformation.

If we substitute $\alpha =\sqrt{3}$ inside the brackets, we get the factorization

$$F_{20}\left(z_1,z_2\right)=\frac{(4\sqrt{3}+7){(3\sqrt{3}-2z_1+3)}^2{(-2z_2+3\sqrt{3}-3)}^2}{4}.$$

The power transformation (82) we do in the polynomial $S_4\left(\mathbf{y}\right)$, we get the polynomial $z_3^4{S}_5\left(\mathbf{z}\right)=S_4\left(\mathbf{y}\right)$.

In $S_5\left(\mathbf{z}\right)$, we substitute, introducing new variables $L_i$,

$$z_1=L_1+\frac{3(\sqrt{3}+1)}{2},z_2=L_2+\frac{3(\sqrt{3}-1)}{2},z_3=L_3.$$

(83)

We get the polynomial $S_6\left(\mathbf{L}\right)=S_5\left(\mathbf{z}\right)$ and for it we calculate Newton’s polyhedron ${\Gamma }_{13}$.

Its graph is shown in Figure 35. The computer computed the polyhedron ${\Gamma }_{13}$ in 87 h and 23 min. It has 9 two-dimensional faces with exterior normals

$$\begin{array}{c}{N}_{2165}=\left(0,-1,0\right),N_{5047}=\left(0,0,1\right),N_{5649}=\left(3,1,-1\right),\\ N_{6083}=\left(1,0,0\right),N_{11169}=\left(1,3,-1\right),N_{13283}=\left(-1,0,0\right),\\ N_{14775}=\left(-1,-1,-1\right),N_{17273}=\left(3,3,-1\right),N_{17669}=\left(0,1,0\right).\end{array}$$

Since $L_1$, $L_2$, and $L_3\to 0$, we select the only normal that has all coordinates negative. This is $N_{14775}=\left(-1,-1,-1\right)$. The corresponding truncation

$$\begin{array}{cc}ftr14775=\hfill & -440301256704(\alpha +1)(384\alpha L_1^2{L}_2^2+672\alpha L_1^2{L}_2{L}_3-315\alpha L_1^2{L}_3^2+2688\alpha L_1{L}_2^2{L}_3\hfill \\ & +5778\alpha L_1{L}_2{L}_3^2+756\alpha L_1{L}_3^3-2880\alpha L_2^3{L}_3-3465\alpha L_2^2{L}_3^2+3024\alpha L_2{L}_3^3+972\alpha L_3^4-192L_1^3{L}_3\hfill \\ & \hspace{1em}\hspace{1em}+640L_1^2{L}_2^2+1344L_1^2{L}_2{L}_3-315L_1^2{L}_3^2+4704L_1{L}_2^2{L}_3+9630L_1{L}_2{L}_3^2+1512L_1{L}_3^3\hfill \\ & \hfill \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}-4992L_2^3{L}_3-5985L_2^2{L}_3^2+5292L_2{L}_3^3+1620L_3^4)\end{array}$$

We do a power transformation.

$$L_1=M_1{M}_3,L_2=M_2{M}_3,L_3=M_3.$$

(84)

Then

$$\begin{array}{c}ftr14775pow=-440301256704(\alpha +1)(384\alpha M_1^2{M}_2^2+672\alpha M_1^2{M}_2+2688\alpha M_1{M}_2^2-2880\alpha M_2^3\hfill \\ \hspace{1em}+640M_1^2{M}_2^2-315\alpha M_1^2+5778\alpha M_1{M}_2-3465\alpha M_2^2-192M_1^3+1344M_1^2{M}_2+4704M_1{M}_2^2\\ \hspace{1em}-4992M_2^3+756\alpha M_1+3024\alpha M_2-315M_1^2+9630M_1{M}_2-5985M_2^2+972\alpha +1512M_1+5292M_2+1620)M_3^4\\ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}=-440301256704\left(\alpha +1\right)M_3^4·F_{40}.\hfill \end{array}$$

The curve $F_{40}=0$ has genus 0, and parameterization

$$\begin{array}{cc}\left[M_1,M_2\right]=[b_1\left(t\right),\hfill & b_2\left(t\right)]=\hfill \\ & =\left\{(5-3\sqrt{3})(1344\sqrt{3}{t}^2+t^3+2352t^2-345047040\sqrt{3}-597639168)/\left(256t^2\right),\right.\hfill \\ & \left.(265-153\sqrt{3})(-t^3+3612672\sqrt{3}t+1380188160\sqrt{3}+6257664t+2390556672)/\left(49152t\right)\right\}\hfill \end{array}$$

(85)

and is shown in Figure 36.

Figure 36 shows the limiting values of $t=-\infty ,-0,+0,+\infty ,$ when the branch goes to infinity. The approximate values of the zeros of the numerators in (85) are. $t_1\approx 481,241,t_2\approx -4623,972,t_3\approx -537,144$ for $b_1\left(t\right)$ and $t_4\approx 3715,095,t_5\approx -3328,446,t_6\approx -386,649$ for $b_2\left(t\right)$.

We do a power transformation (84) to $S_6\left(\mathbf{L}\right),$ and we get

$$-M_3^{4}E\left(\mathbf{M}\right)=S_6\left(\mathbf{L}\right)=-M_3^4\sum _{k=0}{E}_k\left(M_1,M_2\right)M_3^k,$$

where

$$\begin{array}{c}\hfill E_0=440301256704(1+\alpha )(384\alpha M_1^2{M}_2^2+672\alpha M_1^2{M}_2+2688\alpha M_1{M}_2^2-2880\alpha M_2^3+640M_1^2{M}_2^2\\ \hfill -315\alpha M_1^2+5778\alpha M_1{M}_2-3456\alpha M_2^2-192M_1^3+1344M_1^2{M}_2+4704M_1{M}_2^2-4992M_2^3+756\alpha M_1\\ \hfill +3024\alpha M_2-315M_1^2+9630M_1{M}_2-5985M_2^2+972\alpha +1512M_1+5292M_2+1620),\\ \hfill E_1=-2641807540224(1+\alpha )(112\alpha M_1^2{M}_2+448\alpha M_1{M}_2^2-480\alpha M_2^3-105\alpha M_1^2+1926\alpha M_1{M}_2\\ \hfill -1155\alpha M_2^2-32M_1^3+224M_1^2{M}_2+784M_1{M}_2^2-832M_2^3+378\alpha M_1+1512\alpha M_2-105M_1^2+3210M_1{M}_2\\ \hfill -1995M_2^2+648\alpha +756M_1+2646M_2+1080),\\ \hfill E_2=-55037657088(1+\alpha )(7680\alpha M_1^2{M}_2^2+7728\alpha M_1^2{M}_2+30912\alpha M_1{M}_2^2-21600\alpha M_2^3+12800M_1^2{M}_2^2\\ \hfill +4599\alpha M_1^2-32778\alpha M_1{M}_2+50589\alpha M_2^2-1440M_1^3+15456M_1^2{M}_2+54096M_1{M}_2^2-37440M_2^3\\ \hfill -17766\alpha M_1-71064\alpha M_2+4599M_1^2-54630M_1{M}_2+87381M_2^2-52002\alpha -35532M_1-124362M_2-86670).\end{array}$$

(86)

Into the polynomial $E\left(\mathbf{M}\right)$ we substitute

$$M_1=b_1\left(t\right),M_2=b_2\left(t\right)+\epsilon$$

(87)

according to (85). Then the polynomial $E\left(\mathbf{M}\right)$ becomes a polynomial

$$u\left(\epsilon ,M_3\right)=\sum _{p,q\ge 0}{u}_{pq}\left(t\right){\epsilon }^p{M}_3^q,$$

whereby

$$u_{pq}=\frac{1}{p!}·\frac{{\partial }^p{E}_q}{\partial M_2^p}\left(b_1\left(t\right),b_2\left(t\right)\right),$$

(88)

where $M_1=b_1\left(t\right),$ $M_2=b_2\left(t\right)$ according to (85). In particular, from (86)–(88) we obtain

$$u_{00}=E(b_1\left(t\right),b_2\left(t\right))\equiv 0,$$

$$u_{10}\left(M_1,M_2\right)=\frac{\partial E_0}{\partial M_2}=$$

$$\begin{array}{c}\hfill u_{10}\left(M_1,M_2\right)=\frac{\partial E_0}{\partial M_2}=1761205026816\left(7+4\alpha \right)\left(168\alpha M_1^2+336\alpha M_1{M}_2-432\alpha M_2^2+128M_1^2{M}_2\right.\\ \hfill \left.-315\alpha M_2-168M_1^2+336M_1{M}_2-720M_2^2+189\alpha +963M_1-630M_2+189)\right)\\ \hfill =\frac{69984(571\alpha -989){\left(384\alpha t-t^2-516096\alpha +672t-893952\right)}^3{\left(t+672+384\alpha \right)}^3}{t^5},\end{array}$$

$$\begin{array}{c}\hfill u_{01}=E_1\left(b_1\left(t\right),b_2\left(t\right)\right)=\left(\right(729\left(151316\alpha -262087\right)\left(2304\alpha t+t^2+516096\alpha +4032t+893952\right)\\ \hfill {(384\alpha t-t^2-516096\alpha +672t-893952)}^3{(t+672+384\alpha )}^4)/\left(256t^6\right)),\end{array}$$

$$\begin{array}{c}\hfill u_{11}\left(M_1,M_2\right)=\frac{\partial E_1}{\partial M_2}=-5283615080448(5+3\alpha )(112\alpha M_1{M}_2-192\alpha M_2^2+321\alpha M_1-315\alpha M_2\\ \hfill +56M_1^2+224M_1{M}_2-336M_2^2+189\alpha +321M_1-525M_2+378)=\frac{1}{t^4}(209952(571\alpha -989)·\\ \hfill ·(-3333609787974746112-18912215826432t^2+14958900893712384t+t^6+1792t^5+2681856t^4\\ \hfill +1024\alpha t^5+1548288\alpha t^4-4600627200\alpha t^3-1924660508460318720\alpha -10918972882944\alpha t^2\\ \hfill +8636525457702912\alpha t-7968522240t^3\left){\left(t+672+384\alpha \right)}^2\right),\end{array}$$

$$\begin{array}{c}\hfill u_{20}\left(M_1,M_2\right)=\frac{1}{2}\left(\frac{{\partial }^2{E}_0}{\partial M_2^2}\right)=880602513408\left(7+4\alpha \right)\left(336\alpha M_1-864\alpha M_2+128M_1^2-315\alpha \right.\\ \hfill \left.+336M_1-1440M_2-630\right)=(6879707136(\alpha -2)(144\alpha t-t^2-69120\alpha +240t-119808)\times \\ \hfill \times \left(528\alpha t+t^2+963072\alpha +912t+1668096\right)(t+240+144\alpha )(-t+912+528\alpha )/t^4).\end{array}$$

According to [1] (Theorem 1), the solution to the equation $u\left(\epsilon ,M_3\right)=0$ has the form, $\epsilon ={\sum }_{k=1}^{\infty }{c}_k\left(t\right)M_3^k$, where

$$\begin{array}{c}\hfill c_1\left(t\right)=-\frac{u_{01}\left(t\right)}{u_{10}\left(t\right)}=\frac{\left(-265+153\alpha \right)\left(2304\alpha t+t^2+516096\alpha +4032t+893952\right)\left(t+672+384\alpha \right)}{49152t}.\end{array}$$

The denominators in $c_1\left(t\right)$ and in (85) have root $t=0$.

Now let’s go back and approximately from (84) obtain

$$L_1=b_1\left(t\right)M_3,L_2=\left(b_2\left(t\right)+c_1\left(t\right)M_3\right)M_3=b_2\left(t\right)M_3=b_2\left(t\right)M_3+c_1\left(t\right)M_3^2,L_3=M_3.$$

Substitute that into (83).

$$z_1=b_1\left(t\right)M_3+\frac{3(\sqrt{3}+1)}{2},z_2=b_2\left(t\right)M_3+c_1\left(t\right)M_3^2+\frac{3(\sqrt{3}-1)}{2},z_3=M_3.$$

(89)

We substitute the values of (89) into (82) and obtain

$$\begin{array}{cc}\hfill y_1=& b_1\left(t\right)M_3^2+\frac{3(\sqrt{3}+1)}{2}{M}_3,\hfill \\ \hfill y_2=& c_1\left(t\right)M_3^3+b_2\left(t\right)M_3^2+\frac{3(\sqrt{3}-1)}{2}{M}_3,\hfill \\ \hfill y_3=& M_3.\hfill \end{array}$$

(90)

We substitute the values of (90) into (81) and obtain

$$\begin{array}{cc}\hfill x_1=& y_1{y}_3^2=b_1\left(t\right)M_3^4+\frac{3(\sqrt{3}+1)}{2}{M}_3^3,\hfill \\ \hfill x_2=& y_2{y}_3^2=c_1\left(t\right)M_3^5+b_2\left(t\right)M_3^4+\frac{3(\sqrt{3}-1)}{2}{M}_3^3,\hfill \\ \hfill x_3=& y_3=M_3.\hfill \end{array}$$

(91)

We substitute the values of (91) into (80) and obtain

$$\begin{array}{cc}\hfill D_1=& -\frac{1}{8}{x}_1+\frac{\left(3+3\alpha \right)}{8}{x}_3^2=-\frac{1}{8}{b}_1\left(t\right)M_3^4-\frac{3\left(\sqrt{3}+1\right)}{16}{M}_3^3+\frac{\left(3+3\sqrt{3}\right)}{8}{M}_3^2,\hfill \\ \hfill D_2=& \frac{1}{8}{x}_2+\frac{(3-3\alpha )}{8}{x}_3^2=\frac{1}{8}{c}_1\left(t\right)M_3^5+b_2\left(t\right)M_3^4+\frac{3(\sqrt{3}-1)}{16}{M}_3^3+\frac{(3-3\sqrt{3})}{8}{M}_3^2,\hfill \\ \hfill D_3=& M_3.\hfill \end{array}$$

(92)

We substitute the values of (92) into (79) and obtain

$$\begin{array}{cc}\hfill C_1& =\frac{\left(-2+\sqrt{3}\right)\sqrt{3}}{6}{D}_1+\frac{\left(2+\sqrt{3}\right)\sqrt{3}}{6}{D}_2+2=\frac{\left(2\sqrt{3}+3\right)}{48}\left(c_1\left(t\right)M_3^5-4\sqrt{3}{b}_1\left(t\right)M_3^4\right.\hfill \\ \hfill & \left.+7b_1\left(t\right)M_3^4+8b_2\left(t\right)M_3^4-9M_3^3+6\sqrt{3}{M}_3^3+18M_3^2-12\sqrt{3}{M}_3^2-96+64\sqrt{3}\right),\hfill \\ \hfill C_2& =-\frac{\sqrt{3}}{6}{D}_1+\frac{\sqrt{3}}{6}{D}_2+2=\frac{\sqrt{3}}{48}\left(c_1\left(t\right)M_3^5+b_1\left(t\right)M_3^4+8b_2\left(t\right)M_3^4+3\sqrt{3}{M}_3^3-6\sqrt{3}{M}_3^2+32\sqrt{3}\right),\hfill \\ \hfill C_3& =M_3.\hfill \end{array}$$

(93)

We substitute the values of (93) into (65) and obtain

$$\begin{array}{cc}\hfill B_1=& C_1=\frac{\left(2\sqrt{3}+3\right)}{48}(c_1\left(t\right)M_3^5-4\sqrt{3}{b}_1\left(t\right)M_3^4+7b_1\left(t\right)M_3^4+8b_2\left(t\right)M_3^4-9M_3^3\hfill \\ \hfill & +6\sqrt{3}{M}_3^3+18M_3^2-12\sqrt{3}{M}_3^2-96+64\sqrt{3}),\hfill \\ \hfill B_2=& C_2=\frac{\sqrt{3}}{48}\left(c_1\left(t\right)M_3^5+b_1\left(t\right)M_3^4+8b_2\left(t\right)M_3^4+3\sqrt{3}{M}_3^3-6\sqrt{3}{M}_3^2+32\sqrt{3}\right),\hfill \\ \hfill B_3=& C_3^{-1}=\frac{1}{M_3}.\hfill \end{array}$$

(94)

Finally, we substitute the values of (94) into (64) and obtain values $A_1,A_2,A_3$:

$$\begin{array}{cc}\hfill A_1=& B_1{B}_3=\frac{\left(2\sqrt{3}+3\right)}{48}(c_1\left(t\right)M_3^4-4\sqrt{3}{b}_1\left(t\right)M_3^3+7b_1\left(t\right)M_3^3+8b_2\left(t\right)M_3^3\hfill \end{array}$$

$$\begin{array}{cc}\hfill & -9M_3^2+6\sqrt{3}{M}_3^2+18M_3-12\sqrt{3}{M}_3-\frac{96}{M_3}+\frac{64\sqrt{3}}{M_3}),\hfill \end{array}$$

(95)

$$\begin{array}{cc}\hfill A_2=& B_2{B}_3=\frac{\sqrt{3}}{48}(c_1\left(t\right)M_3^4+b_1\left(t\right)M_3^3+8b_2\left(t\right)M_3^3+3\sqrt{3}{M}_3^2-6\sqrt{3}{M}_3+\frac{32\sqrt{3}}{M_3},\hfill \\ \hfill A_3=& B_3=\frac{1}{M_3}.\hfill \end{array}$$

(96)

We need them to drive figures. Figure 37 and Figure 38 show the curves of (95) and (96) at $M_3=0.1$ and $M_3=-0.1$, respectively.

5.4. Third Multiplier in (66)

It is a linear multiplier defined with

$$f_3=C_1+C_2+2.$$

The substitution is

$$D_1=C_1+C_2+2,D_2=C_2,D_3=C_3.$$

Its inverse substitution is

$$C_1=D_1-D_2-2,C_2=D_2,C_3=D_3.$$

(97)

Let’s consider all of this as a coordinate substitution in the polynomial $J\left(\mathbf{C}\right)$. We substitute it into the polynomial $J\left(\mathbf{C}\right)$ and get the polynomial $S_7\left(\mathbf{D}\right)=J\left(\mathbf{C}\right)$. For the polynomial $S_7\left(\mathbf{D}\right)$, we compute Newton’s polyhedron ${\Gamma }_{14}.$

Its graph is shown in Figure 39. It has 7 two-dimensional faces with external normals

$$N_{701}=\left(-2,0,-1\right),N_{1403}=\left(1,1,1\right),N_{1537}=\left(0,0,-1\right),$$

$$N_{1679}=\left(0,-1,0\right),N_{1935}=\left(-2,1,0\right),N_{2009}=\left(0,1,0\right),N_{2085}=\left(-1,0,0\right).$$

Since $D_1\to 0$, $D_2\to const$, and $D_3\to 0$, we select the only normal whose first and third coordinates are negative. This is the normal $N_{701}=\left(-2,0,-1\right)$. It corresponds to a truncated polynomial

$$\begin{array}{c}ftr701=186624{\left(D_2+1\right)}^2{\left(D_2^2+2D_2-2\right)}^2\times \hfill \\ \hfill \times \left(-48D_1{D}_2^2{D}_3^2+64D_1^2{D}_2^2-96D_1{D}_2{D}_3^2+27D_3^4+128D_1^2{D}_2-48D_1{D}_3^2+64D_1^2\right).\end{array}$$

Do a power transformation.

$$D_1=x_1{x}_3^2,D_2=x_2,D_3=x_3.$$

(98)

Then ftr701 after the power transformation (98) is

$$\begin{array}{c}ftr701pow=186624{\left(x_2+1\right)}^2{\left(x_2^2+2x_2-2\right)}^2{x}_3^4\left(64x_1^2{x}_2^2+128x_1^2{x}_2-48x_1{x}_2^2+64x_1^2-96x_1{x}_2-48x_1+27\right)\hfill \\ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}=186624{\left(x_2+1\right)}^2{\left(x_2^2+2x_2-2\right)}^2{x}_3^4·F_{30}\left(x_1,x_2\right),\hfill \end{array}$$

where $F_{30}\left(x_1,x_2\right)=64x_1^2{x}_2^2+128x_1^2{x}_2-48x_1{x}_2^2+64x_1^2-96x_1{x}_2-48x_1+27$. The curve $F_{30}=0$ has genus 0, and parameterization

$$\{x_1,x_2\}=\left\{b_1\left(t\right),b_2\left(t\right)\right\}=\left\{\frac{81{\left(29t+3\right)}^2}{4\left(23283t^2+5466t+499\right)},-\frac{27459t^2+8682t+787}{48\left(3t+2\right)\left(29t+3\right)}\right\}$$

(99)

and is shown in the Figure 40.

The curve $F_{30}=0$ goes to infinity $x_2=\pm \infty$ at

$$t_1=-\frac{3}{29}\approx -0.1034482759,t_2=-\frac{2}{3}\approx -0.6666666667.$$

(100)

We do a power transformation of (98) to the polynomial $S_7\left(\mathbf{D}\right)$ and get the polynomial

$$x_3^{4}W\left(\mathbf{x}\right)=S_7\left(\mathbf{D}\right)=x_3^4\sum _{k=0}{W}_k(x_1,x_2)x_3^k.$$

where

$$\begin{array}{cc}\hfill W_0=& 186624{\left(x_2+1\right)}^2{\left(x_2^2+2x_2-2\right)}^2\left(64x_1^2{x}_2^2+128x_1^2{x}_2-48x_1{x}_2^2+64x_1^2-96x_1{x}_2-48x_1+27\right),\hfill \\ \hfill W_1=& 559872{\left(x_2+1\right)}^2{\left(x_2^2+2x_2-2\right)}^2\left(8x_1{x}_2^2+16x_1{x}_2+8x_1-9\right),\hfill \\ \hfill W_2=& -15552(x_2^2+2x_2-2)(512x_1^3{x}_2^6+6144x_1^3{x}_2^5+19968x_1^3{x}_2^4-2304x_1^2{x}_2^5+24064x_1^3{x}_2^3\hfill \\ \hfill & -10080x_1^2{x}_2^4+6144x_1^3{x}_2^2-13824x_1^2{x}_2^3-7680x_1^3{x}_2-8352x_1^2{x}_2^2+972x_1{x}_2^3-81x_2^4-4096x_1^3\hfill \\ \hfill & -4032x_1^2{x}_2+4860x_1{x}_2^2-324x_2^3-1728x_1^2+5832x_1{x}_2-243x_2^2+1944x_1+162x_2-1296).\hfill \end{array}$$

(101)

Into the polynomial $W\left(\mathbf{x}\right)$ we substitute

$$x_1=b_1\left(t\right),x_2=b_2\left(t\right)+\epsilon$$

(102)

according to (99). Then the polynomial $W\left(\mathbf{x}\right)$ becomes a polynomial

$$u\left(\epsilon ,x_3\right)=\sum _{p,q\ge 0}{u}_{pq}\left(t\right){\epsilon }^p{x}_3^q,$$

whereby

$$u_{pq}=\frac{1}{p!}·\frac{{\partial }^p{W}_q}{\partial x_2^p}\left(x_1,x_2\right),$$

(103)

where $x_1=b_1\left(t\right),$ $x_2=b_2\left(t\right)$ according to (99). In particular, from (101)–(103) we obtain

$$\begin{array}{c}\hfill u_{00}=W(b_1\left(t\right),b_2\left(t\right))\equiv 0,\\ \hfill u_{10}\left(x_1,x_2\right)=\frac{\partial W_0}{\partial x_2}=373248\left(x_2+1\right)\left(x_2^2+2x_2-2\right)\\ \hfill \times \left(256x_1^2{x}_2^4+1024x_1^2{x}_2^3-192x_1{x}_2^4+1152x_1^2{x}_2^2-\right.\\ \hfill \left.-768x_1{x}_2^3+256x_1^2{x}_2-864x_1{x}_2^2-128x_1^2-192x_1{x}_2+81x_2^2+96x_1+162x_2)\right)\\ \hfill =\frac{81(23283t^2+5466t+499){\left(22131t^2+3930t-13\right)}^4}{2048{\left(29t+3\right)}^5{\left(3t+2\right)}^5},\end{array}$$

$$\begin{array}{c}\hfill u_{01}=W_1\left(b_1\left(t\right),b_2\left(t\right)\right)=\frac{{27\left(23283t^2+5466t+499\right)}^2{\left(22131t^2+3930t-13\right)}^5}{8388608{\left(29t+3\right)}^6{\left(3t+2\right)}^8}.\end{array}$$

According to Theorem 1 of [1], the solution of equation $u\left(\epsilon ,x_3\right)=0$ is

$\epsilon ={\sum }_{k=1}^{\infty }{c}_k\left(t\right)x_3^k$, where

$$c_1\left(t\right)=-\frac{u_{01}\left(t\right)}{u_{10}\left(t\right)}=-\frac{\left(23283t^2+5466t+499\right)\left(22131t^2+3930t-13\right)}{12288\left(29t+3\right){\left(3t+2\right)}^3}.$$

The denominator in $c_1\left(t\right)$ has roots (100).

Now let’s go back and approximate from (102) obtain

$$x_1=b_1\left(t\right),x_2=b_2\left(t\right)+c_1\left(t\right)x_3.$$

Substitute that into (98) and we get

$$D_1=b_1\left(t\right)x_3^2,D_2=b_2\left(t\right)+c_1\left(t\right)x_3,D_3=x_3.$$

(104)

We substitute the expression (104) into (97) and obtain

$$\begin{array}{cc}\hfill C_1=& D_1-D_2-2=b_1\left(t\right)x_3^2-b_2\left(t\right)-c_1\left(t\right)x_3-2,\hfill \\ \hfill C_2=& b_2\left(t\right)+c_1\left(t\right)x_3,\hfill \\ \hfill C_3=& x_3.\hfill \end{array}$$

(105)

We substitute (105) into (65) and we get

$$B_1=C_1=b_1\left(t\right)x_3^2-b_2\left(t\right)-c_1\left(t\right)x_3-2,B_2=C_2=b_2\left(t\right)+c_1\left(t\right)x_3,B_3=C_3^{-1}=\frac{1}{x_3}.$$

Finally, according to (64)

$$\begin{array}{cc}\hfill A_1=& B_1{B}_3=\left(b_1\left(t\right)x_3^2-b_2\left(t\right)-c_1\left(t\right)x_3-2\right)/x_3=b_1\left(t\right)x_3-c_1\left(t\right)-(b_2\left(t\right)+2)/x_3,\hfill \end{array}$$

(106)

$$\begin{array}{cc}\hfill A_2=& B_2{B}_3=c_1\left(t\right)+b_2\left(t\right)/x_3,\hfill \\ \hfill A_3=& B_3=1/x_3.\hfill \end{array}$$

(107)

In Figure 41, Figure 42, Figure 43 and Figure 44, are shown the curves of (106) and (107) for values $x_3=-1$, $-0.1$, 1, $0.1$, respectively.

The another branch is symmetric to this one with respect to the line $A_1=A_2$.

5.5. Fourth Multiplier in (66)

5.5.1. Preliminary Calculations

The 4th multiplier is a linear multiplier $f_4=C_1+C_2-1$. Let’s do the substitution $D_1=C_1+C_2-1$, $D_2=C_2$, $D_3=C_3$. Its inverse substitution is

$$C_1=D_1-D_2+1,C_2=D_2,C_3=D_3.$$

(108)

We treat it all as a coordinate transformation in the polynomial $J\left(\mathbf{C}\right)$. We substitute it into the polynomial $J\left(\mathbf{C}\right)$ and get the polynomial $S_8\left(\mathbf{D}\right)=J\left(\mathbf{C}\right)$. For the polynomial $S_8\left(\mathbf{D}\right)$, we compute Newton’s polyhedron ${\Gamma }_{15}$, with graph given in Figure 45. It has 7 two-dimensional faces with external normals

$$\begin{array}{cc}\hfill N_{593}=& \left(-1,0,-1\right),N_{1295}=\left(1,1,1\right),N_{1645}=\left(0,0,-1\right),\hfill \\ \hfill N_{1941}=& \left(-1,0,0\right),N_{1989}=\left(0,3,2\right),N_{2111}=\left(0,-1,0\right),N_{2171}=\left(0,1,0\right).\hfill \end{array}$$

Since $D_1\to 0$, $D_2\to \mathrm{const}$, and $D_3\to 0$, we select the only normal that has negative first and third coordinates. This is $N_{593}=\left(-1,0,-1\right)$ and it corresponds to the truncated polynomial

$$ftr593=186624{\left(2D_2^2-2D_2-1\right)}^4\left(2D_1-3D_3\right)\left(2D_1+3D_3\right).$$

Let’s do the substitution $x_1=2D_1-3D_3$, $x_2=D_2$, $x_3=2D_1+3D_3$. Its inverse substitution is

$$D_1=\frac{1}{4}{x}_1+\frac{1}{4}{x}_3,D_2=x_2,D_3=-\frac{1}{6}{x}_1+\frac{1}{6}{x}_3,$$

(109)

and treat it all as a coordinate change in the polynomial $S_8\left(\mathbf{D}\right)$. Substitute it into the polynomial $S_8\left(\mathbf{D}\right)$ and get the polynomial $S_9\left(\mathbf{x}\right)=S_8\left(\mathbf{D}\right)$. For the polynomial $S_9\left(\mathbf{x}\right)$, we calculate Newton’s polyhedron ${\Gamma }_{17}$.

Its graph is shown in Figure 46. It has 8 two-dimensional faces with external normals

$$\begin{array}{c}\hfill N_{593}=\left(-1,0,-2\right),N_{1923}=\left(-2,0,-1\right),N_{3725}=\left(1,1,1\right),N_{4561}=\left(0,0,-1\right),\\ \hfill N_{4601}=\left(0,1,-2\right),N_{5649}=\left(-1,0,0\right),N_{6363}=\left(0,1,0\right),N_{6485}=\left(0,-1,0\right).\end{array}$$

(110)

Since $x_1\to 0$, $x_2\to const$, and $x_3\to 0$, we select two normals whose first and third coordinates are negative. These are $N_{593}=\left(-1,0,-2\right)$ and $N_{1923}=\left(-2,0,-1\right)$. We will deal with them in separate subsubsections.

5.5.2. The Normal $N_{593}=(-1,0,-2)$

According to result of our program it corresponds to a truncated polynomial

$$ftr593=764411904x_1{\left(2x_2^2-2x_2-1\right)}^3(32x_2^2{x}_3+x_1^2-32x_2{x}_3-16x_3).$$

Making a power transformation

$$x_1=y_1,x_2=y_2,x_3=y_1^2{y}_3,$$

(111)

We get a polynomial

$$\begin{array}{c}ftr593pow=764411904y_1^3{\left(2y_2^2-2y_2-1\right)}^3\left(32y_2^2{y}_3-32y_2{y}_3-16y_3+1\right)=\hfill \\ \hfill \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}=764411904y_1^3{\left(2y_2^2-2y_2-1\right)}^3·F_{50}\left(y_2,y_3\right),\end{array}$$

where $F_{50}=32y_2^2{y}_3-32y_2{y}_3-16y_3+1$. The curve $F_{50}=0$ has genus 0, parameterization

$$\left\{y_2,y_3\right\}=\left\{b_2\left(t\right),b_3\left(t\right)\right\}=\left\{t,-\frac{1}{16(2t^2-2t-1)}\right\}$$

(112)

and is shown in Figure 47.

In (112), the denominator in $b_2\left(t\right)$ has 2 real roots

$$t_1=\frac{1}{2}+\frac{\sqrt{3}}{2}\approx 1.366025404,t_2=\frac{1}{2}-\frac{\sqrt{3}}{2}\approx -0.3660254040.$$

(113)

In fact, here we can also compute the parametric expansion of the $\Omega$ manifold. To do this, we do the power transformation (110) in the polynomial $S_9\left(\mathbf{x}\right)$ and get the polynomial

$$y_1^{3}U\left(\mathbf{y}\right)=S_9\left(\mathbf{x}\right)=y_1^3\sum _{k=0}{U}_k(y_2,y_3)y_1^k.$$

In the polynomials $U_k\left(y_2,y_3\right)$ according to (112) we substitute

$$y_2=b_1\left(t\right)+\epsilon ,y_3=b_2\left(t\right)+\epsilon .$$

We obtain the polynomial $u\left(\epsilon ,y_1\right)=U\left(y_1,y_2,y_3\right)$ with coefficients depending on t through $b_1\left(t\right)$ and $b_2\left(t\right)$. In this polynomial

$$u\left(\epsilon ,y_1\right)=\sum _{k=0}^m{U}_k(b_1+\epsilon ,b_2+\epsilon )y_1^k=\sum _{p,q\ge 0}{u}_{pq}{\epsilon }^p{y}_1^q,$$

where $u_{pq}={\sum }_{p_1+p_2=p\ge 1}\frac{1}{p_1!p_2!}·\frac{{\partial }^p{U}_q}{\partial y_2^{p_1}\partial y_3^{p_2}}$ when $y_i=b_i\left(t\right)$, $i=2,3$, $p_1$, $p_2\ge 0$, $p\ge 1$. Specifically,

$$u_{00}\equiv 0,u_{10}=\frac{\partial U_0\left(y_2,y_3\right)}{\partial y_2}+\frac{\partial U_0\left(y_2,y_3\right)}{\partial y_3}$$

$$=1528823808{\left(2y_2^2-2y_2-1\right)}^2\left(32y_2^4+256y_2^3{y}_3-64y_2^3-384y_2^2{y}_3+38y_2+64y_3+5\right)$$

$$=1528823808\left(2t-3\right)\left(16t^3-8t^2-12t-3\right){\left(2t^2-2t-1\right)}^2\stackrel{\mathrm{def}}{=}H\left(b_1\left(t\right),b_2\left(t\right)\right),$$

$$u_{01}=U_1\left(b_1,b_2\right)=31850496\left(2t+5\right)\left(2t-1\right){\left(2t^2-2t-1\right)}^2\stackrel{\mathrm{def}}{=}G\left(b_1\left(t\right),b_2\left(t\right)\right).$$

The function $u_{10}\left(t\right)$ has real roots

$$\begin{array}{cc}\hfill t_1=& \frac{3}{2}=1.5,\hfill \\ \hfill t_2=& \frac{1}{2}+\frac{\sqrt{3}}{2}\approx 1.3660254042-\mathrm{multiple},\hfill \\ \hfill t_3=& \frac{1}{2}-\frac{\sqrt{3}}{2}\approx -0.3660254040,2-\mathrm{multiple},\hfill \\ \hfill t_4\approx & 1.232176060,\hfill \end{array}$$

(114)

and the function $u_{01}\left(t\right)$ has 2-multiple roots $t_2$, $t_3$ and

$$t_5=-\frac{5}{2}=-2.5,t_6=\frac{1}{2}.$$

By the Implicit Function Theorem [1] (Theorem 1), the equation $u\left(\epsilon ,y_1\right)=0$ has a solution as a power series on $y_1$

$$\epsilon =\sum _{k=1}^{\infty }{c}_k\left(t\right)y_1^k,$$

where $c_k\left(t\right)$ are rational functions of t, which are expressed through the coefficients $u_{pq}\left(t\right)$, which in turn are expressed through $b_1\left(t\right)$ and $b_2\left(t\right)$ according to (110). This expansion is valid for all values of $t,$ except maybe the neighborhood of the roots of (114). In particular,

$$c_1\left(t\right)=-\left(\frac{u_{01}}{u_{10}}\right)=-\frac{G}{H}=-\frac{(2t+5)(2t-1)}{48(2t-3)(16t^3-8t^2-12t-3)},$$

where the denominator has 2 real roots $t_1$ and $t_4$ of (114). Approximately we get

$$\epsilon \approx c_1\left(t\right)y_1.$$

Let’s return to the previous coordinates, which are approximated to be equal for small $\left|y_3\right|$ on the manifold $\Omega$

$$y_2=b_1\left(t\right)+c_1\left(t\right)y_1,y_3=b_2\left(t\right)+c_1\left(t\right)y_1.$$

$$x_1=y_1,x_2=y_2=b_1\left(t\right)+c_1\left(t\right)y_1,x_3=y_1^2{y}_3=b_2\left(t\right)y_1^2+c_1\left(t\right)y_1^3.$$

(115)

We substitute the expressions (115) into the transformation (109) and get

$$\begin{array}{cc}\hfill D_1=& \frac{1}{4}{x}_1+\frac{1}{4}{x}_3=\frac{1}{4}{y}_1+\frac{1}{4}{b}_2\left(t\right)y_1^2+\frac{1}{4}{c}_1\left(t\right)y_1^3,\hfill \\ \hfill D_2=& x_2=b_1\left(t\right)+c_1\left(t\right)y_1,\hfill \\ \hfill D_3=& -\frac{1}{6}{x}_1+\frac{1}{6}{x}_3=-\frac{1}{6}{y}_1+\frac{1}{6}{b}_2\left(t\right)y_1^2+\frac{1}{6}{c}_1\left(t\right)y_1^3.\hfill \end{array}$$

(116)

We substitute the expressions (116) into the transformation (109) and obtain variables defind in (108)

$$\begin{array}{cc}\hfill C_1=& D_1=b_1\left(t\right)+c_1\left(t\right)y_1,\hfill \\ \hfill C_2=& D_1-D_2+1=\frac{1}{4}{y}_1+\frac{1}{4}{b}_2\left(t\right)y_1^2+\frac{1}{4}{c}_1\left(t\right)y_1^3-b_1\left(t\right)-c_1\left(t\right)y_1+1,\hfill \\ \hfill C_3=& D_3=-\frac{1}{6}{y}_1+\frac{1}{6}{b}_2\left(t\right)y_1^2+\frac{1}{6}{c}_1\left(t\right)y_1^3.\hfill \end{array}$$

(117)

Substitute (117) into (65) and obtain

$$\begin{array}{cc}\hfill B_1=& C_1=b_1\left(t\right)+c_1\left(t\right)y_1,\hfill \\ \hfill B_2=& C_2=\frac{1}{4}{c}_1\left(t\right)y_1^3+\frac{1}{4}{b}_2\left(t\right)y_1^2+\frac{1}{4}{y}_1-c_1\left(t\right)y_1-b_1\left(t\right)+1,\hfill \\ \hfill B_3=& C_3^{-1}=\frac{1}{C_3}=\frac{1}{-\frac{1}{6}{y}_1+\frac{1}{6}{b}_2\left(t\right)y_1^2+\frac{1}{6}{c}_1\left(t\right)y_1^3}.\hfill \end{array}$$

(118)

Finally, we substitute the expressions (118) into the transformation (64) and obtain

$$A_1=B_1{B}_3=\frac{b_1\left(t\right)+c_1\left(t\right)y_1}{-\frac{1}{6}{y}_1+\frac{1}{6}{b}_2\left(t\right)y_1^2+\frac{1}{6}{c}_1\left(t\right)y_1^3}=-\frac{6\left(b_1\left(t\right)+c_1\left(t\right)y_1\right)}{y_1-b_2\left(t\right)y_1^2-c_1\left(t\right)y_1^3},$$

(119)

$$\begin{array}{c}\hfill A_2=B_2{B}_3=\frac{\frac{1}{4}{c}_1\left(t\right)y_1^3+\frac{1}{4}{b}_2\left(t\right)y_1^2+\frac{1}{4}{y}_1-c_1\left(t\right)y_1-b_1\left(t\right)+1}{-\frac{1}{6}{y}_1+\frac{1}{6}{b}_2\left(t\right)y_1^2+\frac{1}{6}{c}_1\left(t\right)y_1^3}\\ \hfill =-\frac{\frac{3}{2}{c}_1\left(t\right)y_1^3+\frac{3}{2}{b}_2\left(t\right)y_1^2+\frac{3}{2}{y}_1-6c_1\left(t\right)y_1-6b_1\left(t\right)+6}{y_1-b_2\left(t\right)y_1^2-c_1\left(t\right)y_1^3},\end{array}$$

(120)

$$A_3=B_3=-\frac{6}{y_1-b_2\left(t\right)y_1^2-c_1\left(t\right)y_1^3}.$$

Figure 48 and Figure 49, show the curves (119) and (120) for values $y_3=1$ and $-1$, respectively.

5.5.3. The Normal $N_{1923}=(-2,0,-1)$ from (110)

It corresponds to a truncated polynomial

$$ftr1923=254803968x_3{\left(2x_2^2-2x_2-1\right)}^3\left(4x_2^2{x}_3^2+96x_1{x}_2^2-4x_2{x}_3^2-96x_1{x}_2+x_3^2-48x_1\right).$$

By the power transformation

$$x_1=y_1{y}_3^2,x_2=y_2,x_3=y_3.$$

(121)

We have

$$\begin{array}{c}\hfill ftr1923pow=254803968y_3^3{\left(2y_2^2-2y_2-1\right)}^3\left(96{y_{1}y}_2^2-96y_1{y}_2+4y_2^2-48y_1-4y_2+1\right)=\\ \hfill \hspace{1em}\hspace{1em}\hspace{1em}=254803968y_3^3{\left(2y_2^2-2y_2-1\right)}^3·F_{60}\left(y_1,y_2\right),\end{array}$$

where $F_{60}\left(y_1,y_2\right)=96{y_{1}y}_2^2-96y_1{y}_2+4y_2^2-48y_1-4y_2+1$. The curve $F_{60}=0$ has genus 0, and parameterization

$$\left\{y_1,y_2\right\}=\left\{b_1\left(t\right),b_2\left(t\right)\right\}=\{-\frac{{(2t-1)}^2}{48(2t^2-2t-1)},t\}.$$

(122)

It is shown in Figure 50.

In (122), the denominator in $b_1\left(t\right)$ has 2 real roots $t_1,t_2$ given by (113). In fact, the parametric expansion of the manifold $\Omega$ can also be calculated here. To do this, we do a power transformation (121) in the polynomial $S_9\left(\mathbf{x}\right)$ and get the polynomial

$$-y_3^{3}P\left(\mathbf{y}\right)=S_9\left(\mathbf{x}\right)=-y_3^3\sum _{k=0}{P}_k\left(y_1,y_2\right)y_3^k.$$

Into the polynomials $P_k\left(y_1,y_2\right)$ we substitute

$$y_1=b_1\left(t\right)+\epsilon ,y_2=b_2\left(t\right)+\epsilon =t+\epsilon$$

according to (122).

We obtain the polynomial $u\left(\epsilon ,y_3\right)=P\left(y_1,y_2,y_3\right)$ with coefficients depending on t through $b_1\left(t\right)$ and $b_2\left(t\right).$ In this polynomial

$$u\left(\epsilon ,y_3\right)=\sum _{k=0}^m{P}_k(b_1+\epsilon ,b_2+\epsilon )y_3^k=\sum _{p,q\ge 0}{u}_{pq}{\epsilon }^p{y}_3^q,$$

where $u_{pq}={\sum }_{p_1+p_2=p\ge 1}\frac{1}{p_1!p_2!}·\frac{{\partial }^p{P}_q}{\partial y_1^{p_1}\partial y_2^{p_2}}$ where $y_i=b_i=b_i\left(t\right)$, $i=1,2$, $p_1,p_2\ge 0$, $p\ge 1$. Specifically

$$\begin{array}{cc}\hfill u_{00}\equiv & 0,\hfill \\ \hfill u_{10}=& \frac{\partial P_0\left(y_1,y_2\right)}{\partial y_1}+\frac{\partial P_0\left(y_1,y_2\right)}{\partial y_2}=-509607936{\left(2y_2^2-2y_2-1\right)}^2\times \hfill \\ \hfill & \times \left(768y_1{y}_2^3+96y_2^4-1152y_1{y}_2^2-160y_2^3-48y_2^2+192y_1+114y_2+23\right)=\hfill \\ \hfill =& -1528823808\left(2t-3\right)\left(16t^3-8t^2-12t-3\right){\left(2t^2-2t-1\right)}^2\stackrel{\mathrm{def}}{=}H\left(b_1\left(t\right),b_2\left(t\right)\right),\hfill \\ \hfill u_{01}=& P_1\left(b_1,b_2\right)=-31850496\left(2t+5\right)\left(2t-1\right){\left(2t^2-2t-1\right)}^2\stackrel{\mathrm{def}}{=}G\left(b_1\left(t\right),b_2\left(t\right)\right).\hfill \end{array}$$

(123)

The function $u_{10}\left(t\right)$ has real roots (114). By the Implicit Function Theorem [1] (Theorem 1), the equation $u\left(\epsilon ,y_3\right)=0$ has a solution as a power series on $y_3$

$$\epsilon =\sum _{k=1}^{\infty }{c}_k\left(t\right)y_3^k,$$

where $c_k\left(t\right)$ are rational functions of t, which are expressed through the coefficients $u_{pq}\left(t\right)$, which in turn are expressed through $b_1\left(t\right)$ and $b_2\left(t\right)$ according to (122). This expansion is valid for all values of $t,$ except maybe the neighborhood of the roots of the polynomial (123). In particular,

$$c_1\left(t\right)=-\left(\frac{u_{01}}{u_{10}}\right)=-\frac{G}{H}=-\frac{\left(2t+5\right)\left(2t-1\right)}{48\left(2t-3\right)\left(16t^3-8t^2-12t-3\right)},$$

where the denominator has 2 real roots $t_1$ and $t_4$ of (114). We get an approximation

$$\epsilon \approx c_1\left(t\right)y_3.$$

Let’s return to the previous coordinates, which, for small $\left|y_3\right|$ on $\Omega$, are approximately equal to

$$y_1=b_1\left(t\right)+c_1\left(t\right)y_3,y_2=b_2\left(t\right)+c_1\left(t\right)y_3.$$

(124)

We substitute the expressions (124) into the transformation (121) and get

$$x_1=y_1{y}_3^2=b_1\left(t\right)y_3^2+c_1\left(t\right)y_3^3,x_2=y_2=b_2\left(t\right)+c_1\left(t\right)y_3,x_3=y_3.$$

(125)

We substitute the expressions (125) into the transformation (109) and obtain

$$\begin{array}{cc}\hfill D_1=& \frac{1}{4}{x}_1+\frac{1}{4}{x}_3=\frac{1}{4}{b}_1\left(t\right)y_3^2+\frac{1}{4}{c}_1\left(t\right)y_3^3+\frac{1}{4}{y}_3,\hfill \\ \hfill D_2=& x_2=b_2\left(t\right)+c_1\left(t\right)y_3,\hfill \\ \hfill D_3=& -\frac{1}{6}{x}_1+\frac{1}{6}{x}_3=-\frac{1}{6}{b}_1\left(t\right)y_3^2-\frac{1}{6}{c}_1\left(t\right)y_3^3+\frac{1}{6}{y}_3.\hfill \end{array}$$

(126)

We substitute the expressions (126) into the transformation (108) and get the following results

$$\begin{array}{cc}\hfill C_1=& D_2=b_2\left(t\right)+c_1\left(t\right)y_3,\hfill \\ \hfill C_2=& D_1-D_2+1=\frac{1}{4}{b}_1\left(t\right)y_3^2+\frac{1}{4}{c}_1\left(t\right)y_3^3+\frac{1}{4}{y}_3-b_2\left(t\right)-c_1\left(t\right)y_3+1,\hfill \\ \hfill C_3=& D_3=-\frac{1}{6}{b}_1\left(t\right)y_3^2-\frac{1}{6}{c}_1\left(t\right)y_3^3+\frac{1}{6}{y}_3.\hfill \end{array}$$

(127)

We substitute the expressions (127) into the transformation (65) and obtain

$$\begin{array}{cc}\hfill B_1=& C_1=b_2\left(t\right)+c_1\left(t\right)y_3,\hfill \\ \hfill B_2=& C_2=\frac{1}{4}{c}_1\left(t\right)y_3^3+\frac{1}{4}{b}_1\left(t\right)y_3^2+\frac{1}{4}{y}_3-c_1\left(t\right)y_3-b_2\left(t\right)+1,\hfill \\ \hfill B_3=& C_3^{-1}=\frac{1}{C_3}=\frac{1}{-\frac{1}{6}{b}_1\left(t\right)y_3^2-\frac{1}{6}{c}_1\left(t\right)y_3^3+\frac{1}{6}{y}_3}.\hfill \end{array}$$

(128)

Finally, we substitute the expressions (128) into the transformation of (64) and obtain

$$\begin{array}{cc}\hfill A_1=& B_1{B}_3=\frac{b_2\left(t\right)+c_1\left(t\right)y_3}{-\frac{1}{6}{b}_1\left(t\right)y_3^2-\frac{1}{6}{c}_1\left(t\right)y_3^3+\frac{1}{6}{y}_3}=-\frac{6\left(b_2\left(t\right)+c_1\left(t\right)y_3\right)}{c_1\left(t\right)y_3^3+b_1\left(t\right)y_3^2{-y}_3},\hfill \end{array}$$

(129)

$$\begin{array}{cc}\hfill A_2=& B_2{B}_3=\frac{\frac{1}{4}{c}_1\left(t\right)y_3^3+\frac{1}{4}{b}_1\left(t\right)y_3^2+\frac{1}{4}{y}_3-c_1\left(t\right)y_3-b_2\left(t\right)+1}{-\frac{1}{6}{b}_1\left(t\right)y_3^2-\frac{1}{6}{c}_1\left(t\right)y_3^3+\frac{1}{6}{y}_3}\hfill \end{array}$$

(130)

$$\begin{array}{cc}\hfill =& -\frac{\frac{3}{2}{c}_1\left(t\right)y_3^3+\frac{3}{2}{b}_1\left(t\right)y_3^2+\frac{3}{2}{y}_3-6c_1\left(t\right)y_3-6b_2\left(t\right)+6}{c_1\left(t\right)y_3^3+b_1\left(t\right)y_3^2{-y}_3},\hfill \\ \hfill A_3=& B_3=-\frac{6}{c_1\left(t\right)y_3^3+b_1\left(t\right)y_3^2{-y}_3}.\hfill \end{array}$$

(131)

In Figure 51 and Figure 52, show the curves (129) and (131) for values $y_3=1$ and $y_3=-1$, respectively.

6. Conclusions

In the paper we show that all parametric expansions of variety $\Omega$ near its singularities and infinity can be computed with any accuracy and compute their first terms. We consider a very rich set of cases and find the ways to finish computations in all of them.

We do not intend to explain our results for original problems of Ricci flows and Einstein’s metrics. Let it will be done by authors of [14,15,16,17,18,19,20,21,22], who are specialists in the problem.