Neutral Emden–Fowler Differential Equation of Second Order: Oscillation Criteria of Coles Type

: In this work, we study the asymptotic and oscillatory behavior of solutions to the second-order general neutral Emden–Fowler differential equation ( a ( v ) η ( x ( v )) z ′ ( v )) ′ + q ( v ) F ( x ( g ( v ))) = 0, where v ≥ v 0 and the corresponding function z = x + p ( x ◦ h ) . Besides the importance of equations of the neutral type, studying the qualitative behavior of solutions to these equations is rich in analytical points and interesting issues. We begin by finding the monotonic features of positive solutions. The new properties contribute to obtaining new and improved relationships between x and z for use in studying oscillatory behavior. We present new conditions that exclude the existence of positive solutions to the examined equation, and then we establish oscillation criteria through the symmetry property between non-oscillatory solutions. We use the generalized Riccati substitution method, which enables us to apply the results to a larger area than the special cases of the considered equation. The new results essentially improve and extend previous results in the literature. We support this claim by applying the results to an example and comparing them with previous findings. Moreover, the reduction of our results to Euler’s differential equation introduces the well-known sharp oscillation criterion.


Introduction
Neutral differential equations (NDEs)-in which the largest derivative occurs on the solution with and without delay-are among the most important types of functional differential equations (FDEs).Their importance is due to two main reasons, the first of which is the many practical applications of this type in physics and engineering, and the second is that the study of this type of equation includes many exciting analytical problems.Many applications, such as population dynamics, control, automatic, fluid mixing, and vibrating masses attached to a flexible rod, involve models of NDEs (see Hill [1]).In particular, second-order NDEs are very useful in robotics for building bipedal robots and in biology for explaining the ability of the human body to balance itself [2].
The part of qualitative theory that is concerned with studying the asymptotic and oscillatory features of solutions of FDEs is called oscillation theory.The main objectives of oscillation theory are to characterize the relationship between oscillatory and other fundamental properties of solutions to various classes of differential equations, investigate zero distribution laws, estimate the number of zeros in each interval and the distance between neighboring zeros, and establish conditions for the existence of oscillatory (and non-oscillatory) solutions (see [3,4]).The property of symmetry plays an important role in oscillation theory.In most studies, researchers are interested in obtaining criteria that guarantee the absence of any positive solutions to the examined equations, while the absence of negative solutions is also ensured through symmetry between non-oscillatory solutions.Oscillation theory has been and still is a crucial numerical mathematical technique in many cutting-edge technologies and fields (see [5][6][7]).
The capacity to forecast whether the solutions of a system will oscillate or remain stable is important in many applications.Understanding oscillation, for example, is useful in mechanical and structural engineering when building systems that can resist periodic forces without experiencing resonance, which could lead to structural failure.For modeling phenomena in biological systems, such as brain activity or heart rhythms, where oscillating patterns are suggestive of normal function or diseased states, oscillation criteria are crucial.In electrical engineering, oscillation criteria are essential for designing oscillators and filters in circuits so that they effectively carry out their intended tasks.

General NDE of Second Order
In this work, we provide sufficient conditions to verify that all solutions of the NDE are oscillatory, where v ≥ v 0 and During the study, we assume the following conditions: where For the solution of (1), we mean a function 1) is called oscillatory if it has arbitrarily large zeros and, otherwise, it is said to be non-oscillatory.
In the remainder of the introduction, we review some relevant previous literature.The main results are then divided into three sections.In Section 2, we present some helpful lemmas.Then, the oscillation results section provides new criteria for testing the oscillation of the considered equation in different cases for the functions F and η.In the next section, we apply our results to some special cases of the equation under study.The Conclusions Section comes at the end of the results.

