1.

There are two universal methods for local study of nonlinear equations and systems of different kinds (algebraic, ordinary and partial differential): (a) normal form and (b) truncated equations.

(a) Equations with linear parts can be reduced to their normal forms by local changes of coordinates. For algebraic equation, it is Implicit Function Theorem. For systems of ordinary differential equations (ODE), I completed the theory of normal forms, began by Poincaré (1879) [Poincaré, 1928] and Dulac (1912) [Dulac, 1912] for general systems [Bruno, 1964;1971] and began by Birkhoff (1929) [Birkhoff, 1966] for Hamiltonian systems [Bruno, 1972;1994].

(b) Equations without linear part: I proposed to study properties of solutions to equations (algebraic, ordinary differential and partial differential) by studying sets of vector power exponents of terms of these equations. Namely, to select more simple ("truncated") equations [Bruno, 1962;1989;2000] by means of generalization to polyhedrons the Newton (1678) [Newton, 1964] and the Hadamard (1893) [Hadamard, 1893] polygons.

By means of power transformations [Bruno, 1962;1989;2022b] the normal forms and the truncated equations can be strongly simplified and often solved. Solutions to the truncated equations are asymptotically the first approximations of the solutions to the full equations. Continuing that process, we can obtain then (2.1)

London Journal of Research in Science: Natural and Formal Theorem 2.1. Let f (X, ?, T ) = ?a Q,r (T )X Q ? r , where 0 ? Q ? Z n , 0 ? r ? Z, the sum is finite and a Q,r (T ) are some functions of T = (t 1 , . . . , t m ), besides a 00 (T ) ? 0, a 01 (T ) ? ? 0. Then the solution to the equation f (X, ?, T ) = 0 has the form

where 0 ? R ? Z n , 0 < ?R?, the coefficients b R (T ) are functions on T that are polynomials from a Q,r (T ) with ?Q? + r ? ?R? divided by a

2?R?-1 01

. The expansion b(T, X) is unique. Let g(X, ?, T ) = f (X, ? + b(T, X), T ), (2.2) then g(X, 0, T ) ? 0.

This is a generalization of Theorem 1.1 of [Bruno, 2000, Ch. II] on the implicit function and simultaneously a theorem on reducing the algebraic equation (2.1) to its normal form (2.2) when the linear part a 01 (T ) ? ? 0 is nondegenerate. In it, we must exclude the values of T near the zeros of the function a 01 (T ).

Let the point X 0 = 0 be singular. Write the polynomial in the form

The set S is called the support of the polynomial f (X). Let it consist of points Q 1 , . . . , Q k . The convex hull of the support S(f ) is the set (2.3) which is called Newton's polyhedron.

Its boundary ?Î?"(f ) consists of generalized faces Î?" (d) j , where d is its dimension of 0 ? d ? n -1 and j is the number. Each (generalized) face Î?" ( d ) j corresponds to its:

? boundary subset S (d)

? and normal cone

, (2.4) where P = (p 1 , . . . , p n ) ? R n Let X = (x 1 , . . . , x n ) ? R n or C n , and f (X) be a polynomial. A point X = X 0 , f (X 0 ) = 0 is called simple if the vector (?f /?x 1 , . . . , ?f /?x n ) in it is non-zero. Otherwise, the point X = X 0 is called singular or critical. By shifting X = X 0 + Y we move the point X 0 to the origin Y = 0. If at this point the derivative ?f /?x n ? = 0, then near X 0 all solutions to the equation f (X) = 0 have the form y n = ?b q 1 ,...q n-1 y q

Nonlinear Analysis as a Calculus

At X ? 0 solutions to the full equation f (X) = 0 tend to non-trivial solutions of those truncated equations f (d) j (X) = 0 whose normal cone U (d) j intersects with the negative orthant P ? 0 in R n * .

Remark 1. If in the sum (2.1) all Q belong to a forward cone C:

then in the solution (2.2) of Theorem 2.1 all R belong to the same cone C: [Bruno, 1989, Part I, Chapter 1, § 3].

Let ln X def =(ln x 1 , . . . , ln x n ). The linear transformation of the logarithms of the coordinates

(2.5) [Bruno, 1962;2000: where ? is a nondegenerate square n-matrix, is called power transformation.

By the power transformation (2.5), the monomial X Q tranforms into the monomial Y R , where R = Q (? * ) -1 and the asterisk indicates a transposition.

A matrix ? is called unimodular if all its elements are integers and det ? = ±1. For an unimodular matrix ?, its inverse ? -1 and transpose ? * are also unimodular.

Theorem 2.2. For the face Î?" (d) j there exists a power transformation (2.5) with the unimodular matrix ? which reduces the truncated sum

) is a polynomial. Here S ? Z n . The additional coordinates y d+1 , . . . , y n are local (small).

