Even-Order Neutral Delay Differential Equations with Noncanonical Operator: New Oscillation Criteria

: The main objective of our paper is to investigate the oscillatory properties of solutions of differential equations of neutral type and in the noncanonical case. We follow an approach that simpliﬁes and extends the related previous results. Our results are an extension and reﬂection of developments in the study of second-order equations. We also derive criteria for improving conditions that exclude the decreasing positive solutions of the considered equation


Introduction
The aim of this work is to provide new criteria for testing the oscillation of solutions to the even-order neutral differential equation (NDE) r(f) (υ(f) + p(f)(υ(τ(f)))) (n−1) γ + q(f)υ γ (ρ(f)) = 0, (1) in the noncanonical case, that is, where f ≥ f 0 , n ≥ 4 is an even integer, and γ is a quotient of odd positive integers.
Since the time of Newton, differential equations (DE) are still used to understand and model phenomena.Applied phenomena and problems are constantly increasing as a result of the development of all different branches of science.Delay differential equations (DDE) are DE that take into account the time memory of a solution, and therefore are a more suitable method for modeling different phenomena.However, the problem of finding solutions to the equations resulting from the modeling of phenomena stands as an obstacle to understanding and studying these phenomena.Therefore, the qualitative theory contributes greatly to solving this problem, so that the qualitative properties of equations can be studied without finding their solutions.One of the branches of the qualitative theory is oscillation theory, which is the theory that deals with the oscillatory, non-oscillatory, asymptotic behavior and the distribution of zeros for solutions of DE.
Fourth-order DEs often appear in mathematical models of many biological, chemical, and physical phenomena.Structure deformation, elasticity issues, and soil settlement are examples of these uses.The existence or nonexistence of oscillatory solutions is one of the most crucial topics that arise when investigating mechanical problems.Oscillatory muscle movement is one of the models given by a fourth-order equation with a delay that might arise as a result of a muscle's inertial pregnancy, see [1].Fourth-order equations arise in the theory of numbers, which is unusual, see [2].
NDEs appear as models for many phenomena and problems in applied science, including automatic control, population dynamics, mixing liquids, and vibrating masses attached to an elastic bar (see [3]).In particular, neutral equations of second order are of considerable interest in biology to study the self-balance of the human body and in robotics to construct biped robots (see [4]).The applications of these equations have largely stimulated the study of qualitative properties, in particular the theory of oscillation.
Recently, there has been a great interest in the oscillation theory of DDEs, see .This interest greatly contributed to the development of different new techniques and approaches that significantly improved the oscillation criteria, especially for second-order DE.Furthermore, this development has also resulted in many open and interesting issues.
We find that the study of the oscillatory behavior of the NDE in the canonical case (y 0 (f 0 ) = ∞) has attracted the greatest interest regarding the verification and study by many techniques and methods.However, it is possible to find some studies on the oscillation of even-order delay (nonneutral) DE with a noncanonical operator in the works by Baculikova et al. [5], Zhang et al. [6][7][8], and, recently, Moaaz et al. [9,10].
In 2011, Zhang et al. [6] presented the conditions of oscillation for DDE converge to zero, where α ≤ γ and α is a ratio of odd positive integers.Later, in 2013, Zhang et al. [7] established three independent conditions for oscillation of (2).In [8], under four independent conditions, Zhang et al. obtained some oscillation criteria for (2) when n = 4 and α = γ.The results in [8] differ from the results in [6,7] in that they apply to ordinary DE.By using the comparison theorems, Baculikova et al. [5] investigated the oscillatory propertiws of (2).
The two most common methods for studying oscillation of DDEs are: Riccati substitution and comparison with first-order equations.For NDEs in the noncanonical case, Li and Rogovchenko [18] obtained criteria for oscillation of the NDE (1) using the principle of comparison.They mainly relied on the assumption of the existence of the functions ξ i that fulfill the restrictions: and in addition the delay functions were restricted by the conditions Using these conditions, they were able to compare (1) with three DDEs of first order, and thus obtained three independent conditions to check oscillation.
NDEs of higher order have not received the attention of researchers compared to equations with delay (not neutral) or equations of the second order.The studies that deal with the oscillatory behavior of higher-order neutral equations have many interesting open issues as well as additional restrictions on the functions.
This paper aims to study the oscillatory behavior of solutions to the NDE (1) by using several different methods.The motive of this study is to develop and improve the oscillation criteria for solutions of higher-order DE by: -obtaining oscillation criteria for (1) without requiring the existence of the unknown functions ξ i , and without requiring the condition in (3); -reducing the number of conditions that are sufficient to verify the oscillation of all solutions; -obtaining criteria that apply to ordinary DE.

Preliminaries
In this part, we provide definitions and elementary results in the literature that help us to present the main results.
and υ satisfies (1) on [f * , ∞).To facilitate calculations, we will denote the corresponding function of the solution υ by z Definition 2. A solution of any DDE is called oscillatory if it is neither positive nor negative, ultimately; otherwise, it is called nonoscillatory.

Lemma 1.
[Lemma 2.2.1, [19]] Suppose that w ∈ C n ([f 0 , ∞), (0, ∞)), and w with derivatives up to order κ − 1 of constant sign, and w , where L 2 > 0, L 1 and L 3 are constants.Then, the maximum value of H on R at υ Remark 1.In this work, we consider only the solutions which are not identically zero, eventually.Further, all functional equations and inequalities are supposed to hold eventually.