Related Literature
Next, we revise some related work that has contributed to studying the oscillatory behavior of solutions to FDEs.
One of the earliest seminal contributions to the study of oscillatory behavior was Coles' paper [8].He established a method for testing oscillatory behavior based on the definition of a class of positive and locally integrable functions.
In 2000, Manojlovic [9] investigated the oscillatory features of the equation where α ∈ R + and As a generalization and extension of the results in [10,11], Dzurina and Lacková [12] used the Riccati approach to establish criteria for the oscillation of the NDE (1) where Şahiner [13] studied the oscillatory behavior of solutions to (1) and considered the following three cases for F: Condition (S1) obviously implies (S2), but not the other way around.For instance, (S2) is true for F(x) = x 3 and η(x) = x 2 , even though these functions do not meet (S1).It is necessary for F to be differentiable in (S1) and (S2).It is evident that (S3) does not have this requirement.
Using the generalized Riccati transformation, Erbe et al. [14] presented conditions for the oscillation of the NDE where α ∈ R + .They obtained oscillation criteria of Kamenev type.In all the previously mentioned works, researchers always used the traditional relationship x > (1 − p)z that links x and z.This relationship is essential in studying the oscillatory behavior of NDEs, as it is the means to transform the NDE into a delay inequality.
There has been a significant surge in the study of FDE oscillation theory lately.Researchers have been interested in addressing some analytical problems in the oscillation theory of NDEs.Alsharidi and Muhib [15] studied the oscillatory features of an NDE with a mixed neutral term.Jadlovská et al. [16,17] provided sharp results for testing the oscillation of equations in the canonical case.For the non-canonical case, Muhib [18] presented more effective conditions for oscillation.Guo et al. [19] studied an NDE with positive and negative coefficients.
Moaaz et al. [20] discussed the oscillatory features of the NDE based on the improved relations , n ∈ Z + is even, and Hassan et al. [21] deduced the relation and studied the oscillation of the NDE Recently, Bohner et al. [22] established the relationship and presented improved oscillation criteria for (3), where Moaaz et al. [23] tested the oscillatory features, and presented the relation The study of second-order FDEs is of great importance from a theoretical standpoint, in addition to its various applications.One can observe the reflection of the advancement of the study of the oscillatory behavior of solutions of second-order equations on the study of the oscillation of even-order equations, see [24,25] (for quasi-linear equations), [26,27] (for super-linear equations), and [28] (for equations on Time Scales).
As an extension of the results in [20][21][22][23], we obtain new properties for positive solutions of Equation ( 1) and then employ them to obtain a new relationship between x and z.Using this relationship, we get oscillation conditions for solutions of (1), considering the cases (S1)-(S3).The approach used is an improved extension of the integral averaging technique.The use of improved relationships is directly reflected in the conditions of oscillation, as the new conditions are more efficient in the oscillation test for the studied equations.

Preliminary Results
For the sake of simplification, we provide the next notations: Ψ(v, v 0 ) represents the class of all bounded functions for v ≥ v 0 that are positive and locally integrable and S + represents the class of all eventually positive solutions of (1), and ).Moreover, we need the following conditions: ), we define an integral operator Φ.This operator is defined as in [30], where θ ∈ C([v 0 , ∞), R) and is specified in terms of H(v, s) and α(s).
for any m ≥ 0 is an integer.
Proof.Assume that x ∈ S + .Based on assumption (A2), there is a Hence, we conclude that a(v)η(x(v))z ′ (v) has a fixed sign, eventually.This is equiva- lent to saying that z ′ (v) > 0 or z ′ (v) < 0 for v ≥ v 2 , where v 2 is large enough.However, the case in which z ′ (v) < 0 contradicts condition (2), as shown below: which with (A3) gives Thus, Therefore, the proof is finished.
Using the results of Lemma 2, we conclude that Therefore, the proof ends.
Next, we derive a new relationship between x and z that helps improve the oscillation results.Lemma 4. If x ∈ S + , then Equation (1) can be expressed as the delay form for v ≥ v 1 and any integer m ≥ 0.
Proof.Assume that x ∈ S + .Based on assumption (A2), there is a From Lemma 1, we find that It follows from Lemmas 2 and 3 that Using the previous two inequalities in (10), we obtain Combining this relation with Equation (1), we arrive at Therefore, the proof ends.
Proof.Assume that x ∈ S + .Based on assumption (A2), there is a In addition, we observe that ( 6) and ( 7) hold for v ≥ v 1 .Now, we may verify that z ′ (v) > 0 for v ≥ v 1 .In fact, multiplying (7) by . This goes against (A4); hence, for v ≥ v 1 , z ′ (v) > 0 is required.Next, we find that (9) holds by repeating the procedures from the proofs of Lemmas 3 and 4. The proof is finished.
Lemma 6. Assume that (S3) is satisfied.If x ∈ S + , then (6) holds for some v 1 > v 0 .Additionally, Proof.Assume that x ∈ S + .As per the proofs of Lemmas 3 and 4, we then state that ( 6) and (9) hold for some v 1 ≥ v 0 .Thus, from ( 1) and (S3), we have Note that, It is evident from (13) that ( 12) is true.This brings the proof to a close.

Oscillation Results
In this section, we derive new criteria for the oscillation of solutions to Equation (1) in cases (S1)-(S3).
Therefore, we obtain that is, For v ≥ ℓ ≥ ℓ 1 , we obtain After integrating from ℓ to v and multiplying the above relation by ψ(v), we obtain Hence, Based on condition (16), there is a ℓ 1 ≥ ℓ such that for all v ≥ ℓ 1 .Then, and is obtained by condition (14).As shown, the function S(v) is positive, and we obtain We obtain by the Schwarz inequality.Note that where Using (19) and (20), for every v ≥ ℓ 1 and some ς, 0 ≤ ς < 1, we obtain Integrating ( 21) from ℓ 1 to v, we obtain in contrast to (15).Thus, we have proved the theorem in full.
Next, utilizing Philos-type integral average conditions, we show some new oscillation criteria for (1).
Theorem 2. Assume that (S1) and (N1) are satisfied.Equation ( 1) is oscillatory if there are Proof.Following the steps in the proof of Theorem 1, we observe that ( 18) is true for all v ≥ ℓ ≥ ℓ 0 .When we apply operator Φ(.; v, ℓ) to ( 18), we obtain Completing squares of ω in the above inequality, we obtain So, Thus, we have We derive a contradiction to condition (22) by dividing both sides of the inequality mentioned above and taking the lim sup in it as v → ∞.The proof of this theorem is therefore finished.Theorem 3. Assume that (S1) and (N1) are satisfied.Equation ( 1) is oscillatory if there are and where ρ 1 and ρ 2 satisfy Proof.As in the proof of Theorem 2, we find that ( 23) and ( 24) hold.By dividing ( 24) by H(v, ℓ), we may determine that by the use of ( 25) and ( 26).This implies that 1 From ( 23), we obtain Together with (25), this yields for v ≥ ℓ ≥ ℓ 0 and where c is a constant.Currently, we claim that In the case that ( 29) is not satisfied, there is a sequence Observe that (30) and can be used for a large enough n.For a sufficiently large n, this and (30) indicate that Hence, following the Schwarz inequality.Our result is based on ( 31) and (32).Contrary to (30), the right hand side of ( 33) is bounded by (26).Therefore, using (28), we obtain We obtain a contradiction with (27).The proof is now complete.