The article [Bruno, Azimov, 2023] specifies an algorithm for computing the unimodular matrix ? of Theorem 2.2.

be a face of the Newton polyhedron Î?"(f ). Let the full equation f (X) = 0 is changed into the equation g(Y ) = 0 after the power transformation of Theorem 2.2. Thus ?(d) j (y 1 , . . . , y d ) = g(y 1 , . . . , y d , 0, . . . , 0).

London Journal of Research in Science: Natural and Formal

Let the polynomial ?j be the product of several irreducible polynomials

where 0 < l k ? Z. Let the polynomial h k be one of them. Three cases are possible:

Case 1. The equation h k = 0 has a polynomial solution y d = ?(y 1 , . . . , y d-1 ). Then in the full polynomial g(Y ) let us substitute the coordinates

for the resulting polynomial h(y 1 , . . . , y d-1 , z d , y d+1 . . . , y n ) again construct the Newton polyhedron, separate the truncated polynomials, etc. Such calculations were made in [Bruno, Batkhin, 2012] and were shown in [Bruno, 2000, Introduction].

Case 2. The equation h k = 0 has no polynomial solution, but has a parametrization of solutions

Then in the full polynomial g ( Y ) we substitute the coordinates

where ? j = const, ? |? j | ? = 0, and from the full polynomial g(Y ) we get the polynomial

where

If in the expansion (5.7) l k = 1, then a 01 ? ? 0. By Theorem 2.1, all solutions to the equation h = 0 have the form i.e., according to (2.7) the solutions to the equation g = 0 have the form

Such calculations were proposed in [Bruno, 2018a].

If in (5.7) l k > 1, then in (2.8) a 01 (T ) ? 0 and for the polynomial (2.8) from Y ?? , ? we construct a Newton polyhedron by support S(h) = {Q ?? , r : a Q ?? ,r (T ) ? ? 0}, separate the truncations and so on.

Case 3. The equation h k = 0 has neither a polynomial solution nor a parametric one. Then, using Hadamard's polyhedron [Bruno, 2018a;2019a], one can compute a piece-wise approximate parametric solution to the equation h k = 0 and look for an approximate parametric expansion.

Similarly, one can study the position of an algebraic manifold in infinity.

Here we consider an ordinary differential equation of the form

where x is independent variable, y is the dependent variable, y ? = dy/dx and f is a polynomial of its arguments. Near x 0 = 0 or ? we look for solutions of equation (3.1) in the form of asymptotic series

III. SINGLE ODE [BRUNO, 2004] 3.1. Setting of the problem:

where b k are functions of log x and ?s k > ?s k+1 with

We set X = (x, y). By a differential monomial a(x, y) we mean the product of an ordinary monomial To every differential monomial a (X) one assigns its (vector) exponent Q(a) = (q 1 , q 2 ) ? R 2 by the following rules. For a monomial of the form (3.4) let r 1 ,r 2 ); for a derivative of the form (3.5) let Q d l y/dx l = (-l, 1).

When differential monomials are multiplied, their exponents are summed as vectors:

The set S(f ) of exponents Q (a i ) of all the differential monomials a 2 (X) in a differential sum of the form (3.6) is called the support of the sum f (X). Obviously, S(f ) ? R 2 . The closure Î?"(f ) of the convex hull of the support S(f ) is referred to as the polygon of the sum f (X). The boundary ?Î?"(f ) of the polygon Î?"(f ) consists of vertices Î?"

Let R 2 * be the plane conjugate to the plane R 2 so that the inner (scalar) product

is defined for any

and the truncated sum (3.7). All these constructions are applicable to equation (3.1), where f is a differential sum.

Let x ? 0 or x ? ? and suppose that a solution of the equation (3.1) has the form

where c r is a coefficient, c r = const ? C, c r ? = 0, the exponents r and ? are in R, and ?? < 0. Then we say that the expression

gives the power-law asymptotic form of the solution (3.8).

Thus, corresponding to any face

* and the truncated equation

Then the solution (3.7), (3.10) satisfies the equation

London Journal of Research in Science: Natural and Formal g(X) = 0

Substituting y = cx r into g(X), we see that g (x, cx r ) does not depend on x, c and is a polynomial in r, that is,

