Next Article in Journal
Response of Shock Isolators with Piecewise Linear Asymmetric Damping
Next Article in Special Issue
On a Hierarchy of Vector Derivative Nonlinear Schrödinger Equations
Previous Article in Journal
Symmetric Time-Variant IBLF-Based Tracking Control with Prescribed Performance for a Robot
Previous Article in Special Issue
Existence and General Decay of Solutions for a Weakly Coupled System of Viscoelastic Kirchhoff Plate and Wave Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Results on the Oscillation of Solutions of Third-Order Differential Equations with Multiple Delays

1
Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
2
Department of Mathematics, College of Science, Qassim University, P.O. Box 6644, Buraydah 51452, Saudi Arabia
3
Department of Mathematical Science, College of Science, Princess Nourah bint Abdulrahman University, P.O. Box 105862, Riyadh 11656, Saudi Arabia
*
Author to whom correspondence should be addressed.
Submission received: 15 September 2023 / Revised: 6 October 2023 / Accepted: 11 October 2023 / Published: 15 October 2023
(This article belongs to the Special Issue Differential Equations and Applied Mathematics)

Abstract

:
This study aims to examine the oscillatory behavior of third-order differential equations involving various delays within the context of functional differential equations of the neutral type. The oscillation criteria for the solutions of our equation have been obtained in this study to extend and supplement existing findings in the literature. In this study, a technique that relies on repeatedly improving monotonic properties was used in order to exclude positive solutions to the studied equation. Negative solutions are excluded based on the symmetry between the positive and negative solutions. Our results are important because they become sharper when applied to a Euler-type equation as compared to previous studies of the same equation. The significance of the findings was illustrated through the application of these findings to specific cases of the investigated equation.

1. Introduction

Delay differential equations (DDEs), a subclass of functional differential equations, take into account the system’s reliance on the past to produce predictions for the future that are more precise and effective. One of the most important roles that the concept of delay in systems is thought to play is modeling the length of time needed to complete certain unseen activities. The predator–prey model demonstrates a delay when the birth rate of the predator takes into account both present and past numbers of predators and prey. With the rapid development of communication technologies, transmitting measured signals to a remote control center has become much simpler. However, the biggest obstacle for engineers is the time that it takes for the signal to reach the controller after a measurement has been taken. In order to minimize the possibility of experimental instability and potential harm, this lag must be considered during the planning phase. Modeling such phenomena, as well as others, requires the use of DDEs (see [1,2,3,4,5]).
Neutral delay differential equations (NDDEs) are encountered in several kinds of phenomena, such as electric transmission line problems, which are utilized for interconnecting switching circuits in high-speed computers, the study of vibrating masses connected to elastic bars, the solution of variational problems involving time delays or in the theory of automatic control, and neuro-mechanical systems where inertia is a significant factor (see [6,7,8,9,10]). The reader is directed to consult the references [11,12,13,14,15] for comprehensive insights into the methodologies, techniques, and findings relating to the investigation of oscillatory behavior in third-order NDDEs. Furthermore, the aforementioned studies [16,17,18,19,20] primarily center their attention on the examination of DDEs with odd orders.
This study focuses on third-order NDDEs expressed in linear form with several delays
κ 2 ( ) κ 1 ( ) z ( ) + i = 1 Ω q i y τ i = 0 ,
where 0 , z ( ) : = y ( ) + p ( ) y ( σ ( ) ) , and Ω is a positive natural number. We suppose throughout this paper that the following hypotheses are fulfilled:
(A1)
κ 1 C 2 ( [ 0 , ) , 0 , ) , κ 2 C ( [ 0 , ) , 0 , ) and
0 1 κ 1 ( ρ ) d ρ = 0 1 κ 2 ( ρ ) d ρ = ;
(A2)
p , q i C ( [ 0 , ) , [ 0 , ) ) , q i 0 , q i ( ) does not vanish identically for each i = 1 , 2 , , Ω and 0 p ( ) p 0 < 1 ;
(A3)
τ i , σ C ( [ 0 , ) , R ) , τ i ( ) , σ ( ) , lim σ ( ) = , and lim τ i ( ) = , for each i = 1 , 2 , , Ω .
A function y C 3 ( [ y , ) , R ) , y 0 , is said to be a solution of (1), which has the property z , κ 1 z , and κ 2 κ 1 z belong to C 1 [ y , ) and satisfies (1) on [ y , ) . Furthermore, we consider only solutions y of (1) that satisfy
sup { | y ( ) | : L } > 0 , for all L y .
If a solution y is neither eventually positive nor eventually negative, then it is said to be oscillatory. Otherwise, it is said to be non-oscillatory. The equation itself is termed oscillatory if all of its solutions oscillate.
The previous studies on the oscillatory characteristics of neutral differential equations with odd orders primarily concentrated on establishing a suitable criterion for verifying whether the solutions exhibit oscillatory behavior or approach zero, as referenced in [21,22,23,24]. In the following, we provide some background details regarding the study of various classes of neutral differential equations.
In 2010, Baculíková and Džurina [11,25] investigated the asymptotic properties of the third-order NDDE
κ z γ + q f y τ = 0 .
They obtained conditions that test the convergence of all non-oscillatory solutions to zero. In [11], they used comparisons with first-order equations, while in [25], they obtained Hille and Nehari criteria. Thandapani and Li [26] found some fulfilling conditions that confirm that every solution of (3) either converges to zero or is oscillatory by using the Riccati transformation. In [27], Baculíková and Džurina examined the oscillation of the NDDE
κ z + q y τ = 0 .
They obtained results based on the comparison theorems, which allowed them to reduce the problem of the oscillation in a third-order equation coupled to a first-order equation.
As an improvement over and completion of previous studies, Džurina et al. [28] established conditions to ensure that all solutions of linear NDDE
κ 2 κ 1 z + q y τ = 0 ,
using a comparison with first-order delay equations. Moaaz et al. [29] investigated the oscillatory behavior of the NDDE
κ z α + q y α τ = 0 .
Utilizing the iterative technique, they established criteria of an iterative nature and found a criterion for the nonexistence of the so-called Kneser solutions.
Recently, Jadlovská et al. [30] studied the oscillation of the NDDE
κ 2 κ 1 y + q y τ = 0 .
Their results tested the convergence of all non-oscillatory solutions to zero. Their results are also sharp when applied to the Euler-type DDE, and they improved all previous results with regard to the criterion that tests the convergence of all non-oscillatory solutions to zero.
Our paper investigates the oscillatory properties of a third-order NDDE with multiple delays. The main motivation of this study is to extend the results of [28] to equations with multiple delays with respect to the convergence of non-oscillatory solutions to zero. Moreover, we create standards that guarantee the oscillation of all solutions of the studied equation by establishing a standard that excludes so-called Kneser solutions. Applying our results to a particular case of the considered equation supported the findings.