Oscillation Theorems for Case (S2)
Theorem 4. Assume that (S2), (N1), and (N2) are satisfied.Equation ( 1) is oscillatory if there are and where Proof.Assume the contrary, that x ∈ S + .From Lemma 5, we find that ( 6) and (9) hold for some ℓ > v 0 .Taking into account the function ω(v) as stated by (17), we arrive at ) and (N2), we can now determine that Hence, The remaining portions of the proof follow Theorem 1.
Likewise, as in Case (S1), we can obtain the following two theorems, so their proof has been omitted.Theorem 5. Assume that (S2), (N1), and (N2) are satisfied.Equation ( 1) is oscillatory if there are Proof.We begin with inequality (37) and continue as in the Theorem 2 proof.
For the sake of completeness, we declare an analogous theorem to Theorem 3 below.This may be obtained by following the same technique as in the proof of Theorem 3. Theorem 6. Assume that (S2), (N1), and (N2) are satisfied.Equation ( 1) is oscillatory if there are and where ρ 1 and ρ 2 satisfy 3.3.Oscillation Theorems for Case (S3) and where Proof.Assume the contrary, that x ∈ S + .From Lemma 6, we find that ( 6) and ( 12) hold for some ℓ > v 0 .The function ω 1 (v) is defined by for all v > ℓ 0 .After applying (12) and differentiating (45), we obtain ) The type of inequality ( 46) is the same as that of inequality (18).Therefore, we can finish the Theorem 7 proof using a similar process.
Likewise, as in Case (S1), we can obtain the following two theorems, so their proof has been omitted.

Special Cases
In this section, to facilitate the application of the results, we present some results by identifying special cases of functions Ψ and H. Corollary 1. Assume that (S1) and (N1) are satisfied.Equation ( 1) is oscillatory if there are and Proof.Taking the function Then, Equation (1) oscillates according to Theorem 1.
We define where n is an integer and n > 1.Furthermore, As such, we derive the following oscillation criteria as a result of Theorems 2 and 3.
Corollary 5. Assume that (S1) and (N1) are satisfied.Equation ( 1) is oscillatory if there exists a , and an integer n > 1 such that where ρ 1 and ρ 2 satisfy In the following example, we apply our results to a form of (1) and also compare the results with previous results in the literature.
On the other hand, using Corollary 3, we find ρ(v) = µv, and Therefore, Equation (53) oscillates if (55) holds.We note that Corollary 3 provides a more efficient criterion than Corollary 1. Figure 1 shows one of the numerical solutions to (53).Remark 2. Consider the special case of Equation (53) when δ = β = 0.In this case, condition (55) reduces to q 0 > 1 4µ p 0 . (56) In 2018, Grace et al. [32] investigated the oscillatory behavior of solutions of the neutral differential equation (NDE) a(v) z ′ (v) where α is a ratio of odd natural numbers.They used the Riccati method and the comparison method with a first-order equation.Their results are essentially an improvement on results in the literature that preceded their work.For testing the oscillation of (x(v) + p 0 x(λv)) ′′ + q 0 v 2 x(µv) = 0, the best results they obtained (in Example 3 [32]) were where ϵ = 1 1 + (1 − p 0 )q 0 µ .

Conclusions
The investigation of the oscillatory features of delay equations is rich in analytical issues.Therefore, there are always attempts to improve the relationships, inequalities, and techniques used in oscillation theory.In this study, we investigated the oscillatory properties of solutions to a class of neutral-type FDEs.We deduced a new relationship between x and z, which is one of the vital relations in the oscillation theory of NDEs.Next, we utilized the Riccati method to derive Coles-type oscillation criteria.Moreover, we introduced many conditions for the oscillation of the examined equation by considering cases (S1)-(S3).It is worth noting that our results are a direct improvement on previous relevant results.This is clear by setting M = 0, so we find that and the relation (11) reduces to the traditional relationship x > (1 − p)z.In future work, we look forward to using the same approach to study oscillatory behavior in the non-canonical case as well as the oscillation of higher-order equations.Moreover, we will work in the future to obtain oscillation standards with fewer restrictions than those imposed here.