where ?(r) is the characteristic polynomial of the differential sum f (0) j (X). Hence, in a solution (3.9) of the equation (3.10) the exponent r is a root of the characteristic equation (3.11) and the coefficient c r is arbitrary. Among the roots r i of the equation (3.11), one must single out only those for which one of the vectors ? (1, r), where ? = ±1, belongs to the normal cone U (0)

In this case the value of ? uniquely determined. The corresponding expressions of the sum with an arbitrary constant c r are candidates for the role of truncated solutions of the equation (3.1). Moreover, by (3.3)

Complex roots r to characteristic equation (3.11) may bring to exotic expansions of solutions (3.2), where coefficients b k are power series in x ?i with real ? ? R and i 2 = -1.

(1)

j , this condition uniquely determines the exponent r of the truncated solution (3.9) and the value ? = ±1 in (3.3). To find the coefficient c r , one must substitute the expression (3.9) into the truncated equation (3.10). After cancelling a factor which is a power of x, we obtain an algebraic defining equation for the coefficient c r ,

r ? = 0 of this equation is an expression of the form (3.9) which is a candidate for the role of a truncated solution of the equation (3.1). Moreover, by (3.3), if in the normal cone U

(1) j one has p 1 < 0, then x ? 0, and if p 1 > 0, then x ? ?. From the polygon Î?" of the initial equation (3.1) we take a vertex or an edge Î?" (d) j . Then we found a power solution y = b 1 x P 1 of the truncated equation f (d) j (X) = 0, as it was described above, put

We construct the polygon 1 Î?" for the new equation, take a vertex or an edge 1 Î?" (e) k , solve the truncated equation

and obtain the second term b 2 x P 2 of expansion (3.2) and so on.

London Journal of Research in Science: Natural and Formal

We construct the polygon 1 Î?" for the new equation, take a vertex or an edge 1 Î?" (e)

k , solve the truncated equation

and obtain the second term b 2 x P 2 of expansion (3.2) and so on.

In [Bruno, 2004] there are some properties, that simplify computation. Thus, we can obtain the 4 types of expansions (3. [Bruno, 2006;2018b]; 4. Exotic, when all b k are power series in x i? [Bruno, 2007].

Except expansions (3.2) of solutions y(x) of equation (3.1), there are exponential expansions

where b k (x) and ?(x) are power series in x [Bruno, 2012a,b].

Also there are solutions in the form of transseries [Bruno, 2019b]. These results were applied to 6 Painlevè equations [Bruno, 2015;2018b,c;Bruno, Goruchkina, 2010]. Written as differential sums they are: Equation P 5 :

Here a, b, c, d and ?, ?, ?, ? are complex parameters. If all they are nonzero, then polygons for these equations are shown in Figures 1,2,3.

Nonlinear Analysis as a Calculus

Then there exists such power series ?(x) with integral increasing exponents, that after substitution

has only resonant terms b m x m , where

and m ? ?.

So here the eigenvalue

3) truncated differential sum f (0) 1 (X) have eigenvalues ? 1 , . . . , ? l , 0 ? l ? n; 4) the most left point of the support S(f ) in the axis q 2 = 0 be (?, 0). Evidently ? ? Z.

Supports and polygons for equations P 3 (left), P 4 (right).

Nonlinear Analysis as a Calculus

Theorem 3.3. Let 1) f x, y, y ? , . . . , y (n) be a polynomial in x, y, y ? , . . . , y (n) ; 2) its Newton polygon Î?"(f ) have a vertex Î?" (0) 1 = (v, 1) at the right side of its boundary ?Î?";

3) truncated differential sum f (0) j (X) have eigenvalues ? 1 , . . . , ? l , 0 ? l ? n; 4) the most right point of the support S(f ) in the axis q 2 = 0 be (?, 0). Evidently ? ? Z.

Then there exists such power series ?(x) with integral decreasing exponents, that after substitution (3.13), the differential sum (3.14) for identities (3.15) has only resonant terms b m x m , where equality (3.16) is true, and m ? ?. f (0) j (X) has no integral eigenvalue ? k ? ?-v (for Theorem 3.2) or ? k ? ? -v (for Theorem 3.3), then the initial equation f (X) = 0 has formal solution y = ?(x). If the truncated sum f (0) j (X) contains the derivation y (n) , then the series ?(x) converges according to Theorem 3.4 in [Bruno, 2004].