2. Preliminary Results

For convenience, we define the following:
τ : = min τ i ; i = 1 , 2 , , Ω , q i ( ) : = min q i ( ) , q i ( σ ( ) ,
τ ˜ : = max τ i ; i = 1 , 2 , , Ω ,
J 0 z = z , J 1 z = κ 1 z , J 2 z = κ 2 κ 1 z , J 3 z = κ 2 κ 1 z ,
M 1 ( ) : = 0 d ρ κ 1 ( ρ ) , M 2 ( ) : = 0 d ρ κ 2 ( ρ ) ,
and
M 12 ( ) : = 0 M 2 ( ρ ) κ 1 ( ρ ) d ρ .
Lemma 1. 
Ref. [28] Lemma 1—suppose that there is a constant l > 0 such that
lim κ 1 M 1 κ 2 M 2 = l .
Then,
M 12 ( ) ϵ 1 + l M 1 ( ) M 2 ( ) ,
eventually for all ϵ ( 0 , 1 ) .
To proceed with proving our results we need to define the following limits:
lim inf M 12 ( ) M 12 ( τ ( ) ) : = λ * ,
lim inf κ 2 ( ) M 2 ( ) M 12 ( τ ( ) ) i = 1 Ω q i ( ) ( 1 p ( τ i ( ) ) : = β * ,
and
lim inf M 2 β * ( ) M 12 ( ) 0 M 2 1 β * ( ρ ) κ 1 ( ρ ) d ρ : = k * , for β * ( 0 , 1 ) .

3. Main Results

In this section, we provide sufficient conditions to ensure the oscillation of all solutions of the studied equation. For the following results, we assume that λ * , β * , k * 0 , .
Lemma 2. 
Suppose that y is a positive solution of (1). Then, J 3 z ( ) 0 , and there are only two categories:
C l a s s ( 1 ) : z > 0 , J 1 z > 0 , J 2 z > 0 , C l a s s ( 2 ) : z > 0 , J 1 z < 0 , J 2 z > 0 .
Proof. 
The proof is straightforward; hence, we omit the details. □
Notation 1. 
By x 1 we mean that solution x with corresponding function has class (1) properties, while by x 2 we mean that the solution x with corresponding function has class (2) properties.

3.1. Class (1)

In this section, we present some characteristics of solutions that belong to class 1 . We also obtain criteria that rule out the existence of solutions with class 1 properties.
Remark 1. 
From the definition of M 12 , β * , and k * , the following can be concluded:
(R1)
For any β ( 0 , β * ) , there is a β 0 such that
κ 2 ( ) M 2 ( ) M 12 ( τ ( ) ) i = 1 Ω q i ( ) ( 1 p ( τ i ( ) ) β , for all β .
Lemma 3. 
Suppose that y 1 and β * > 0 . Then, eventually,
(a)
the functions J 2 z ( ) , J 1 z / M 2 ( ) , and z ( ) / M 12 ( ) converge to zero;
(b)
J 1 z / M 2 is decreasing;
(c)
z / M 12 is decreasing.
Proof. 
Suppose that y 1 . Since z ( ) y ( ) and z ( ) > 0 , we have y ( ) ( 1 p ( ) ) z ( ) , and so y ( τ i ) ( 1 p ( τ i ) z ( τ i ) . Thus, (1) becomes
J 3 z ( ) = i = 1 Ω q i y τ i i = 1 Ω q i ( ) ( 1 p ( τ i ( ) ) z τ i ( ) z τ ( ) i = 1 Ω q i ( ) ( 1 p ( τ i ( ) ) .
a Using the facts J 2 z ( ) > 0 and J 3 z ( ) 0 , it is obvious that J 2 z ( ) ϱ 0 as . Assume the contrary that ϱ 0 > 0 . Hence, it follows that J 2 z ( ) ϱ 0 > 0 . Therefore,
J 1 z ( ) 1 1 κ 2 ( ρ ) J 2 z ( ρ ) d ρ ,
and so
z ( ) 1 1 κ 1 ( ρ ) J 1 z ( ρ ) d ρ 1 1 κ 1 ( ρ ) 1 ρ 1 κ 2 ( u ) J 2 z ( u ) d u d ρ ϱ 0 1 1 κ 1 ( ρ ) 1 ρ 1 κ 2 ( u ) d u d ρ = ϱ 0 1 M 2 ρ κ 1 ( ρ ) d ρ > δ ϱ 0 M 12 .
for all δ 0 , 1 . Hence, (7) reduces to
J 3 z ( ) ϱ 0 M 12 ( τ ) i = 1 Ω q i ( ) 1 p ( τ i ( ) .
Integrating this inequality from 1 and using (6), we obtain for 2 = max 1 , β ,
J 2 z ( ) ϱ 0 1 M 12 ( τ ( ρ ) ) i = 1 Ω q i ( ρ ) ( 1 p ( τ i ( ρ ) ) ) d ρ β ϱ 0 1 1 κ 2 ( ρ ) M 2 ( ρ ) d ρ = β ϱ 0 ln M 2 ( ) M 2 ( 1 ) .
Then, lim J 2 z ( ) = , a contradiction. Consequently, ξ = 0 . The rest of the properties in (a) are proven directly by using L’Hopital’s principle.
b Using (8) and the fact that J 3 z ( ) 0 , we arrive at
J 1 z ( ) J 2 z ( ) 1 1 κ 2 ( ρ ) d ρ = M 2 ( ) J 2 z ( ) .
and thus,
J 1 z M 2 = J 2 z M 2 J 1 z M 2 2 κ 2 < 0 .
c Since J 1 z / M 2 tend to zero and the function is decreasing, we find
z ( ) 1 M 2 ( ρ ) κ 1 ( ρ ) J 1 z ( ρ ) M 2 ( ρ ) d ρ J 1 z ( ) M 2 ( ) 1 M 2 ( ρ ) κ 1 ( ρ ) d ρ > J 1 z ( ) M 2 ( ) M 12 ( ) .
Therefore,
z M 12 = J 1 z M 12 z M 2 M 12 2 κ 1 < 0 .
We have reached the end of the proof. □
Remark 2. 
From the definition of M 12 , β * , and k * , the following can be concluded:
(R2)
Assume that β * ( 0 , 1 ) , we can conclude that k * 1 . For any k ( 1 , ) , there is k 0 such that
M 2 β ( ) M 12 ( ) 0 M 2 1 β ( ρ ) κ 1 ( ρ ) d ρ k , k .
Lemma 4. 
Suppose that y 1 and β * > 0 . Then, eventually,
(C01)
J 1 z / M 2 1 β * ( ) is decreasing;
(C02)
J 1 z ( ) / M 2 1 β * ( ) converges to zero;
(C03)
z > k M 12 / M 2 J 1 z , and z / M 12 1 / k is decreasing.
Proof. 
Assume that y 1 . From Lemma 3, we have that (a), (b), and (c) hold.
C 01 We define
ω ( ) : = J 1 z ( ) M 2 ( ) J 2 z ( ) .
Thus, from (b), ω ( ) > 0 for all 1 . Hence, we obtain
ω ( ) = J 1 z ( ) M 2 ( ) J 2 z ( ) = M 2 ( ) J 3 z ( ) ,
which with (1) gives
ω ( ) = M 2 ( ) i = 1 Ω q i ( ) y ( τ i ( ) ) .
Since z ( ) y ( ) and z ( ) > 0 , we have y ( ) ( 1 p ( ) ) z ( ) . Thus,
ω ( ) M 2 ( ) i = 1 Ω q i ( ) ( 1 p ( τ i ( ) ) ) z τ i ( ) z τ ( ) M 2 ( ) i = 1 Ω q i ( ) ( 1 p ( τ i ( ) ) ) .
Using (6) and (c), we obtain
ω ( ) β z ( τ ( ) ) κ 2 ( ) M 12 ( τ ( ) ) β z ( ) κ 2 ( ) M 12 ( ) β J 1 z κ 2 ( ) M 2 ( ) .
for 2 1 . Integrating this inequality from 2 to , we have
ω ( ) β 2 J 1 z ( ρ ) κ 2 ( ρ ) M 2 ( ρ ) d ρ β J 1 z ( ) M 2 ( ) 2 1 κ 2 ( ρ ) d ρ > β J 1 z ( ) .
From the definition of ω , we obtain 1 β J 1 z ( ) > M 2 ( ) J 2 z ( ) , and so
J 1 z M 2 1 β = J 2 z M 2 1 β J 1 z M 2 M 3 2 β κ 2 < 0 , 3 .
C 02 From (14), we have β < 1 . From (13), we find that
ω ( ) β 2 J 1 z ( ρ ) κ 2 ( ) M 2 ( ρ ) d ρ β J 1 z ( ) M 2 1 β ( ) 2 1 κ 2 ( ) M 2 β ( ρ ) d ρ β 1 β J 1 z ( ) M 2 1 β ( ) M 2 1 β ( ) M 2 1 β ( 2 ) ϵ 1 β 1 β J 1 z ( )
for all ϵ 1 0 , 1 . Then, eventually,
ω ( ) β * + c 2 J 1 z ( ) ,
for c 2 > 0 , by choosing β β * / ( 1 + β * ) , β * . Hence,
1 β * J 1 z ( ) > 1 β * c 2 J 1 z ( ) > M 2 ( ) J 2 z ( ) .
Thus, J 1 z / M 2 1 β * c 2 is decreasing. Now, if we assume lim J 1 y / M 2 1 β * > 0 ,we find
J 1 z ( ) M 2 1 β * c 2 ( ) = J 1 z ( ) M 2 1 β * ( ) M 2 c 2 ( ) as ,
a contradiction. Therefore, lim J 1 y / M 2 1 β * = 0 .
C 03 From the previous facts we can conclude that
z ( ) = z ( 3 ) + 3 M 2 1 β * ( ρ ) κ 1 ( ρ ) J 1 z ( ρ ) M 2 1 β * ( ρ ) d ρ z ( 3 ) + J 1 z ( ) M 2 1 β * ( ) 3 M 2 1 β * ( ρ ) κ 1 ( ρ ) d ρ = z ( 3 ) + J 1 z ( ) M 2 1 β * ( ) 0 M 2 1 β * ( ρ ) κ 1 ( ρ ) d ρ J 1 z ( ) M 2 1 β * ( ) 0 3 M 2 1 β * ( ρ ) κ 1 ( ρ ) d ρ > J 1 z ( ) M 2 1 β * ( ) 0 M 2 1 β * ( ρ ) κ 1 ( ρ ) d ρ k M 12 ( ) M 2 ( ) J 1 z ( ) , 5 3 .
Therefore, eventually,
z M 12 1 / k < 0 .
We have reached the end of the proof. □
Lemma 5. 
If β * 1 , then 1 = .
Proof. 
Assume that y 1 . Since J 1 z / M 2 1 β * ( ) is decreasing and J 2 z > 0 , we obtain β * < 1 , a contradiction. Thus, 1 = . □
We can improve the previous properties by defining the following sequences:
β n = β * for n = 0 β 0 k n 1 λ * 1 1 k n 1 1 β n 1 for n = 1 , 2 , ,
and
k n = lim inf M 2 β n ( ) M 12 ( ) 0 M 2 1 β n ( ρ ) κ 1 ( ρ ) d ρ , for n = 0 , 1 , 2 , .
Remark 3. 
If β i < 1 and k i [ 1 , ) for i = 0 , 1 , , n , then β n + 1 can be clearly determined. In this case, the following inequality holds:
β 1 = β 0 k 0 1 β 0 λ * 1 1 k 0 > β 0 .
Thus,
k 1 = lim inf M 2 β 0 ( ) 0 M 2 1 β 0 ( ρ ) κ 1 ( ρ ) d ρ M 12 ( ) = lim inf M 2 β 1 ( ) 0 M 2 1 β 0 β 1 β 0 ( ρ ) κ 1 ( ρ ) d ρ M 12 ( ) lim inf M 2 β 0 ( ) 0 M 2 1 β 0 ( ρ ) κ 1 ( ρ ) d ρ M 12 ( ) = k ,
Applying induction to n, it is straightforward to demonstrate that
β n + 1 β n = l n > 1 ,
where
0 : = k 0 λ * 1 1 k 0 1 β 0 , n : = k n 1 β n 1 λ * 1 1 k n 1 1 k n 1 β n k n 1 , n N ,
with
k n k n 1 .
Lemma 6. 
Ref. [28] Lemma 2—suppose that (4) holds, β * > 0 , and β i < 1 , i = 0 , 1 , , n . We can conclude that
k n ( ) β n l 1 + l + 1 > 1 , n N 0 .
Lemma 7. 
Suppose that y 1 and β * > 0 . Then, eventually,
(Cn1)
J 1 z / M 2 1 β n is decreasing;
(Cn2)
lim J 1 z ( ) / M 2 1 β n ( ) converges to zero;
(Cn3)
z > ε n k n M 12 / M 2 J 1 z and y / M 12 1 / ε n k n is decreasing for any ε n ( 0 , 1 ) .
Proof. 
Assume that y 1 . We will employ an induction argument on n. For n = 0 , the conclusion directly follows from Lemma 4 with ε 0 = k / k * . Next, assuming that C n 1 C n 3 hold when n 1 for n 1 , we need to demonstrate that these conditions also hold for n + 1 .
Based on Lemma 4, the proof is exactly similar to the proof of Lemma 5 in [28]; therefore, it was omitted. □
Corollary 1. 
If β i < 1 for i = 0 , 1 , , n 1 , and β n 1 , then 1 = .
From the previous results and taking into account (18), the sequence β n has the limit
lim β n n = β ζ : = β * k ζ λ * 1 1 / k ζ 1 β ζ ( 0 , 1 ) ,
where
k ζ = lim inf M 2 β ζ ( ) M 3 ( ) 0 M 2 1 β ζ ( ρ ) κ 1 ( ρ ) d ρ .
Theorem 1. 
If (20) does not possess a root on ( 0 , 1 ) , then 1 = .
Corollary 2. 
If
β * > α ,
then 1 = , where
α : = max β ζ ( 1 β ζ ) k ζ 1 λ * 1 / k ζ 1 : 0 < β ζ < 1 .

3.2. Class (2)

Lemma 8. 
Assume y 2 and
0 1 κ 1 ( s ) s 1 κ 2 ( ρ ) ρ i = 1 Ω q i ( u ) d u d ρ d s = .
Then, lim z ( ) = 0 .
Proof. 
Assume that y 2 . Since z ( ) > 0 and z ( ) < 0 , we have z ( ) ν 0 as , where ν 0 0 . Assume that ν 0 > 0 ; then, we have for all ϵ > 0 , ν 0 < z ( ) < ν 0 + ϵ , eventually. By choosing 0 < ϵ < ν 0 ( 1 p 0 ) p 0 , it is easy to verify that
y ( ) = z ( ) p ( ) y ( σ ( ) ) > ν 0 p 0 z ( σ ( ) ) > ν 0 p 0 ( ν 0 + ϵ ) > h ( ν 0 + ϵ ) > h z ,
where h = ν 0 p 0 ( ν 0 + ϵ ) ν 0 + ϵ > 0 . Then, (1) becomes
J 3 z ( ) = i = 1 Ω q i ( ) y ( τ i ( ) ) h i = 1 Ω q ( i ) z ( τ ) h ν 0 i = 1 Ω q ( i ) .
Integrating from to , we obtain
J 2 z ( ) ν 0 h i = 1 Ω q i ( ρ ) h ν 0 d ρ ,
In other words,
J 1 z ( ) ν 0 h κ 2 ( ) i = 1 Ω q i ( ρ ) d ρ .
Integration (23) from to gives
z ( ) ν 0 h κ 1 ( ) 1 κ 2 ( u ) u i = 1 Ω q i ( ρ ) d ρ d u ,
and hence,
z ( ) z ( 2 ) ν 0 h 2 1 κ 1 ( x ) x 1 κ 2 ( u ) u i = 1 Ω q i ( ρ ) d ρ d u d x as ,
which contradicts the positivity of z. Then, the proof of this lemma is complete. □
In the following theorem, we establish certain conditions that guarantee the absence of Kneser solutions, which are solutions whose corresponding function satisfies the properties in class (2). In the following, we need the conditions
τ i ( σ ) = σ ( τ i ) , and σ σ 0 > 0 .
Theorem 2. 
Suppose that there is a function ζ C 0 , , 0 , satisfying τ ˜ < ζ and σ 1 ζ < . If the DDE
w + σ 0 σ 0 + p 0 K ζ , τ ˜ w σ 1 ζ i = 1 Ω q ˜ i = 0 ,
is oscillatory, then 2 = , where
K ς , ϱ : = ϱ ς 1 κ 1 s s ς 1 κ 2 u d u d s .
Proof. 
Assume that y 2 . This implies that
z > 0 , J 1 z < 0 , and J 2 z > 0 .
From (1), we see that
0 p 0 σ κ 2 σ κ 1 σ z σ + p 0 i = 1 Ω q i σ x ( τ i ( σ ) ) p 0 σ 0 J 3 z σ + p 0 i = 1 Ω q i σ x ( τ i ( σ ) ) = p 0 σ 0 J 3 z σ + p 0 i = 1 Ω q i σ x σ τ i .
Combining (1) and (26), we obtain
0 J 3 z + p 0 σ 0 J 3 z σ + i = 1 m q i x ( τ i ( ) ) + p 0 i = 1 Ω q i σ x σ τ i J 3 z + p 0 σ 0 J 3 z σ + i = 1 Ω q ˜ i t x ( τ i ( ) ) + p 0 x σ τ i .
From definition of z, we have
z τ i = x τ i + p τ i x σ τ i x τ i + p 0 x σ τ i .
By using the latter inequality in (27), we obtain
0 J 3 z + p 0 σ 0 J 3 z σ + i = 1 Ω q ˜ i t z ( τ i ( ) ) .
Since z is decreasing, then
0 J 3 z + p 0 σ 0 J 3 z σ + z τ ˜ i = 1 Ω q ˜ i
That is
0 J 2 z + p 0 τ 0 J 2 z σ + z τ ˜ i = 1 Ω q ˜ i .
On the other hand, it follows from the monotonicity of J 2 z that
J 1 z ϱ J 1 z ς J 1 z ϱ = ϱ ς J 1 z s d s = ϱ ς J 2 z s κ 2 s d s J 2 z ς ϱ ς 1 κ 2 s d s
Integrating (29) from ϱ to ς , we have
z ϱ J 2 z ς ϱ ς 1 κ 1 s s ς 1 κ 2 u d u d s .
Thus, we have
z τ ˜ J 2 z ζ K ζ , τ ˜ ,
which, by virtue of (28), yields that
J 2 z + p 0 σ 0 J 2 z σ + K ζ , τ ˜ J 2 z ζ i = 1 Ω q ˜ i 0 .
Now, set
w t = J 2 z + p 0 σ 0 J 2 z σ > 0 .
From the fact that J 2 z is non-increasing, we have
w t J 2 z σ 1 + p 0 σ 0 ,
or equivalently,
J 2 z ζ σ 0 σ 0 + p 0 w σ 1 ζ .
Using (32) in (31), we see that w is a positive solution of the differential inequality
w + σ 0 σ 0 + p 0 K ζ , τ ˜ w σ 1 ζ i = 1 Ω q ˜ i 0 .
In view of [31] Theorem 1, we have that (24) also has a positive solution, a contradiction. Thus, the proof is complete. □
Corollary 3. 
Suppose that there is a function ζ C 0 , , 0 , satisfying τ ˜ < ζ and σ 1 ζ < . If
lim inf t σ 1 ζ t t K ζ s , τ s i = 1 Ω q ˜ i s d s > σ 0 + p 0 σ 0 e ,
then 2 = .
Proof. 
The results in [32] guarantee the oscillation of Equation (24) under condition (33). □

3.3. Oscillatory Theorems and Examples

We obtain the criteria in the following theorems by directly combining the results in the previous two subsections. Assuming that the solution is positive means that it belongs to one of two categories: 2 or 2 . Therefore, when it is confirmed that categories 2 and 2 are empty, this means that there are no positive solutions, and accordingly, all solutions are oscillatory (this is based on the principle of symmetry between positive and negative solutions).
Theorem 3. 
Suppose β * > 1 , and (22) holds. Then, every solution of (1) either converges to zero or is oscillatory.
Theorem 4. 
Suppose (21) and (22) hold. Then, every solution of (1) either converges to zero or is oscillatory.
Theorem 5. 
Suppose β * > 1 , and (33) holds. Then, every solution of (1) is oscillatory.
Theorem 6. 
Suppose (21) and (33) hold. Then, every solution of (1) is oscillatory.
The following example demonstrates the significance of the results obtained.
Example 1. 
Consider
y ( ) + p 0 y σ 0 + i = 1 Ω q 0 3 y ( τ i ) = 0 ,
where 0 p 0 < 1 , and τ i , σ 0 ( 0 , 1 ) . Clearly: κ 1 ( ) = κ 2 ( ) = 1 , σ = ρ 0 , τ = τ 0 = min τ i , i = 1 , 2 , , Ω , p = p 0 , q i = q 0 / 3 , and
M 1 ( ) , M 2 ( ) , M 12 ( ) 2 / 2 .
Then, we can compute the value of β * as follows:
β * = lim inf M 2 ( ) κ 2 M 12 ( τ ) i = 1 Ω q i 1 p τ i = lim inf τ 0 2 2 2 Ω q 0 3 1 p 0 = 1 2 Ω τ 0 2 1 p 0 q 0 ,
For β * 1 , we have
q 0 > 2 Ω τ 2 ( 1 p 0 ) .
Moreover
0 1 κ 1 ( v ) 1 κ 2 ( u ) ρ φ ( ρ ) d ρ d u d v = 1 ρ Ω q 0 ρ 3 ( 1 p 0 ) d ρ d u d v = .
Thus, the assumption of Theorem 3 is satisfied, and then, every solution of (34) either converges to zero or is oscillatory.
Example 2. 
Consider the third-order neutral delay differential equation
e y ( ) + p 0 y ( σ ( ) + i = 1 Ω q 0 y τ i = 0 , > 1 .
where q 0 > 0 . It is easy to verify that
κ 1 ( ) = 1 , κ 2 ( ) = e , σ = σ 0 , τ = τ 0 = min τ i , i = 1 , 2 , , Ω , p = p 0 , q i = q 0 , and
M 1 ( ) , M 2 ( ) e , M 12 ( ) e .
Then
λ * = lim t inf M 12 ( ) M 12 ( τ ( ) ) = lim inf e 1 τ = ,
and
β * = lim inf M 2 ( ) κ 2 M 12 ( τ ) i = 1 Ω q i 1 p τ i = lim inf e τ Ω q 0 1 p 0 > 0 ,
where 0 p 0 < 1 . Moreover,
0 1 κ 2 ( u ) u i = 1 Ω q i ( ρ ) d ρ d u = 1 e u u Ω q 0 d ρ d u = 1 Ω q 0 e u e d u = .
Hence, if (36) holds, all assumptions of Theorem 3 are satisfied, and then every solution of (35) either converges to zero or is oscillatory.
Remark 4. 
Consider the differential equation
y + 0.5 y 6 + q 0 3 x 0.5 + x 0.6 + x 0.7 = 0 ,
where Ω = 3 , p 0 = 0.5 , σ = / 6 and τ = min 0.5 , 0.6 , 0.7 = 0.5 . Then, every solution of (37) either converges to zero or is oscillatory.
q 0 > 2 3 0.5 2 1 1 2 = 5.3 .

4. Conclusions

This study focuses on the oscillatory characteristics of solutions to the third-order neutral equation with several delays. Although there have been numerous studies related to this subject, we have discovered enough evidence in these studies to assert that any non-oscillatory solution will lead to zero. In this study, we introduce new standards which guarantee that all solutions to Equation (1) are oscillatory. Our results expand and improve upon those found in the literature [28]. For a certain type of general third-order delay differential equation, we propose new oscillation criteria in the event that the functions κ i are of the same kind, using an iterative technique in Theorems 5 and 6. In a particular case, a single condition ensures that Equation (1) oscillates. It is noteworthy that our criteria are relevant even when τ ( ) = , as they do not require τ ( ) to be a non-decreasing function.
In future research endeavors within this particular domain, we are enthusiastic about the potential to expand our investigation to include quasi-linear third-order neutral differential equations of the form:
κ 2 ( ) κ 1 ( ) z ( ) α + i = 1 Ω q i y α τ i = 0 .

Author Contributions

Conceptualization N.O., O.M., G.A. and E.M.E.; Methodology, N.O., O.M., G.A. and E.M.E.; Investigation, N.O., O.M., G.A.; Writing– original draft, N.O.; Writing – review and editing, O.M., G.A. and E.M.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R45), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Acknowledgments

Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R45), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Conflicts of Interest

This work does not have any conflict of interest.

References

  1. Hale, J.K. Functional differential equations. In Analytic Theory of Differential Equations; Springer: Berlin/Heidelberg, Germany, 1971; pp. 9–22. [Google Scholar]
  2. Gopalsamy, K. Stability and Oscillations in Delay Differential Equations of Population Dynamics; Springer: Dordrecht, The Netherlands, 1992. [Google Scholar]
  3. Agarwal, R.P.; Bohner, M.; Li, W.-T. Nonoscillation and Oscillation Theory for Functional Differential Equations; CRC Press: Boca Raton, FL, USA, 2004; Volume 267. [Google Scholar]
  4. Cimen, E.; Ekinci, Y. Numerical method for a neutral delay differential problem. Int. Math. Comput. Sci. 2017, 1, 1–11. [Google Scholar]
  5. Cimen, E.; Uncu, S. On the solution of the delay differential equation via Laplace transform. Commun. Math. Appl. 2020, 11, 379. [Google Scholar]
  6. Alqudah, M.A.; Ravichandran, C.; Abdeljawad, T.; Valliammal, N. New results on Caputo fractional-order neutral differential inclusions without compactness. Adv. Differ. Equ. 2019, 2019, 528. [Google Scholar] [CrossRef]
  7. Bainov, D.D.; Mishev, D.P. Oscillation Theory for Neutral Differential Equations with Delay; CRC Press: Boca Raton, FL, USA, 1991. [Google Scholar]
  8. Logeswari, K.; Ravichandran, C. A new exploration on existence of fractional neutral integro-differential equations in the concept of Atangana–Baleanu derivative. Phys. Stat. Mech. Its Appl. 2020, 544, 123454. [Google Scholar] [CrossRef]
  9. Ravichandran, C.; Valliammal, N.; Nieto, J.J. New results on exact controllability of a class of fractional neutral integro-differential systems with state-dependent delay in Banach spaces. J. Franklin Inst. 2019, 356, 1535–1565. [Google Scholar] [CrossRef]
  10. Valliammal, N.; Ravichandran, C. Results on fractional neutral integro-di¤erential systems with state-dependent delay in Banach spaces. Nonlinear Stud. 2018, 25, 159–171. [Google Scholar]
  11. Baculikova, B.; Dzurina, J. On the asymptotic behavior of a class of third order nonlinear neutral differential equations. Cent. Eur. J. Math. 2010, 8, 1091–1103. [Google Scholar] [CrossRef]
  12. Candan, T. Asymptotic properties of solutions of third-order nonlinear neutral dynamic equations. Adv. Differ. Equ. 2014, 2014, 35. [Google Scholar] [CrossRef]
  13. Chatzarakis, G.E.; Dzurina, J.; Jadlovska, I. Oscillatory and asymptotic properties of third-order quasilinear delay differential equations. J. Inequal. Appl. 2019, 2019, 23. [Google Scholar] [CrossRef]
  14. Chatzarakis, G.E.; Grace, S.R.; Jadlovska, I. Oscillation criteria for third-order delay differential equations. Adv. Differ. Equ. 2017, 2017, 330. [Google Scholar] [CrossRef]
  15. Jadlovská, I.; Chatzarakis, G.E.; Džurina, J.; Grace, S.R. On sharp oscillation criteria for general third-order delay differential equations. Mathematics 2021, 9, 1675. [Google Scholar] [CrossRef]
  16. Moaaz, O.; Awrejcewicz, J.; Muhib, A. Establishing new criteria for oscillation of odd-order nonlinear differential equations. Mathematics 2020, 8, 937. [Google Scholar] [CrossRef]
  17. Baculikova, B.; Dzurina, J. On the oscillation of odd order advanced differential equations. Bound. Value Probl. 2014, 2014, 214. [Google Scholar] [CrossRef]
  18. Moaaz, O.; Dassios, I.; Muhsin, W.; Muhib, A. Oscillation theory for non-linear neutral delay differential equations of third order. Appl. Sci. 2020, 10, 4855. [Google Scholar] [CrossRef]
  19. Muhib, A.; Abdeljawad, T.; Moaaz, O.; Elabbasy, E.M. Oscillatory properties of odd-order delay differential equations with distribution deviating arguments. Appl. Sci. 2020, 10, 5952. [Google Scholar] [CrossRef]
  20. Masood, F.; Moaaz, O.; Askar, S.S.; Alshamrani, A. New Conditions for Testing the Asymptotic Behavior of Solutions of Odd-Order Neutral Differential Equations with Multiple Delays. Axioms 2023, 12, 658. [Google Scholar] [CrossRef]
  21. Jiang, Y.; Li, T. Asymptotic behavior of a third-order nonlinear neutral delay differential equation. J. Inequal. Appl. 2014, 2014, 512. [Google Scholar] [CrossRef]
  22. Graef, J.; Tunc, E.; Grace, S. Oscillatory and asymptotic behavior of a third-order nonlinear neutral differential equation. Opuscula Math. 2017, 37, 839–852. [Google Scholar] [CrossRef]
  23. Thandapani, E.; Tamilvanan, S.; Jambulingam, E.; Tech, V. Oscillation of third order half linear neutral delay differential equations. Int. J. Pure Appl. Math. 2012, 77, 359–368. [Google Scholar]
  24. Tunc, E. Oscillatory and asymptotic behavior of third-order neutral differential equations with distributed deviating arguments. Electron. J. Differ. Equ. 2017, 16, 1–12. [Google Scholar]
  25. Baculikova, B.; Dzurina, J. Oscillation of third-order neutral differential equations. Math. Comput. Model. 2010, 52, 215–226. [Google Scholar] [CrossRef]
  26. Thandapani, E.; Li, T. On the oscillation of third-order quasi-linear neutral functional differential equations. Arch. Math. 2011, 47, 181–199. [Google Scholar]
  27. Baculikova, B.; Dzurina, J. Oscillation theorems for third order neutral differential equations. Carpathisn J. Math. 2012, 28, 199–206. [Google Scholar] [CrossRef]
  28. Dzurina, J.; Grace, S.R.; Jadlovska, I. On nonexistence of Kneser solutions of third-order neutral delay differentia equations. Appl. Math. Lett. 2019, 88, 193–200. [Google Scholar] [CrossRef]
  29. Moaaz, O.; Mahmoud, E.E.; Alharbi, W.R. Third-Order Neutral Delay Differential Equations: New Iterative Criteria for Oscillation. J. Funct. Spaces 2020, 2020, 8. [Google Scholar] [CrossRef]
  30. Baculikova, B.; Dzurina, J. Oscillation of third-order nonlinear differential equations. Appl. Math. Lett. 2011, 24, 466–470. [Google Scholar] [CrossRef]
  31. Philos, C. On the existence of nonoscillatory solutions tending to zero at for differential equations with positive delays. Arch. Math. 1981, 36, 168–178. [Google Scholar] [CrossRef]
  32. Kitamura, Y.; Kusano, T. Oscillation of first-order nonlinear differential equations with deviating arguments. Proc. Amer. Math. Soc. 1980, 78, 64–68. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Omar, N.; Moaaz, O.; AlNemer, G.; Elabbasy, E.M. New Results on the Oscillation of Solutions of Third-Order Differential Equations with Multiple Delays. Symmetry 2023, 15, 1920. https://0-doi-org.brum.beds.ac.uk/10.3390/sym15101920

AMA Style

Omar N, Moaaz O, AlNemer G, Elabbasy EM. New Results on the Oscillation of Solutions of Third-Order Differential Equations with Multiple Delays. Symmetry. 2023; 15(10):1920. https://0-doi-org.brum.beds.ac.uk/10.3390/sym15101920

Chicago/Turabian Style

Omar, Najiyah, Osama Moaaz, Ghada AlNemer, and Elmetwally M. Elabbasy. 2023. "New Results on the Oscillation of Solutions of Third-Order Differential Equations with Multiple Delays" Symmetry 15, no. 10: 1920. https://0-doi-org.brum.beds.ac.uk/10.3390/sym15101920

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop