D -Stability of the Initial Value Problem for Symmetric Nonlinear Functional Differential Equations

: This paper presents a method of establishing the D -stability terms of the symmetric solution of scalar symmetric linear and nonlinear functional differential equations. We determine the general conditions of the unique solvability of the initial value problem for symmetric functional differential equations. Here, we show the conditions of the symmetric property of the unique solution of symmetric functional differential equations. Furthermore, in this paper, an illustration of a particular symmetric equation is presented. In this example, all theoretical investigations referred to earlier are demonstrated. In addition, we graphically demonstrate two possible linear functions with the required symmetry properties.


Introduction
Time-reversal symmetry is one of the fundamental symmetries discussed in natural science. Consequently, it arises in many physically motivated dynamical systems, in particular, in classical and quantum mechanics [1]. There are various approaches which have attracted our interest that utilize the modeling of reversible dynamical systems, such as the study of time-reversal symmetry in nonequilibrium statistical mechanics [2], the time reversibility investigation in nonequilibrium thermodynamics [3], and the evaluation of certain properties of the Liapunov spectrum for the driven Lorentz gas [4].
The investigation of the impact of symmetries in functional differential equations is a topical problem. It is known that the well-posed problem or stability of solutions is one of the most important properties of functional differential equations. There are many publications that study this problem (herein, we only cite some of them: see, for example, [5][6][7][8][9][10] and the references therein). Using the results from [5,6,[11][12][13][14], we investigate the D-stability of the symmetric solutions for functional differential equations. Here, the general symmetry of the solution (see (3)) is studied, which is a generalization of the periodicity of solutions. A similar problem is studied, for example, in [11,[15][16][17].
The main focus of our research was to obtain the conditions of D-stability of functional differential equations. A similar topic was studied by Azbelev et al. in [5,6]. However, much more research is devoted to differential equations with delays (see, for example, [7][8][9]18,19] and the references therein).
Moreover, the big amount of the theorems in this field is based on Schauder's theorem or Banach's theorem of contractive mapping for the unique solvability or the Bohl-Perron theorem to prove the exponential stability, which imposes serious restrictions. In the present paper, an alternative research method that can be applied to a wide range of differential equations is described. It allows one to obtain not only unique solvability conditions but also a representation of the unique solution in the form of a functional series. Then, the symmetry conditions of the unique solution can be easily obtained. The principles of the obtained results for linear functional differential equations, the conditions of D-stability, and the symmetry of nonlinear equations were established. Thanks to the symmetry conditions, the examination of the equations is now possible on the full R axis, using the theory of boundary value problems. We see the limitations of this method only in its application to positive operators.
The present paper was mainly motivated by several papers that deal with the conditions of stability and solvability of differential equations with delay. For example, the solution estimate for a differential equation with delay x (t) + ∑ m k=1 b k (t)x((h k (t)) = 0, t ≥ t 0 , x(t) = φ(t), t ≤ 0, k = 1, . . . , m was obtained in [7]. New explicit conditions of exponential stability for the nonautonomous linear equation x (t) + ∑ m k=1 a k (t)x((h k (t)) = 0 were studied in [8]. In [9], they deal with the influence of "mistakes" in coefficients and delays on a solutions' behavior in the delay differential neutral system . . , n, t ∈ [0, ∞). Our investigation of functional differential equations with symmetric properties starts in [11,12]. They present the applicability of the theory of boundary value problems in the investigation of functional differential equations with symmetric properties on R. Then, the unique solvability conditions of linear functional differential equations were established in [13,14,20], which represent the unique solution in view of the functional series.

Problem Formulation
A class of symmetric solutions of the scalar nonlinear functional differential equations is considered here where t ∈ R, f : R m+2 → R is continuous, m ≥ 1, µ i , ν i , τ i : R → R are measurable functions, and p i , g i , q ∈ L(R, R), i = 1, 2, . . . , m. (1), we understand an absolutely continuous function x : R → R on every compact interval which satisfies (1) almost everywhere.

Definition 1 ([5]). By the solution of the Equation
The main focus of this paper was to obtain the D-stability conditions (see Definition 4) for the unique solution for Equation (1), with initial value condition and with the symmetric property where ψ is a monotonously increasing C 1 -function. The value of alpha in Formula (2) is selected by Condition (8). A typical example is a periodicity with ψ(t) = t + p. Because in reality there is no periodic symmetry, since there are always small noises, we demonstrate our theory in an example with a perturbation of periodic symmetry in Section 8. Of course, we develop our results for the more general ψ(t) introduced above. The paper is structured as follows: the main result of this investigation is presented in Section 7, in which we establish the D-stability conditions for the trivial solution of the nonlinear scalar initial value Problems (1) and (2) with the symmetric property (3). Notation and definitions can be found in Sections 3 and 5, respectively. In Section 4, we summarize the conditions from [11,12] under which the symmetric functional differential Equation (1) setting on R can be studied using boundary value theory. In Section 6.1, we establish the conditions for a unique solution for the initial value linear symmetric scalar problem, where we also find conditions under which the unique solution is symmetrical (see Section 6.2).

•
R is the space of real numbers with norm | · |;
From the problem formulation, it is clear that a restriction y = x I ψ of every solution x on the interval I ψ := [t 0 , ψ(t 0 )] satisfies a two-point boundary value condition Introducing the following notation and assumptions are necessary for further investigation: 1.
The increasing function ψ generates an increasing numerical sequence 2. Every point from Sequence (9) from (1) divides R on a specified quantity of intervals 3.
Assume that the number j is a number of the interval (10).

Definition 2 ([12]
). For every t ∈ R, we define the number l(t) by a number of the interval (10) that contains the point t.
Taking into account the definition of the function l : R → Z, then the next lemma is true.
Then, the function x : R → R defined by (11) is a solution to the problem ( (1) and (3)).

Definitions of the D-Stability
Taking into account (12)- (14), and Lemma 3, let us apply Equation (1) for t ∈ I ψ Definition 3 ([5]). The linear Equation is D-stable if Problem (18), (8) have a unique solution y ∈ D for arbitrary q ∈ L, α ∈ R and this solution continuously depends on {q, α}; this mean that for arbitrary ε > 0, there exists δ > 0 such that Let us consider the nonlinear functional differential Equation where operator N : D → L and N (0) = 0.

Definition 4 ([5]
). The nonlinear equation N y = 0 has a D-stable property in the vicinity of the trivial solution y ≡ 0 if there exist such δ 0 > 0 that for every pair {q, α} ∈ L × R, satisfying conditions q L < δ 0 , |α| < δ 0 , the problem (19) has a unique solution x ∈ D and for arbitrary ε > 0 there exists such

A Unique and Symmetric Solution for the Initial Value Problem for Linear Scalar Functional Differential Equations
Remark 1. Taking into account Lemmas 2 and 3, we study the initial-value problem (2) for nonlinear Equation (1) with the symmetric property (3) as a Problem (15), (8) in which operators ξ i , κ i , σ i and functions µ i , ν i , τ i have properties (12)- (14) and (4)-(7) on interval I ψ , respectively.

Unique Solvability
Here, we consider symmetric linear functional differential Equation (16) as a homogenous initial value problem and a nonhomogenous initial value problem where operators ξ i and κ i , i = 1, 2, . . . , m are defined by (12) and (13) with properties (4) and (5) and For the following investigation, we use the well-known result from the general theory of boundary value problems for the functional differential equation. (22) is uniquely solvable if the corresponding homogeneous Problem (21) for linear Equation (20) only has a trivial solution.

Lemma 4. The nonhomogeneous Problem (23) for linear Equation
In Section 6.1, theorems are established regarding the conditions of the existence of a trivial solution for the homogenous Problem (21) for linear Equation (20) and regarding a unique solution for the corresponding nonhomogeneous Problem (23) for linear Equation (22). Furthermore, in Section 6.2, we find conditions in α from (2) sufficient for the symmetric property of the unique solution from Section 6.1.
If we take (12) and (13) into account, then the next obvious Lemma is true.

Lemma 5. If measurable functions
have the properties (4), (5), and (24), and operators ξ i , κ i are defined by (12) and (13), then the nonhomogeneous Equation (22) is equivalent to the following equation: and the homogeneous Equation (20) is equivalent to the following equation: Let us specify an arbitrary absolutely continuous u 0 ∈ C(I ψ , R) and construct a sequence of functions u k , k = 0, 1, . . . defined by the recurrence relation where functions {ν i , µ i }, i = 1, 2, have the properties (25) and functions In the case where r = 1 and b 1 = 1, Equality (28) can be rewritten as where functions {ν i , µ i }, i = 1, 2, have properties (25). Obviously, Formula (30) defines the standard iteration sequence used in studies of the uniqueness of the trivial solution of the integral functional equation which, obviously, is equivalent to the homogeneous problem ( (27) and (21)).

Definition 5.
Consider the following equation: An operator l is said to be positive if for all t ∈ I ψ , the relation y(t) ≥ 0 implies that (ly)(t) ≥ 0 for almost every t ∈ I ψ .
have the properties (25) and are fulfilled, then the operator is positive.
Proof. Taking into account Definition 5 and (32), one can conclude that operator (33) is positive.
In Theorems 1, 3, and Corollaries 1-3, the conditions regarding the existence of trivial solution of the linear Cauchy problem (21) for the Equation (20) are established, and correspondingly, on a unique solution for the linear Cauchy problem (23) for the Equation (22) with a variable deviation (not only a delay deviation). Theorem 1. Suppose that (32) is true and functions µ i and ν i , i = 1, 2, . . . , m, have the properties (24) and (25). Assume also that one can specify some integers r and m, m ≥ r ≥ 1, a real number , and certain absolutely continuous vector-functions u 0 , u 1 , . . . , u r−1 satisfying Conditions (28) and (29), and the relation such that the functional differential inequality is satisfied for almost every t from I ψ , where, by definition for r ∈ N, m ≥ r, and k = 0, 1, . . . , m − 1.
Then, the homogenous linear initial value problem (21) for the Equation (27) only has a trivial solution and the nonhomogeneous linear Cauchy problem (23) for the Equation (26) is uniquely solvable for an arbitrary α ∈ R and an arbitrary function q ∈ L(I ψ , R). The unique solution of the problem (23) for the Equation (26) is representable in the form of a functional series is uniformly convergent on I ψ . Moreover, if α and the function q are such that for all t from I ψ , one has then, at every point t of the interval I ψ , the unique solution (37) of a problem (23) for the Equation (26) is non-negative.
satisfying Condition (29), and relation (34) such that the functional differential inequality is satisfied for almost every t from I ψ , where, by definition T r,m (k) is defined by (36) for r ∈ N, m ≥ r, and k = 0, 1, . . . , m − 1.
Then, the homogeneous linear Cauchy problem (21) for the Equation (31) (with q = 0) only has a trivial solution and the inhomogeneous linear Cauchy problem (23) for the Equation (31) is uniquely solvable for an arbitrary α ∈ R and an arbitrary function q ∈ L(I ψ , R). The unique solution of the problem (23) for the Equation (31) is representable in the form of a functional series uniformly convergent on I ψ , namely, whereq(t) is defined by (38). Moreover, if α and the function q are such that for all t from I ψ , one has (39), then, at every point t of the interval I ψ , the unique solution (41) of a problem (23) for the Equation (31) is non-negative. Now, let us put operator l by (33). Obviously, operator l in (33) is positive. Using Formula (33), it is easy to verify that, in this case, (35) coincides with Inequality (40). Applying Theorem 2, we obtain Theorem 1.

Proof.
To prove Corollary 1, we need Theorem 1. First of all, note that Sequence (43) is a particular case of (35) with a suitable choice of coefficients. More precisely, u 0 = z 0 , u 1 = z 1 , r = 2, m = 2 and We can apply Theorem 1. Indeed, let us set β k = β, k = 0, 1, 2 where β is a certain positive constant. Since, by (44), b 1 + b 2 = ρ, it is obvious that the relations (42) are satisfied. Inequality (35) in our case has the form (43).
a natural number m, non-negative integers k and r ≥ 1, and real numbers ρ, ρ > 1, c ∈ (0, 1) such that, for almost every t from the interval I ψ , the following inequality is satisfied: where u k is constructed by (30). Then, the assertion of Theorem 1 is true.
Theorem 4 (Theorem 2 [14]). Suppose that in the linear functional differential equation (31) with initial condition (23), operator l : C[I ψ , R] → L[I ψ , R] is positive and, moreover, there exists an absolutely continuous function u 0 : I ψ → R with property (45), a natural number m, non-negative integers k and r ≥ 1, and real numbers ρ, ρ > 1, c ∈ (0, 1) such that, for almost every t from the interval I ψ , the following inequality is satisfied: Then, the assertion of Theorem 2 is true.
It is easy to see that (46) is Inequality (47) with operator l defined by (33). Condition (32) ensures the positivity of the operator from Equations (20) and (22). Applying Theorem 4, we adhere to the conditions of Theorem 3.

Proof.
To prove Theorem 2, we need Theorem 3. Here, we have function u 0 defined by In this case, according to (49), Taking into account the recurrent sequence (30), we obtain (48) from (46). Applying Theorem 3 with k = 0 and r = 1, we obtain Corollary 2.
is fulfilled, where γ ∈ (0, 1) is a certain constant. Then, the assertion of Theorem 1 is true.
Proof. For γ ∈ (0, 1) and r ∈ N, we choose number c and ρ such that This can obviously be done. In view of (32) for µ i and ν i , i = 1, 2, . . . , m, with properties (24) and (25) for such c and ρ, the relation (51) yields Let us define u 0 by (49). In view of (50), (u 0 (t) ≥ 0) implies that the condition (46) is satisfied for the values of c and ρ specified above and the positivity of the operator (33) (in view of Lemma 6). It is easy to verify that Inequality (51) is a particular case of (46) for k = 0, r = 1, and, in this case, according to (30), we have Applying Theorem 3, we obtain Corollary 3.
There are many more other conditions on the unique solvability of the initial value problem for functional differential equations, for example, in [13,14,[20][21][22][23] and the references therein.

Symmetric Solution
The objective of our investigation was to find not only a unique solution for the linear problem (23) for the Equation (22) but also a symmetric solution (solution with Property (3)).
It is a well-known representation of the solution (37) of the problem (23) for the Equation (26) from where is the solution (37) of Equation (26) with q(t) = 0 and α = 1; and is the solution (37) of Equation (26) with α = 0 and arbitrary q(t) ∈ L. Assume that y 0 (ψ(t 0 )) = 1 (which means that this is not case of resonance), then from (8) and (52), we have where y 0 and Q(t) are defined by (53) and (54), respectively. If we have the resonance case y 0 (ψ(t 0 )) = 1, then this is Fredholm's alternative, and a solution exists only for q(t) such that Q(ψ(t 0 )) = 0, which is not unique and the kernel is 1-dimensional. Now, we are ready to establish the next Theorem. and then, the problem (23) for the Equation (26) has a unique symmetric solution represented by (52) with (53), (54), and (55).

D-Stability of Nonlinear Functional Differential Equations
In this section, we assume that operator N from (19) has the representation where L : D → L is a linear bounded operator from (16) and F : D → L is a nonlinear operator from (17) and has the property F (0) = 0 and q ∈ L.
Proof. In view of Lemmas 1 and 3, the initial value problem ( (1) and (2)) is equivalent to the following equation: g i (s)y(ν i (s))+ + f ((σ 1 y)(s), (σ 2 y)(s), . . . , (σ m y)(s), y(s), s) + q(s) ds + α t ∈ I ψ . (63) From the conditions of our theorem, we have that for some positive (K + M) < 1 ψ(t 0 )−t 0 there exists such δ > 0 that for arbitrary y i ∈ D, y i D ≤ δ, i = 1, 2, the inequalities (59)-(61) are fulfilled. Taking into account (59)-(61) and the local Banach theorem regarding contractive mapping in a sphere there exists a unique solution y ∈ D for Equation (63) and, in view of Lemmas 1 and 3, correspondingly, there exists a unique solution to the problem ((1) and (2)). Furthermore, for the solutions y i , i = 1, 2, of this problem with the condition where we have that It follows the continuous dependence of the solution y of Equation (63) on a norm of the space D from r with r D < δ(1 − (ψ(t 0 ) − t 0 )(K + M)) and, consequently, for the problem (1) from q and α with q L < δ, |α| < δ. Thus, the theorem is proved.

Theorem 7.
Assume that the conditions from Theorem 1, Theorem 3, Corollary 1, Corollary 2, or Corollary 3 are true, and for any M > 0, there is a δ > 0 such that inequality (60) holds for all y 1 D ≤ δ, y 2 D ≤ δ, and (62) is also satisfied. Then, the problem ( (1) and (2)) is D-stable in the vicinity of the trivial solution (see Definition 4). If the assumptions of Theorem 5 are true, then the symmetric problem ( (1) and (3)) is D-stable in the vicinity of the trivial solution.

Unique Solvability of the Initial Value Problem for Linear Scalar Functional Differential Equations
Here, we establish conditions sufficient for the unique solvability of the initial value problem for the following linear scalar functional differential equation: The next Theorem is true.
Theorem 8. Suppose that (70) is true, and for a certain constant γ ∈ (0, 1) and d satisfying condition d > 0, the following inequality holds almost everywhere on [0, 2π]: Then, the linear initial value problem ((75) and (76)) is uniquely solvable for an arbitrary α ∈ R and an arbitrary function q ∈ L([0, 2π], R). The unique solution of the problem ((75) and (76)) is representable in the form of a functional series Moreover, if α and function q is such that, for all t from [0, 2π], one has t 0 q(v)dv ≥ −α, then, at every point t of the interval [0, 2π], the unique solution (78) of the problem ((75) and (76)) is non-negative.

Remark 4.
Note that function p(t) cannot be a constant function because p(t) = constant does not fulfil the symmetric property (7).

Symmetric Property of the Unique Solution of the Initial Value Problem for Linear Scalar Functional Differential Equations
Here, we investigate a unique solution (78) with symmetric property (72). Assuming that u 0 (2π) = 1 (which means that this is not a resonance case), then u 0 (2π) is defined by (53): then the next Theorem is true.  75) and (76)) has a unique symmetric solution represented by (78).
Proof. It is easy to see that Theorem 9 is a particular case of Theorem 5.

D-Stability of Nonlinear Functional Differential Equations
Remark 5. In further investigations, we only study the D-stability of the trivial solution of Equation (1), as more general cases come down to this one.
Taking into account Definition 4, the next obvious Corollary 5, from Corollary 4 and Theorem 8, is true.
Then, the nonhomogeneous problem (2) for the Equation (68) is D-stable in the vicinity of the trivial solution. If the assumptions of Theorem 5 are true, then the symmetric problem (69) for the Equation (68) is D-stable in the vicinity of the trivial solution.
Proof. It is easy to see that Theorem 10 is a particular case of the Theorem 7.

Conclusions
The aim of the present manuscript was to establish conditions on the D-stability of the initial value problem for a class of nonlinear symmetric functional differential equations. Furthermore, we construct conditions on the unique solvability of the initial value problem for a class of linear symmetric functional differential equations. In view of the symmetric properties and our previous results, we study the aforementioned problem on R using the theory of boundary value problems. The unique solution is represented in the form of a functional series. The conditions of the symmetric property of the unique solution of the symmetric functional differential equations are also given. In addition, we present the application of the obtained results in an example.
The question of the D-stability of functional differential equations with symmetries in R n remains open, i.e., considering the D-stability of symmetric solutions of higher dimensional functional differential equations satisfying Ax(t) = x(ψ(t)) for a matrix A.