Simplified Criteria for Oscillation
When studying the oscillatory properties of solutions of DDEs, it is known that positive solutions must be classified according to the sign of their derivatives.Now, we assume that υ is an eventually positive solution of (1).Then, from the definition of z, we have that z(f) > 0, eventually.From the DE in (1) and taking into account that q(f) > 0, we have that r(f) z (n−1) (f) γ is a nonincreasing function.Furthermore, according to Lemma 1, we obtain the following three cases, eventually: (1): In order to facilitate the presentation of results, we will use the following notations: The set of all eventually positive solutions of (1); The set of all solutions υ ∈ U + with z satisfying case (i) above, i = 1, 2, 3; and and for k = 0, 1, . . ., n − 2, eventually.
Proof.Assume, for the sake of contradiction, that (1) has a nonoscillatory solution.
Without loss of generality, we let υ ∈ U + .Thus, υ belongs to one of the classes U + i , i = 1, 2, 3. Using Theorem 2.1 in [20], we find that the condition ( 8) is in opposition to υ ∈ U + 2 .So, υ belongs to either U + 1 or U + 3 .First, let us assume that υ ∈ U + 3 .From Lemma 4, (4) and ( 5) hold.Using (5) with k = n − 2, we obtain that Hence, from the definition of z and the above inequality, we see that Then, from this inequality and (1), we obtain Now, by integrating (10) twice from f 1 to f, we have Taking into account the monotonicity of r(f) z (n−1) (f) γ , we obtain that Then, from (4) with k = n − 2, we arrive at which in view of (11) implies that From (9) and the above inequality, we obtain that z (n−2) (f) → −∞ when f → ∞.This is a contradiction with the positivity of z (n−2) (f).
Proof.Assume, for the sake of contradiction, that (1) has a nonoscillatory solution.
Without loss of generality, we let υ ∈ U + .As in the proof of Theorem 1, using the condition (8), we have that υ belongs to either U + 1 or U + 3 .First, let us assume that υ ∈ U + 3 .In a similar fashion as in the proof of Theorem 1, we can arrive at (10).Then, integrating (10) from f 1 to f, we find Since z (f) < 0 and ρ(f) ≤ f, we find that From ( 4), with k = n − 2, in Lemma 4, we have that Combining ( 16) and ( 17), we conclude that This implies that which contradicts (14).Now, let us assume that υ ∈ U + 1 on [f 1 , ∞), where f 1 ≥ f 0 .From ( 14), there exists a

and thus
Letting f → ∞, we obtain that (13).As in the proof of Theorem 1, we arrive at (12).
Proof.Assume, for the sake of contradiction, that (1) has a nonoscillatory solution.
Without loss of generality, we let υ ∈ U + .As in the proof of Theorem 1, using the condition (8), we have that υ belongs to either U + 1 or U + 3 .First, let us assume that υ ∈ U + 3 .In a similar fashion as in the proof of Theorem 1, we can arrive at (15).From Lemma 4, (4) holds.By (4) with k = n − 3 and (15), we find that Hence, we note that w := z is a positive solution of the delay differential inequality It follows from Theorem 1 in [21] that there is also a positive solution of the DDE By Theorem 2 in [22], we find that the condition (18) guarantees oscillation of all solutions of (19), which is a contradiction.
Finally, let us assume that υ ∈ U + 1 .Since the condition ( 13) is necessary for the validity of (18), the proof of this part is exactly the same as in Theorem 1.Thus, the proof is complete.

Improved Criteria Ensure That U
In what follows, we present a finer estimate for the ratio z(ρ(f))/z(f) that will allow us to improve the previous results and apply the new results when (20) holds.
Lemma 5. Assume that υ ∈ U + 3 .If there exist a f 1 ≥ f 0 and a constant β ≥ 0 such that Proof.Assume that υ ∈ U + 3 .Proceeding as in the proof of Theorem 2, we arrive at (16).From Lemma 4, we have that (4) holds for k = n − 3, that is, From ( 16) and ( 23), we obtain Now, we see that It follows from ( 21) that d df z(f) This completes the proof.
Proof.Assume the contrary, that υ ∈ U + 3 for f ≥ f 1 ≥ f 0 .Using Lemma 5, we obtain that (22) holds, and then As in the proof of Theorem 2, we arrive at (15).From ( 15) and ( 26), we get that By replacing ( 16) by (27), and completing the proof as in Theorem 2, we arrive at a contradiction with (25).This completes the proof.
In what follows, using the comparison principle with a second-order DDE, we infer a criterion which ensures that U + 3 = ∅.

Conclusions
The oscillatory properties of a class of NDEs were studied.To extend the evolution of the study of the second-order DDEs, we presented criteria with only two conditions that ensure the oscillation of the considered equation.Due to the importance of the conditions that exclude solutions in class U + 3 , we improve these conditions by finding a better estimate of the ratio z(ρ(f))/z(f).We establish more than one criterion for oscillation, and these new criteria are distinguished by taking into consideration the impact of delay functions, as well as not needing to define additional functions and conditions.Extension of our results to neutral equations of odd order will be interesting as one of the future issues.There are several studies concerned with the oscillatory properties of solutions of fractional DE, see [26][27][28][29].It is also interesting to extend the approach taken in this work to study the oscillatory properties of fractional DEs.

Author Contributions:
Conceptualization, O.M and F.M.; Formal analysis, O.M. and B.A.; Investigation, F.M. and D.A.; Methodology, F.M. and D.A.; writing-original draft preparation, D.A.; writing-review and editing, O.M. and B.A. All authors have read and agreed to the published version of the manuscript.Funding: Princess Nourah bint Abdulrahman University Researchers Supporting Project number PNURSP2022R216, Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.Institutional Review Board Statement: Not applicable.Informed Consent Statement: Not applicable.Data Availability Statement: Not applicable.