Remark 2. If the truncated sum f (0) j (X) has integral eigenvalue ? k ? ?v (for Theorem 3.2) or ? k ? ? -v (for Theorem 3.3), then the initial equation f (X) = 0 Supports and polygons for equations P 5 (left), P 6 (right).

Nonlinear Analysis as a Calculus

We will consider such a generalization of the power function cx r which preserves their main properties. The real number

where arg x = const ? [0, 2?), is called the order of the function ?(x) on the ray when x ? 0 or x ? ?. The order p ? (?) is not defined on the ray arg x = const, where the limit point x = 0 or x = ? is a point of accumulation of poles of the function ?(x).

In Subsections 3.2-3.4 it was shown that as x ? 0 (? = -1) or as x ? ? (? = 1) solutions y = ?(x) to the ODE f (x, y) = 0, where f (x, y) is a differential sum, can be found by means of algorithms of Plane PG, if

where n is the maximal order of derivatives in f (x, y). Here we introduce algorithms, which allow calculate solutions y = ?(x) with the property

Note, that in Plane PG we had ? ? = -1, i. e. ? +?? ? = 0. So, new interesting possibilities correspond to ? + ?? ? > 0.

We consider the ODE

where f (x, y) is a differential sum. To each differential monomial a i (x, y), we assign its (vector) power exponent Q(a i ) = (q 1 , q 2 , q 3 ) ? R 3 by the following rules: power exponent of the product of differential monomials is the sum of power exponents of factors:

The set S(f ) of power exponents Q(a i ) of all differential monomials a i (x, y) presented in the differential sum f (x, y) is called the space support of the sum f (x, y). Obviously, S(f ) ? R 3 . The convex hull Î?"(f ) of the support S(f ) is called the polyhedron of the sum f (x, y). The boundary ?Î?"(f ) of the polyhedron Î?"(f ) consists of the vertices Î?" (0) j , the edges Î?"

Support and polyhedron for equation P 1 . The approach allows to obtain solutions with expansions (3.2), where coefficients b k (x) are all periodic or all elliptic functions [Bruno, 2012c,d;Bruno, Parusnikova, 2012].

