On a Class of Orthogonal Polynomials as Corrections in Lienard Differential System: Applications

In this paper we demonstrate some specialized modules for investigating the dynamics of differential models, an integral part of a planned much more general Web-based application for scientific computing. As “corrections” in the Lienard differential system is presented a class of orthogonal polynomials (also known as “shell polynomials”). We will note that some specifics of the amplitudes of these polynomials open up the possibility of modeling signals from the field of antennafeeder techniques. Algorithms and modules have been consistently used for: automatic generation of a theorem on the number and type of limit cycles (in the light of Melnikov’s considerations); study of the Hamiltonian of the system and “level curves”; for the study of catastrophic surfaces (in the light of Zeeman’s considerations), etc. Similar studies have been carried out for associated polynomials. Numerical examples, illustrating our results using CAS MATHEMATICA are given.


Introduction
122 years ago, Hilbert [1] proposed 23 mathematical problems, of which the second part of the 16th one is to find the maximal number of limit cycles and their relative locations for polynomial vector fields. To date, several thousand articles and studies have been devoted to this problem. The study of dynamical systems includes bifurcation theory with branch catastrophe theory [2]. See also model in the light of Zeeman's approach [3]. Arnold [4] discussed the catastrophes of the ADE classification, because of their relation with the Lie groups. For more details of the existing important results on the generalized polynomial Lienard differential systems and the limit-cycle bifurcations of some generalized polynomial Lienard systems, see  (where the reader can to discover a substantial additional bibliography). Some of our previous research [31][32][33] on this issue encouraged us to begin developing specialized modules as part of a much more general Web-based application for scientific computing. In Section 2, we demonstrate some algorithms and modules for: automatically generating a theorem for the number and type of limit cycles (in the light of Melnikov's considerations); study of the Hamiltonian of the system and "level curves"; for the study of catastrophic surfaces (in the light of Zeeman's considerations); generation and simulation of antenna factor, etc. Concluding remarks are placed in Section 3.
In this article we consider the Lienard system [6] where 0 ≤ < 1; F(x) are specially chosen polynomials, and Poly(x) are the polynomials P n (x)-see Figure 1. The level curves are studied. Some applications of the polynomials are also given.

Main Results
Let a = 1. Then Let Poly(x) (in (1)) coincides with polynomials P n (x). Without going into details, we will note some more interesting level curves.  The case (2): n = 7. The Hamiltonian of system (1) ( = 0) is The level curves L h i = {H(x, y) = h i } are depicted at Figure 3. The model in the light of Melnikov's considerations The Melnikov polynomial [5] for the system is defined as It is known [7,8] that the system for sufficiently small = 0 has at most n limit cycles asymptotic to circles of radii r j , j = 1, 2, . . . , n as → 0 if and only if the nth degree polynomial M(r 2 , n) has n positive roots r 2 = r 2 j , j = 1, 2, . . . , n. The case n = 7. Consider the model The following is valid.
Proposition 1. The Lienard-type system for n = 7, and for all sufficiently small = 0.
Proof. For the Melnikov polynomial in r 2 (see Figure 4) we have: Evidently, for example µ = 0.0735880318901 we have a simple limit cycle and cycle with multiplicity two. The case n = 11. Consider the model The following is valid.
Proposition 2. The Lienard-type system for n = 11, and for all sufficiently small = 0.
Proof. For the Melnikov polynomial in r 2 (see Figure 5) we have: Evidently, for example µ = 0.0012953041893 we have three simple limit cycles and limit cycle with multiplicity two. Numerical methods for finding zeros of polynomials can be found in [38][39][40].

Some simulations
The simulations on the Lienard-type system:   We will note that some specifics of the amplitudes of these polynomials open up the possibility of modeling signals from the field of antenna-feeder technology.
It is easy to take into account that the change of the variable t with t = b cos θ + c (θ is the azimuthal angle and c is the phase difference) in the y(t)-component of the solution of the system (6) leads to radiation diagrams [41,42]. Consider the following model in the light of Zeeman's approach: with c > 0 and G(x) = x 9 − 4609 2520 x 7 + 10, 219 10, 080 The catastrophe surfaces (x, y, p) = G(x) − y (q = 5, 10, 15) for the model is depicted on Figure 8. Consider the model (7) with c > 0 and x 3 − qx.

Some simulations
The simulations on the Lienard-type system: where with x 0 = 0.8, y 0 = 0.6, b = 0.9, c = 0.7809 are depicted on Figure 11. All experiments and algorithms were carried out using our own module in CAS Mathematica.
We define the normalized factor as follows

Calculator Software Application
We are developing a high-scalable, cloud-based software calculator using serverless architecture [43]. The serverless architecture enables automatic scaling of the system during high load. Furthermore, it can be used to parallelize suitable computations for higher efficiency. Where possible, we employ various optimization techniques for highperformance calculations, including multi-processor and multi-threading calculations, and hardware intrinsics [44][45][46]. The system is exposing the implemented algorithms using industry-standard application programming interface using HTTP and REST, with data being serialized in JSON and XML formats.
We are developing mobile, native and web-based clients (intellectual property at this stage) to provide an end-user experience for researchers. The following capabilities has been implemented: the user sets: the polynomials Poly i (x) (see (1)), which can be arbitrary orthogonal polynomials and their associated polynomials (such as associated Hermite polynomials, associated Gegenbauer polynomials, associated Legendre polynomials, associated Lommel polynomials, q-Lommel polynomials associated with the Jackson q-Bessel function, continuous and bivariate q-Hermite polynomials, extended Gegenbauer polynomials and their q-analogues, associated Jacoby polynomials, Chebyshev and Gegenbauer polynomials of higher kind etc.); the function F(x).
The API can be used by reporting and analytics systems like PowerBI and Excel to further investigate the results [47].
The application provides an opportunity for research in two directions-the study of the dynamics of differential systems and the generation of special classes of radiation diagrams.
Some of the algorithms used in this paper

1.
In the Melnikov polynomials M(r 2 , n), the coefficients are fractional numbers, and with a high degree of the polynomial, the user is faced with solving a "classic problem with imprecise data". This requires as a first step the use of an algorithm to approximately find the multiplicities of the zeros of the polynomial, then proceed to an algorithm with a user-fixed rate of convergence to determine all the zeros of the polynomial.

2.
Specialized algorithm for detailed Hamiltonian study of system (1) and visualization of "level curves" (assuming implementation of software tools in a user-selected computer-algebraic system for scientific calculations).

3.
Algorithm for matching the initial approximations when solving the differential system (1), given its interesting specificity and behavior of the solution in confidential time intervals. 4.
Algorithm for control and visualization of the "antenna factor" (with a possibly user-set value of the lateral radiation).

Concluding Remarks
Some of our previous research on this issue encouraged us to start developing specialized modules, part of a much more general Web-based application for scientific computing. We mention the above algorithms because some of them are hidden from the user. Determining the number and type of cycles, as we have already mentioned, is a complicated task (with inaccurate data-Algorithm 1). Our proposed module automatically generates theorems in light of Melnikov's considerations (see e.g., Propositions 1 and 2). This is very important for the user to take further steps in the detailed study of the dynamic model (for example, "level curves"-Algorithm 2) with the corrective corrections set by him in the Lienard differential system of the type of arbitrary high-order orthogonal polynomials or their associated or appropriate classes of Morse-type polynomials, etc. We will explicitly note that the user does not have to be a mathematician! What information about it would be, for example, the standard comment from existing computer algebra platforms: "The solution of the differential system lacks consistency under your chosen set of initial approximations". After all, we have to provide the user with a satisfactory solution to the task (typical example-the hidden Algorithm 3). Another algorithm hidden from the user is checking the conditions in Lienard's theorem for the existence of a limit cycle at all! Another example of a hidden algorithm is the recurrent generation of the polynomials under study at a user-fixed degree n. Algorithm 4 is extremely complex (and at this stage-insufficiently specified and developed). The use of y(θ)-the solution component of the corresponding Lienard differential system as an antenna factor is very complicated. As far as the Dolph-Chebyshev technique for synthesis of power pattern end filter prototypes is well known, we note that by analogy we can define hypothetical transmitting functions based on "van Doorn polynomials" (object of consideration in this article). The experts have a word! Of course, research on the "van Doorn array" can be carried out in the light of Soltis considerations [48], but this is the subject of future development and will be omitted here. We fully understand that the construction of such an ambitious Web-based platform for scientific computing can only be realized with the active participation of specialists from various branches of scientific knowledge.