Expansions of solutions to more complicated equations such as hierarchies Painlevé see in [Anoshin, Beketova, (et al.), 2023; Bruno, For P 1 -P 5 with all parameters nonzero, their polyhedrons are shown in Figures 4,5,6,7,8 correspondingly.

Here we consider the system ?i = f i (X), i = 1, . . . , n, (4.1) where

Let the point X 0 = 0 be a singular point. Then the system (4.1) has the linear part

where A is a square n-matrix. Let ? = (? 1 , . . . , ? n ) be a vector of its eigenvalues.

Theorem 4.1 ([Bruno, 1964;1971, 1972]). There exists an invertible formal change of coordinates

where ? i (Y ) are power series from Y = (y 1 , . . . , y n ) without free terms, which reduces the system (4.1) to normal form

IV. AUTONOMOUS ODE SYSTEM

Support and polyhedron for equation P 2 . Here y i g i (Y ) are power series on Y without free terms.

and

Then the number k of linearly independent Q ? N satisfying the equation ( 4.3) is called multiplicity of resonance.

Theorem 4.2. Let k be the multiplicity of resonance of the system (4.1). Then there exists a power transformation

with unimodular matrix ? which reduces the normal form (4.2), ( 4.3) to the system

in which the first k coordinates form a closed subsystem without a linear part, and the remaining nk coordinates are expressed via them by means of integrals.

Thus, if ? ? = 0, then the original system (4.1) of order n can be reduced to a system of order k, but without the linear part. Support and polyhedron for equation P 3 .

Let's write the system (4.1) as (4.4) and put A Q = (a 1Q , . . . , a nQ ).

The set

is called the support of the system (4.4). Its convex hull Î?" (2.3) is its Newton's polyhedron. Its boundary ?Î?" consists of generalized faces Î?"

According to [Bruno, 2000, Chapt. 1, §2] let d > 0 and Q be the interior point of a face Î?" (d) j , that is, Q does not lie in a face of smaller dimension. If d = 0, then 4.2: Newton's polyhedron [Bruno, 1962;2000]. Support and polyhedron for equation P 4 .

R ? Z n , which reduce the system (4.4) to the form

where the system

corresponds to the truncated system (4.5). be singular for the truncated system (4.7). Near the point (4.8), the local coordinates are

Let at the point Z = (z 1 , . . . , z n ) = 0 the eigenvalues of the matrix of the linear part of the system (4.7) are ? = ?1 , . . . , ?n , where ?1 , . . . , ?d are the eigenvalues of the subsystem of the first d equations.

Theorem 4.4. There exists an invertible formal change of coordinates [Bruno, 2022b]: where W = (w 1 , . . . , w n ) which reduces the system (4.6) to the generalized normal form

where

(4.10)

The system (4.9), (4.10) is reduced to a system of lower order by the power transformation of Theorem 4.2 .

Let X = X 0 be a singular point of the sys tem (4.1). Two cases are possible: Case 1. ? ? = 0, then by Theorem 4.1 we reduce the system to a normal form, then by Theorem 4.2 we reduce the normal form to a subsystem of order k < n without linear part and obtain the problem of studying its singular points.

Case 2. ? = 0, then we compute the Newton polyhedron and separate truncated systems in which the normal cone U (d) j intersects the negative orthant of P ? 0. Each of them is reduced to the form (4.6), (4.7) by the transformation of Theorem 4.3. For each singular point (4.8), we apply Theorem 4.4 and obtain a subsystem of smaller order.

Continuing this branching process, after a finite number of resolution of singularities we come to an explicitly solvable system from which we can understand the nature of solutions of the original system. But Theorem 4.3 can be applied to the original system (4.1), i.e. to each of the generalized faces Î?" (d) j of its Newton polyhedron Î?". Then to each singular point (4.8) we apply Theorems 4.4, 4.2 and reduce the order of the system. Here also through a finite number of steps of the singularity resolution we come to an explicitly solvable system. This allows us to study the singularities of the original system in infinity. This is the basis of the integrability criterion in [Bruno, Enderal, 2009;Bruno, Enderal, Romanovski, 2017].

The normal form can be computed in the neighborhood of a periodic solution or invariant torus [Bruno, 1972, II, §11], [Bruno, 2022a].

See [Bruno, Batkhin, 2023] for similar computations for a system of partial differential equations.

and is defined by one Hamiltonian function H(x, y), where x = (x 1 , . . . , x n ), y = (y 1 , . . . , y m ). Here the normal form of the system (4.11) corresponds to the normal form of one Hamiltonian function. See details in [Bruno, Batkhin, 2021].

Let X = ( x 1 , . . . , x n ) ? C n or R n independent variables and y ? C or R be a dependent one. Consider Z = (z 1 , . . . , z n , z n+1 ) = (x 1 , . . . , x n , y).

Differential monomial a(Z) is called a product of an ordinary monomial

where c = const, and a finite number of derivatives of the following form

Vector power exponent Q(a) ? R n+1 corresponds to the differential monomial a(Z), it is constructed according to the following rules:

The product of monomials corresponds to the sum of their vector power exponents:

Differential sum is the sum of differential monomials V. ONE PARTIAL DIFFERENTIAL EQUATION 5.1. Support [Bruno, 2000 Ch. be 6-8]:

Let the support S(f ) of the differential sum (5.1) consists of one point E n+1 = (0, . . . , 0, 1). Then the substitution y = cX P , P = (p 1 , . . . , p n ) ? R n

(5.2) in the differential sum f (Z) gives the monomial c? P (P )X P

where ? P (P ) is a polynomial of P which coefficients depend on P .

Monomial (5.2) will be called resonant for f (Z) if for it ? P (P ) = 0.

Let µ k be the maximal order of the derivative over x k in f (Z), k = 1, . . . , n.

where ?(P ) is the characteristic polynomial of the sum of f (Z) and its coefficients do not depend on P . But if the inequalities (5.3) are not satisfied, then ? P (P ) ? = ?(P ).

For a differential sum f (Z) we denote by f k (Z) the sum of all differential monomials of f (Z) which have n + 1 coordinate q n+1 of vector power exponents Q = (q 1 , . . . , q n , q n+1 ) equal to k:

Consider the PDE f (Z) = 0.

(5.4)

5.2. Resonant monomials:

1. f 0 (Z) = ?(X) is a power series from X without a free term, 2. f 1 (Z) = a(Z)+b(Z), where S(a) = E n+1 = (0, . . . , 0, 1), S(b) ? Z n+1 + \0 × {q n+1 = 1}.

Then there exists a substitution y = ? + (X), where (X) is a power series from X without a free term, which transforms the equation (5.4) to the normal form g(X, ?) = 0,

(5.5)

where g 0 (X) = c P X P is a power series without a free term, P ? Z n + containing only resonant monomials c P X P for sum a(Z). is the formal solution to the equation (5.4).

If in equation (5.4) differential sum does not contain derivatives, then

Closure of a convex hull where the space R n+1 * is conjugate to the space R n+1 , ??, ?? is the scalar product, and truncated sum

Consider the equation

where f is the differential sum. In the solution of equation (5.6)

where ? is a series on the powers of x k and their logarithms, the series ? corresponds to its support, polyhedron, normal cones u i and truncations. The logarithm ln x i has a zero power exponent on x i . The truncated solution y = ? corresponds to the normal cone

Theorem 5.2. If the normal cone u intersects with the normal cone (5.2), then the truncation y = ?(X) of the solution (5.3) satisfies the truncated equation

(5.8)

To simplify the truncated equation (5.8), it is convenient to use a power transformation. Let ? be a square real nondegenerate block matrix of dimension n + 1 of the form

where ? 11 and ? 22 are square matrices of dimensions n and 1, respectively. We denote ln Z = (ln z 1 , . . . , ln z n+1 ), and by the asterisk * we denote transposition.

Variable change. ln W = (ln Z) ? (5.9) is called the power transformation.

Theorem 5. 3 ([Bruno, 2000]). The power transformation (5.5) reduces a differential monomial a(Z) with a power exponent Q(a) into a differential sum b(W ) with a power exponent Q(b):

London Journal of Research in Science: Natural and Formal

Corollary 5.3.1. The power transformation (5.9) reduces the differential sum (2.1) with support S(f ) to the differential sum g(W ) with support S(g) = S(f )? -1 * , i.e.

Theorem 5.4. For the truncated equation

there is a power transformation (5.9) and monomial Z T that translates the equation above into the equation

where g(W ) is a differential sum whose support has n + 1d zero coordinates.

Let z j be one of the coordinates x k or y. Transformation ? j = ln z j is called logarithmic.

Theorem 5.5. Let f (Z) be a differential sum such that all its monomials have a jth component q j of the vector exponent of degree Q = (q 1 , . . . , q m+n ) equal to zero, then the logarithmic transformation (5.1) reduces the differential sum f (Z) into a differential sum from z 1 , . . . , ? j , . . . , z n .

A truncated equation 5.7. Calculating asymptotic forms of solutions: f (n) j (Z) = 0 is taken. If it cannot be solved, then a power transformation of the Theorem 5.4 and then a logarithmic transformation of the Theorem 5.5 should be performed. Then a simpler equation is obtained. In case it is not solvable again, the above procedure is repeated until we get a solvable equation. Having its solutions, we can return to the original coordinates by doing inverse coordinate transformations. So the solutions written in original coordinates are the asymptotic forms of solutions to the original equation (5.2).

In [Bruno, Batkhin, 2023] method of selecting truncated equations was applied to systems of PDE.

Traditional approach to PDE see in [Oleinik, Samokhin, 1999;Polyanin, Zhurov, 2021].

Here we provide a list of some applications in complicated problems of (c) Mathematics, (d) Mechanics, (e) Celestial Mechanics and (f) Hydromechanics.

(c) In Mathematics: together with my students I found all asymptotic expansions of five types of solutions to the six Painlevé equations (1906) [Bruno, 2018c;Bruno, Goruchkina, 2010] and also gave very effective method of determination of integrability of ODE system [Bruno, Enderal, 2009;Bruno, Enderal, Romanovski, 2017].

(d) In Mechanics: I computed with high precision influence of small mutation oscillations on velocity of precession of a gyroscope [Bruno, 1989] and also studied values of parameters of a centrifuge, ensuring stability of its rotation [Batkhin, Bruno, (et al.), 2012].

(e) In Celestial Mechanics: together with my students I studied periodic solutions of the Beletsky equation ( 1956) [Bruno, 2002;Bruno, Varin, 2004], describing motion of satellite around its mass center, moving along an elliptic orbit. I found new families of periodic solutions, which are important for passive orientation of the satellite [Bruno, 1989], including cases with big values of the eccentricity of the orbit, inducing a singularity. Besides, simultaneously with [Hénon, 1997], I found all regular and singular generating families of periodic solutions of the restricted three-body problem and studied bifurcations of generated families. It allowed to explain some singularities of motions of small bodies of the Solar System [Bruno, Varin, 2007]. In particular, I found orbits of periodic flies round planets with close approach to the Earth [Bruno, 1981].

(f) In Hydromechanics: I studied small surface waves on a water [Bruno, 2000, Chapter 5], a boundary layer on a needle [Bruno, Shadrina, 2007], where equations of a flow have a singularity, and an one-dimensional model of turbulence bursts [Bruno, Batkhin, 2023].