Next Article in Journal
On the Effects of Tokamak Plasma Edge Symmetries on Turbulence Relaxation
Previous Article in Journal
The Problem of Moments: A Bunch of Classical Results with Some Novelties
Previous Article in Special Issue
Fourth-Order Emden–Fowler Neutral Differential Equations: Investigating Some Qualitative Properties of Solutions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimizing the Monotonic Properties of Fourth-Order Neutral Differential Equations and Their Applications

1
Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
2
Section of Mathematics, International Telematic University Uninettuno, CorsoVittorio Emanuele II, 39, 00186 Rome, Italy
3
Department of Statistics and Operations Research, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
*
Authors to whom correspondence should be addressed.
Symmetry 2023, 15(9), 1744; https://doi.org/10.3390/sym15091744
Submission received: 26 May 2023 / Revised: 30 August 2023 / Accepted: 8 September 2023 / Published: 11 September 2023

Abstract

:
We investigate the oscillation of the fourth-order differential equation for a class of functional differential equations of the neutral type. We obtain a new single-oscillation criterion for the oscillation of all the solutions of our equation. We establish new monotonic properties for some cases of positive solutions of the studied equation. Moreover, we improve these properties by using an iterative method. This development of monotonic properties contributes to obtaining new and more efficient criteria for verifying the oscillation of the equation. The results obtained extend and improve previous findings in the literature by using an Euler-type equation as an example. The importance of the results was clarified by applying them to some special cases of the studied equation. The fourth-order delay differential equations have great practical importance due to their wide applications in civil, mechanical, and aeronautical engineering. Research on this type of equation is still ongoing due to its remarkable importance in many fields.

1. Introduction

Differential equations have been a significant area of pure and applied mathematics since their establishment in the middle of the 17th century. Despite their extensive study in the past, it remains an important field for research with the arrival of new connections with other branches of mathematics, the fruitful interaction with applied fields, the interesting reformulation of fundamental issues and theories in various eras, the new perspectives in the twentieth century, and so on. Ordinary differential equations (ODEs) have several applications in mathematics and other fields, but when they are used to explain certain phenomena, including natural phenomena, we find that they contain delay times in their modeling, which leads to the so-called delay differential equations (DDEs). DDEs are a type of differential equation that takes into account time delays in the dynamics of the system. This indicates that the delay differential equation can directly represent any event that happened in the past, which gives it the ability to capture and analyze the behavior of systems where time delays play a critical role. Therefore, it is easy to see how these equations are utilized in physics, engineering, biology, and other sciences (see references [1,2]). A subtype of delay differential equations is known as neutral delay differential equations (NDDEs), where the highest-order derivative of the unknown function appears on the solution both with and without delay, and the development of NDDEs involves past values of the time and state variables. The delay differential equation solution requires information about the state at a certain time in the past in addition to the current state. There are numerous applications for neutral delay differential equations (NDDEs) in science and engineering. They are employed in the modeling of systems with delayed feedback, including control systems, neural networks, chemical reactions, and populations, as highlighted in references [3,4]. One of the fundamental goals of oscillation theory is to find sufficient conditions to ensure that all differential equation solutions oscillate. The first monograph that dealt with oscillation theory was that of Ladas et al. [5], which covered the results until 1984. There has been a lot of research done in the last few years on the oscillation and the oscillatory properties of differential equations (see references [6,7,8,9,10]). In recent years, there have been numerous studies on the oscillation and non-oscillation of solutions to various kinds of neutral functional differential equations (see references [11,12]). Numerous authors have examined the oscillations of fourth-order differential equations, and a number of methods for generating oscillatory criteria for these equations [13,14].
In this paper, we pay particular attention to the oscillatory behavior of solutions to the fourth-order neutral differential equation
z 4 s + q s x τ s = 0 ,
where s s 0 , z s = x s + p s x σ s is called the corresponding function of the solution x. We will assume the following conditions:
(H 1 )
p , q C s 0 , , 0 , , 0 p s p 0 < 1 ;
(H 2 )
τ , σ C s 0 , , R , τ ( s ) s , σ ( s ) s and lim s τ ( s ) = , lim s σ ( s ) = .
Via a solution of (1), we mean a function x C s x , , R for s x s 0 , which has the property z C 4 s x , , and satisfies (1) on s x , . We only take into account the solutions x of (1) that satisfy Sup x s : s T > 0 for all T s x oscillatory.
Definition 1.
If solution x for (1) is ultimately positive or negative, it is said to be non-oscillatory; if not, it is said to be oscillatory. If all of an equation’s solutions oscillate, the equation itself is said to be oscillatory.
One of the most important motivations for conducting this research is the importance of neutral differential equations, which have many uses in technology and natural science. They are often employed, for instance, in the study of distributed networks with lossless transmission lines (see reference [15]), therefore, their qualitative characteristics are crucial. Besides the importance of fourth-order differential equations in mathematical representations of several physical, biological, and chemical phenomena, fourth-order differential equations are frequently encountered. Their applications include, for example, elasticity issues, structure distortion, or soil settlement (see reference [16]). Complementary to the motives behind this paper is the fact that one of the conditions for oscillation is to find a condition in the form of a Kneser-type oscillation. The Kneser oscillation theorem states that a second-order linear differential equation of the form
x s + q s x s = 0 ,
is oscillatory if
lim   inf   s s 2 q s > 1 4 ,
where one finds that the above condition ensures the oscillatory behavior of all solutions while there is a positive solution in the case
lim sup s s 2 q s < 1 4 .
Thus, such conditions are more accurate and effective for the oscillation test. Therefore, the aim of this study was to extend the results obtained in the second- and fourth-order delay equations to neutral. Moreover, there are a number of related results that inspired our study in particular; Chatzarakis et al. [17] analyzed the oscillation behavior of the following fourth-order differential equations
r s x s + p s x τ s α + a b q s , ν f x σ s , ν d ν = 0 ,
under the canonical case s 0 r 1 / α ζ d ζ = .
Li et al. [11] studied the oscillatory behavior of the fourth-order nonlinear differential equation:
r s z s ( 4 ) + q s x σ s = 0 .
Bazighifan et al. [18] study the oscillatory properties of solutions of the following equation:
r s z s β + i = 1 j q i s x k τ i s = 0 ,
for s s 0 , under the canonical case s 0 r 1 / α ζ d ζ = .
In this work, we investigate the oscillatory behavior of solutions of fourth-order differential equations with neutral-delay arguments. We establish new monotonic properties for some cases of positive solutions of the studied equation and improve these properties by using an iterative method. This development of monotonic properties contributes to obtaining new and more efficient criteria for verifying the oscillation of the equation. The solutions of any equation are classified as positive, negative, and oscillatory. Most of the techniques used in studying oscillation to find oscillation standards are based on excluding positive and negative solutions. In this paper, we are interested in finding conditions that exclude positive solutions only, and this is based on the fact that every negative value of a positive solution to the studied equation is also considered a solution, or what is called symmetry between positive and negative solutions. As usual, Euler-type differential equations are used to highlight the improvement over the previous results from the literature. We will organize our paper as follows. In Section 2.1, we introduce the essential notations and the base of the method established in the sequel. In Section 2.2, we introduce the we introduce a number of lemmas that iteratively enhance the monotonic properties of the positive solutions. In Section 2.3, we present the main results and our main oscillations, that is, a single-oscillation criterion for (1) based on a series of lemmas. In the end, we highlight the importance of our results by comparing them with previous results in the literature.
Lemma 1
([19]). Let F C m s 0 , , R + . If F m is eventually of one sign for all large s, say, s 1 s 0 , then there exists a s x s 0 and an integer l, 0 l m , with m + l even for F m s 0 , or m + l odd for F m s 0 such that
l 0 implies that F k s > 0 for s s x , k = 0 , 1 , l 1 ,
and l m 1 implies that 1 l + k F k s > 0 for s s x , k = l , l + 1 , m 1 .

2. Main Results

In this section, we will establish some important lemmas that we will use in the proof to illustrate the main results of the research.
Notation 1.
Firstly, we will display the important notation used in this paper. Our results are dependent on the necessity of positive β * stated by
β * = lim inf s τ 3 s s q s 1 p τ s 3 ! ,
also, let us define
γ * = lim inf s τ s s 3 q s 1 p τ s 3 ! , δ * = lim inf s s τ s ,
where β * , γ * are positive because of (H 1 ) and (H 2 ). We will use, in our proof, the statement that there is a sufficiently large s 1 s 0 , such that
τ 3 s s q s 1 p τ s 3 ! β , τ s s 3 q s 1 p τ s 3 ! γ and s τ s δ ,
on s 1 , , where for arbitrary but fixed
β 0 , β * , γ 0 , γ * and δ 1 , δ * ,
for δ * > 1 , and δ = δ * for δ * = 1 .
In the following lemma, we classify the signs of the derivatives of non-oscillatory solutions to study the oscillatory features of solutions.
Lemma 2.
Assume that x is a positive solution of (1), then there are eventually only two possible cases for z
C a s e 1 z s > 0 , z s > 0 , z s > 0 , z s > 0 , z 4 s < 0 , C a s e 2 z s > 0 , z s > 0 , z s < 0 , z s > 0 .
Proof. 
Let x be a positive solution of (1), we get z 4 s 0 from (1). By using Lemma 1, we obtain case 1 , case 2 , and their derivatives. □
Notation 2.
We will eventually refer to the class of positive solutions whose corresponding function to Case 1 by ρ 1 , and whose corresponding function to Case 2 by ρ 2 .
Lemma 3.
Suppose that x is a positive solution of (1), then
z 4 s + q s 1 p τ s z τ s 0 .
Proof. 
Suppose that x is a positive solution of (1), it follows that there exists s 1 s 0 such that x s > 0 , x τ s > 0 and x σ s > 0 for s s 1 . From the definition of z, we obtain
x s z s p s x σ s z s p s z σ s 1 p s z s ,
with which with (1), we obtain (5). The proof is achieved. □

2.1. The Properties of the Solution in ρ 1

We will proceed to the first lemma, which analyses and provides details regarding the behavior of the positive solutions ρ 1 .
Lemma 4.
Let β * > 0 and x is a positive solution of (1) belonging to the class ρ 1 . Then, eventually:
(A 1 )
lim s z i s / s 3 i converges to 0 for i = 0 , 1 , 2 , 3 ;
(A 2 )
z s / s is decreasing;
(A 3 )
z s / s 2 is decreasing;
(A 4 )
z s / s 3 is decreasing;
Proof. 
Suppose x to the class ρ 1 , then for s 1 s 0 there is x s > 0 , x τ s > 0 and x σ s > 0 for s s 1 .
(A 1 ): Since we have z s , which is a non-increasing positive function, then z s 0 as s , if l > 0 then z s l > 0 , so
z s s s 1 3 3 ! ,
for s s 2 s 1 . From (4) and (5), we see that
z 4 s 3 ! β τ 3 s s z τ s .
From (7) into (8), we get
z 4 s β τ s s 1 3 τ 3 s s .
It is obvious that there exists s 3 > s 2 such that τ s s 1 3 1 2 τ 3 s for s s 3 , so we find from (9)
z 4 s β 2 s ,
for s s 3 . By integrating the above inequality from s 3 to s, we obtain
z s 3 z s + β 2 ln s s 3 + β 2 ln s s 3 as s ,
We find that there is a contradiction, therefore = 0 . We see when z ρ 1 that z s ,   z s as s , and also z s > 0 for i = 2 , is increasing such that z s / s 0 as s . Then, according to L’Hôpital’s rule, we find that (A 1 ) is satisfied.
(A 2 ): Since z s is non-increasing in ρ 1 , we see that
z s = z s 1 + s 1 s z ζ d ζ z s 1 + z s s s 1 t z s ,
where by (A 1 ) there is s 4 > s 3 such that z s 1 s 1 z s for s s 4 . So
z s s = z s s z s s 2 < 0 ,
for s s 4 , then z s / s is decreasing, which proves (A 2 ).
(A 3 ): From (A 1 ) and (A 2 ), z s / s decreases and tends to zero. Then, we find
z s = z s 4 + s 4 s z ζ d ζ z s 4 + z s s s 2 2 s 4 2 2 = z s 4 + z s s 2 z s s 4 2 t > z s s 2 , s s 5 ,
for s 5 > s 4 . Hence
z s s 2 = z s s 2 z s s 3 < 0 ,
for s s 5 . We arrive at (A 3 ).
(A 4 ): Likewise, since z s / s 2 is decreasing and tends to zero, we obtain
z s = z s 5 + s 5 s z ζ d ζ z s 5 + z s s 2 s 3 3 s 5 3 3 = z s s 3 , s s 6 ,
for s 6 > s 5 , so
z s s 3 = z s s 3 z s s 4 < 0 ,
for s s 6 . That proves (A 4 ).
As a result, the proof of the lemma is complete. □
In this lemma, we will establish some additional properties of the behavior of positive solutions in ρ 1
Lemma 5.
Assume x is a solution a positive of (1) belonging to the class ρ 1 and let β * > 0 . Then, for s large enough and every β 0 , β * :
(A 5 )
z s / s 1 β is decreasing;
(A 6 )
β < 1 ;
(A 7 )
lim s z i s / s 3 i β = 0 , i = 0 , 1 , 2 ;
(A 8 )
z s / s 2 β is decreasing;
(A 9 )
z s / s 3 β is decreasing.
Proof. 
Assume x is a positive solution of (1) to the class ρ 1 , then for s 1 s 0 there is x s > 0 , x τ s > 0 and x σ s > 0 for s s 1 , and (A 1 )–(A 2 ) of Lemma 4 holds for s s 1 .
(A 5 ): Define the positive function
ω s = z s s z s ,
by differentiating ω s and employing (8), in addition to having z s / s 3 decrease in (A 4 ), we get
ω s = s z 4 s 3 ! β z τ s τ 3 s > 3 ! β z s s 3 .
From (A 4 ) and (A 3 ) respectively, in the above inequality, we see
ω s > 2 β z s s 2 > β z s s .
By integrating the above inequality from s 1 to s, and using (A 1 ) and (A 2 ), we obtain
ω s > ω s 1 + β s 1 s z ζ ζ d ζ ω s 1 + β z s s s 1 s d ζ ,
for s s 2 , which is
ω s β z s ,
it follows from (11) that 1 β z s s z s for s s 2 , and hence
z s s 1 β < 0 , s s 2 .
We observe that z s / s 1 β is decreasing, thus (A 5 ) holds.
(A 6 ): Since z s is increasing, and from (A 5 ) we find that β < 1 , (A 6 ) thus holds.
(A 7 ): For i = 2 , to display that
lim s z s s 1 β = 0 ,
it suffices to prove that there is ε > 1 such that, for large s,
z s s 1 ε β < 0 ,
which if z s s 1 β c > 0 , then
z s s 1 β β ε 1 c s β ε 1 , as s .
We notice that there is a contradiction. We find that by using (14) for any η 2 β / 2 , 1 , then there is s 3 > s 2 large enough to obtain
z s = z s 2 + s 2 s z ζ ζ 1 β ζ 1 β d ζ z s 2 + z s s 1 β s 2 s ζ 1 β d ζ = z s 2 + z s s 1 β s 2 β s 2 2 β 2 β > k z s s 2 β ,
for s s 3 , by employing this in (13), we have
ω s > 2 β z s s 2 > 2 β k z s 2 β s .
Integrating the above inequality from s 3 to s and from (A 1 ) this forms
ω s > ω s 3 + 2 β k z s 2 β s s s 3 > 2 β μ 2 β z s ,
for s s 4 , from (11) in the above inequality; for s 4 > s 3 we arrive at
1 2 β k 2 β z s > s z s ,
so, from this, we can see that (16) is satisfied with
ε = 2 k 2 β > 1 .
For i = 0 , 1 the other limits in (A 7 ) are obtained from (15) and from using L’Hôpital’s rule.
(A 8 ): From (16) into (17), we see
z s z s 2 + z s s 2 β z s s 2 2 β s 1 β 2 β > z s s 2 β ,
for s s 5 , and for s 5 > s 4 we arrive at
z s s 2 β = s z s 2 β z s s 3 β < 0 ,
for s s 5 , when it is obvious that z s / s 2 β is decreasing, then (A 8 ) holds.
(A 9 ): We notice from (A 7 ) and (A 8 ) that z s / s 2 β is decreasing and tends to zero, then
z s z s 5 + s 5 s z ζ ζ 2 β ζ 2 β d ζ z s 5 + z s s 2 β s 3 β s 5 3 β 3 β = z s 5 + z s s 3 β z s s 5 3 β s 2 β 3 β > z s s 3 β ,
for s s 6 , and for s 6 > s 5 , we obtain
z s s 3 β = s z s 3 β z s s 4 β < 0 ,
for s s 6 . It follows that (A 9 ) holds, thus Lemma 5 is proved. □
The next lemma is a result of (A 9 ).
Lemma 6.
Assume that β * > 0 and δ * = then ρ 1 = .
Proof. 
Suppose the opposite is true and assume x is a positive solution of (1) to the class ρ 1 ; then, for s 1 s 0 there is x s > 0 , x τ s > 0 and x σ s > 0 for s s 1 . From (4) in (12), and taking into consideration (A 9 ), we get
ω s = s z 4 s 3 ! β z τ s τ 3 β s τ β s > 3 ! β z s s 3 s τ s β ,
which is
ω s > 3 ! β δ β z s s 3 .
Furthermore, from (A 9 ) and (A 8 ), respectively, we obtain
ω s > 3 ! β δ β z s 3 β s 2 > 3 ! β δ β z s 3 β 2 β s .
By integrating the above inequality from s 1 to s and from (A 1 ) and (A 2 ) in Lemma 5, with the definition of ω s in (11) for s 2 s 1 , we arrive at
1 3 ! β δ β 3 β 2 β z s > s z s ,
for s s 2 , where δ can be arbitrarily large, we chose it in such a way that
δ β > 3 β 2 β 3 ! β .
This indicates that z s > s z s , where z s and z s are positive, this causes a contradiction and completes the proof of Lemma 6. □
From (A 10 ) and Lemma 6 we can assume that δ * < , so ρ 1 .
The following lemma can be considered as an iterative version of Lemma 5.
Let us define a sequence β n (that is needed in the next lemma) as follows:
J 0 = J * , J m = 3 ! J 0 δ * J m 1 3 J m 1 2 J m 1 1 J m 1 ,
for m N , where J denotes β or γ .
By induction, it is simple to demonstrate that, if J i < 1 for i = 1 , 2 , . . . , m , then J m + 1 holds, such that
J m + 1 J m = J m > 1 ,
where the definitions of J m is as follows:
J 0 = J 1 J 0 = 3 ! δ * J 0 3 J m 2 J m 1 J m > 1 ,
J m = J m + 1 J m = δ * J m 3 J m 1 2 J m 1 1 J m 1 δ * J m 1 3 J m 2 J m 1 J m > 1 , m N .
We need to specify the sequence ε J n as follows:
ε J 0 = J J * < 1 ,
ε J m = ε J 0 δ ε J m 1 J m 1 3 J m 1 2 J m 1 1 J m 1 δ * J m 1 3 ε J m 1 J m 1 2 ε J m 1 J m 1 1 ε J m 1 J m 1 , m N .
The value of ε J n is arbitrary and determined by the value of β , where β is defined in (4). It is simple to show that
lim J J * ε J 0 = 1 , lim δ δ * J J * ε J m + 1 = 1 .
Lemma 7.
Assume x is a positive solution of (1) belonging to the class ρ 1 , and let β * > 0 . Then, for s large enough and ε β m ( 0 , 1 )
(A 10 )
m z s / s 1 β ˜ m is decreasing;
(A 11 )
m β ˜ m < 1 ;
(A 12 )
m lim s z i s / s 3 i β ˜ m = 0 , i = 0 , 1 , 2 ;
(A 13 )
m z s s / s 2 β ˜ m is decreasing;
(A 14 )
m z s / s 3 β ˜ m is decreasing,
where β ˜ m = ε β m β m
Proof. 
Assume x is a positive solution of (1) belonging to the class ρ 1 ; then, for s 1 s 0 there is x s > 0 , x τ s > 0 and x σ s > 0 for s s 1 . This Lemma will be proved by induction on m. For m = 0 , it holds from Lemma 5 with β = β ˜ 0 . After that, assume that (A 10 ) m (A 13 ) m hold for m 1 and s s m s 1 . We will display that they all hold m + 1 .
(A 10 ) m + 1 ; by using (4) and (A 8 ) in (12) we see
ω s 3 ! β ˜ 0 z τ s τ 3 s = 3 ! β ˜ 0 z τ s τ 3 β ˜ m s τ β ˜ m s 3 ! β ˜ 0 z s s 3 s τ s β ˜ m 3 ! β ˜ 0 δ β ˜ m z s s 3 ,
from (A 14 ) m and (A 13 ) m we obtain
ω s > 3 ! β ˜ 0 δ β ˜ m z s 3 β ˜ m s 2 > 3 ! β ˜ 0 δ β ˜ m z s 3 β ˜ m 2 β ˜ m s .
By integrating the above inequality from s m to s and from (A 10 ) m and (A 12 ) m , we find that there exists s m > s m such that
ω s ω s m + 3 ! β ˜ 0 δ β ˜ m 3 β ˜ m 2 β ˜ m s m s z ζ ζ 1 β ˜ m ζ β ˜ m d ζ ω s m + 3 ! β ˜ 0 δ β ˜ m z s 1 β ˜ m 3 β ˜ m 2 β ˜ m s 1 β ˜ m s 1 β ˜ m s m 1 β ˜ m > 3 ! β ˜ 0 δ β ˜ m 1 β ˜ m 3 β ˜ m 2 β ˜ m z s ,
which is
ω s > β ˜ m + 1 z s ,
from the definition of ω s , it follows that
1 β ˜ m + 1 z s > s z s ,
and we obtain that
z s s 1 β ˜ m + 1 < 0 .
That is the prove of (A 10 ) m + 1 .
(A 11 ) m + 1 : We have z s increasing and, from (A 10 ) m + 1 , we arrive at proving (A 11 ) m + 1 .
(A 12 ) m + 1 : To prove this case, it suffices to prove that there is ε > 1 , as done in the case m = 0 , such that for i = 2
z i s s 3 i ϵ β ˜ m + 1 < 0 .
From (23), we find that there is s m > s m sufficiently large such that
z s = z s m + s m s z ζ ζ 1 β ˜ m + 1 ζ 1 β ˜ m + 1 d ζ z s m + z s s 1 β ˜ m + 1 s m s ζ 1 β ˜ m + 1 d ζ = z s m + z s s 1 β ˜ m + 1 s 2 β ˜ m + 1 s m 2 β ˜ m + 1 2 β ˜ m + 1 ,
which is
z s > η 2 β ˜ m + 1 z s s , s s m ,
for any η ( 0 , 1 ) . If we merge the above inequality with (21) we obtain
ω s > 3 ! β ˜ 0 δ β ˜ m η 3 β ˜ m 2 β ˜ m + 1 z s s ,
by integration from s m to s and from (A 12 ) m we receive
ω s > ω s m + 3 ! β ˜ 0 δ β ˜ m η 3 β ˜ m 2 β ˜ m + 1 s m s z ζ ζ 1 β ˜ m ζ β ˜ m d ζ > 3 ! β ˜ 0 δ β ˜ m η 1 β ˜ m 3 β ˜ m 2 β ˜ m + 1 z s = η 2 β ˜ m 2 β ˜ m + 1 β ˜ m + 1 z s = ε β ˜ m + 1 z s ,
for s s m > s m . Since β ˜ m < β ˜ m + 1 , we can choose η such that ε > 1 , from (A 2 ) and the definition of ω we see that (16) is satisfied; the other limits are the same as those for m = 0 .
(A 13 ) m + 1 : By using (A 12 ) m + 1 in (25), we obtain
2 β ˜ m + 1 z s > z s s .
Then, (A 13 ) m + 1 is satisfied.
(A 14 ) m + 1 : From (A 7 ) and (A 8 ), we obtain
z s > z s m + s m s z ζ ζ 2 β ˜ m + 1 ζ 2 β ˜ m + 1 d ζ > z s 3 β ˜ m + 1 s 2 β ˜ m + 1 s 3 β ˜ m + 1 s m 3 β ˜ m + 1 > z s s 3 β ˜ m + 1 ,
which indicates (A 14 ) m + 1 holds and submits the lemma’s proof. □
The following lemma can be easily deduced from the aforementioned arguments.
Lemma 8.
Suppose that δ * < and
lim inf s τ 3 s s q s > ϱ 0 ,
where
ϱ 0 = max c ( 1 c ) ( 2 c ) ( 3 c ) δ * c : 0 < c < 1 1 p 0 .
Then ρ 1 = ⌀.
Proof. 
Suppose the opposite is true, that x ρ 1 . We claim that
β m 1 < 1 , m N .
From the case (A 11 ) m we have β ˜ m < 1 . Since ε β m 0 , 1 can be picked arbitrarily, set ε β m > 1 / β m , where β m is given by (19). Then
1 > β ˜ m = ε β m β m β m 1 > β m 1 ,
which supports the claim. From (28) we conclude that the sequence β m m = 0 is increasing and bounded from above, which is defined by (18), which means that
lim s β m = c ,
where c 0 , 1 is a root of the following equation:
c ( 1 c ) ( 2 c ) ( 3 c ) δ * c 1 p 0 = 3 ! β * .
However, from (26) we find that (29) has no positive solutions. As a result, ρ 1 = , which is the end of the proof and the Lemma. □

2.2. The Properties of the Solution in ρ 2

In this section, we will show results similar to the previous results of the section of solutions in the class ρ 2 .
Lemma 9.
Assume x is a positive solution of (1) relating to the class ρ 2 and let γ * > 0 . Then, for s large enough:
(A 15 )
lim s z i s / s 1 i converges to 0 for i = 0 , 1 ;
(A 16 )
z s / s is decreasing.
Proof. 
Suppose x is a positive solution of (1) belonging to the class ρ 2 , then, for s 1 s 0 , there is x s > 0 , x τ s > 0 and x σ s > 0 for s s 1 .
(A 15 ): As a result of the fact that z s is decreasing positive function, z s 0 as s , if l > 0 then z s l > 0 , so
z s s s 1 > s / 3 ,
for s s 2 s 1 . By (4) and (5) we obtain
z 4 s 3 ! γ τ s s 3 z τ s .
From (30) in the above inequality we obtain
z 4 s 2 γ / s 3 ,
we have obtained this by integrating twice from s to
z s γ s .
Integrating from s 2 to s we arrive at
z s 2 z s + γ ln s 2 s as s ,
there is, as we found, a contradiction. Therefore = 0 . We can show that (A 15 ) is satisfied by applying L’Hôpital’s rule.
(A 16 ): we have z s non-increasing and by (A 15 ) we see that
z s = z s 1 + s 1 s z ζ d ζ z s 1 + z s s s 1 t z s ,
for s s 3 > s 1 , where s 3 is sufficiently large such that z s 1 s 1 z s > 0 for s s 3 . Then
z s s = s z s z s s 2 ,
for s s 3 , the completion of the proof follows. □
Lemma 10.
Assume x is a positive solution of (1) belonging to the class ρ 2 and let γ * > 0 . Then, for s large enough and any γ 0 , γ * :
(A 17 )
z s / s 1 γ is decreasing;
(A 18 )
γ < 1 ;
(A 19 )
lim s z s / s 1 γ = 0 ;
(A 20 )
z s / s γ is non-decreasing.
Proof. 
Suppose x is a positive solution of (1) belonging to the class ρ 2 , then for s 1 s 0 there is x s > 0 , x τ s > 0 and x σ s > 0 for s s 1 .
(A 17 ): From (A 16 ) we have z s / s , which is decreasing; in (31) we obtain
z 4 s 3 ! γ τ s s 3 z τ s 3 ! γ z s s 4 ,
by integrating from s to twice such that z s is increasing, it goes as follows:
z s γ z s s 2 .
Let us define a positive function
μ s = z s s z s ,
by differentiating and from (32) we find
μ s = s z s γ z s s .
Integrating, again, s 1 to s, and we have z s / s decreasing and tending to zero, and we see
μ s μ s 1 + γ s 1 s z ζ ζ d ζ μ s 1 + γ z s s s s 1 > z s γ ,
for s s 2 , where s 2 > s 1 is sufficiently large such that μ s 1 z s s 1 / s > 0 for s s 2 . From the definition of μ s we obtain
1 γ z s s z s ,
and
z s s 1 γ = s z s 1 γ z s s 2 γ < 0 ,
s s 2 , then (A 17 ) holds.
(A 18 ): This simply implies from (A 17 ), and from the case that x is increasing.
(A 19 ): The proof is identical to the proof for class ρ 1 , and it is sufficient to show that
z s s 1 ε γ < 0 .
For ε > 1 , we can derive from (35) into (34), finding that there exists s 3 s 2 , such that
μ s μ s 2 + γ s 2 s z ζ ζ 1 γ ζ γ d ζ μ s 2 + γ z s 1 γ s 1 γ s 1 γ s 2 1 γ > η γ 1 γ z s , s s 3 ,
for any η 1 γ , 1 , from that we deduce that
1 η γ 1 γ z s s z s .
It is now obvious that (36) holds with ε = η / 1 γ > 1 .
(A 20 ): By integrating (32) from s to we see
z s γ s z s s 2 γ z s s ,
and so
z s s γ 0 ,
which is proof that z s / s γ is increasing. The lemma’s proof is now accomplished. □
Lemma 11.
Assume that γ * > 0 and δ * = then ρ 2 = .
Proof. 
Suppose x is a positive solution of (1) belonging to the class ρ 2 , then for s 1 s 0 there is x s > 0 , x τ s > 0 and x σ s > 0 for s s 1 . From (4) into (31), and taking into account (A 17 ), we arrive at
z 4 s 3 ! γ s 3 z τ s τ 1 γ s τ γ s 3 ! γ z s s 4 s τ s γ 3 ! γ δ γ z s s 4 .
By twice integrating from s to , with the assumption that x is increasing, we obtain
z s γ δ γ z s s 2 .
From the above inequality into (33)
μ s = s z s γ δ γ z s s .
Using integration as in (34), and replacing γ by γ δ γ , we obtain
1 γ δ γ z s s z s ,
where we can choose δ in such a way that it can be arbitrarily large, so that δ γ > 1 / γ , this shows z s = s z s —a contradiction. This illustrates the lemma. □
Now, from Lemma 10, we obtain an iterative.
Lemma 12.
Assume x is a positive solution of (1) belonging to the class ρ 2 and let γ * > 0 . Then, for s large enough and any ε γ m 0 , 1 :
(A 21 )
m z s / s 1 γ ˜ m is decreasing;
(A 22 )
m γ ˜ m < 1 ;
(A 23 )
m lim s z s / s 1 γ ˜ m = 0 ;
(A 24 )
m z s / s γ ˜ m is non-decreasing;
where γ ˜ m = ε γ m γ m .
Proof. 
Assume x is a positive solution of (1) belonging to the class ρ 2 , then for s 1 s 0 there is x s > 0 , x τ s > 0 and x σ s > 0 for s s 1 by induction on m. For m = 0 , it holds from Lemma 5 that γ = γ ˜ 0 . After that, assume that (A 21 ) m (A 24 ) m hold for m 1 and s s m s 1 . We will show that (A 21 ) m + 1 holds.
(A 21 ) m + 1 : From (31) and (A 21 ) m we find
z 4 s 3 ! γ ˜ 0 s 3 z τ s τ 1 γ ˜ m s τ γ ˜ m s 3 ! γ z s s 4 s τ s γ ˜ m 3 ! γ ˜ 0 δ γ ˜ m z s s 4 .
Integrate the above inequality from s to , and taking into account that z s / s γ ˜ m is increasing, yield
z s 3 ! γ ˜ 0 δ γ ˜ m s 2 s z ζ ζ 4 γ ˜ m ζ γ ˜ m d ζ 3 ! γ ˜ 0 δ γ ˜ m z ζ s γ ˜ m s 2 s d ζ ζ 4 γ ˜ m 3 ! γ δ γ ˜ m z s 3 γ ˜ m s 3 .
Repeating this process, we obtain
z s 3 ! γ δ γ ˜ m z s 3 γ ˜ m 2 γ ˜ m s 2 ,
By incorporating this into (33), we have
μ s = s z s 3 ! γ δ γ ˜ m z s 3 γ ˜ m 2 γ ˜ m s ,
By integrating from s m to s and from (A 21 ) m and (A 23 ) m , we find
μ s m μ s + 3 ! γ 0 δ γ ˜ m 3 γ ˜ m 2 γ ˜ m s m s z ζ ζ 1 γ ˜ m ζ γ ˜ m d ζ 3 ! γ 0 δ γ ˜ m z s 3 γ ˜ m 2 γ ˜ m 1 γ ˜ m s 1 γ ˜ m s 1 γ ˜ m s m 1 γ ˜ m γ ˜ m + 1 z s , s s m ,
That is the proof of (A 21 ) m + 1 . The other parts’ proofs of the lemma are the same as those in the case where m = 0 . □
Lemma 13.
Suppose that δ * < and
lim inf s τ s s 3 q s > ϱ 0 ,
where
ϱ 0 = max c ( 1 c ) ( 2 c ) ( 3 c ) α c 3 1 p 0 ,
for 0 < c < 1 , then ρ 2 = ⌀.
Proof. 
The proof is similar to the proof of Lemma 8 and defines γ m as in (18). □
Now, in the next theorem, we offer the fundamental result in this work by combining the results from the previous two sections.

2.3. Oscillation Results

Theorem 1.
Suppose that
lim inf s τ 3 s s q s > 0 f o r δ * = , ϱ 0 f o r δ * < ,
where
ϱ 0 = max c ( 1 c ) ( 2 c ) ( 3 c ) δ * c 1 p 0 ,
where 0 < c < 1 , then (1) is oscillatory.
Proof. 
From (41) observe that β * > 0 , and since
lim inf s τ 3 s s q s lim inf s τ s s 3 q s ,
we find that γ * > 0 . Now, if δ * = , then Lemmas 6 and 11 imply that ρ 2 = ρ 1 = . For δ * < , from Lemmas 8 and 13, the same conclusion is derived. This illustrates the theorem. □
Corollary 1.
Let τ s = α s with 0 < α 1 . If
lim inf s s 4 q s > max c ( 1 c ) ( 2 c ) ( 3 c ) α c 3 : 0 < c < 1 1 p 0 ,
then (1) is oscillatory.

2.4. Application and Discussion

In the next section, we provide an example to highlight our study results.
Example 1.
Now, consider the fourth-order Euler delay differential equation
x s + p 0 x σ s 4 s + q 0 δ * 3 s 4 x 1 δ * s = 0 ,
for s > 1 , where p 0 > 0 , q 0 > 0 and δ * 1 , By applying condition (41), we obtain
lim inf s s δ * 3 s q 0 δ * 3 s 4 > max c ( 1 c ) ( 2 c ) ( 3 c ) δ * c 1 p 0 ,
which is,
q 0 > max c ( 1 c ) ( 2 c ) ( 3 c ) δ * c 1 p 0 ,
where 0 < c < 1 . Thus, by applying Theorem 1, we can guarantee that all solutions of Equation (42) are oscillatory if condition (43) is satisfied.
Remark 1.
If we consider the special case p 0 = 0.5 and δ * = 2 , the condition (43) reduces to
q 0 > 1.785 .
By checking the result of the oscillation constants for Equation (42) in references [20] and [21], respectively, with p 0 = 0.5 and δ * = 2 , we see
q 0 > 96 e ln 2 50.951 ,
and
q 0 > 82 9 9.111 .
Example 2.
Consider the NDE as the following:
x s + 3 4 x 5 s 2 4 s + 4 3 q 0 s 4 x s 4 = 0 ,
where s > 1 , p 0 = 3 4 , q s = 4 3 q 0 s 4 and δ * = 4 , σ s = 5 s 2 . To check the oscillation of Equation (47), we will apply condition (41) of Theorem 1 in the previous section and see that
lim inf s s 4 3 s 64 q 0 s 4 > max c ( 1 c ) ( 2 c ) ( 3 c ) 4 c 1 p 0 ,
which is
q 0 > 2.561 .
Applying results in both [22] and [23] to Equation (47), we get, respectively,
q 0 > 832 ,
and
q 0 > 675 14 e ln 10 8.627 .
We notice that the condition (48) improves the condition (49) and (50). It also improves results (45) and (46).
Remark 2.
It can be easily observed that condition (44) and (48) improve conditions (45), (46), (49), and (50). In addition to this improvement, there is something that distinguishes our results from other results in [20] and [21] is that their results require constraints τ s < s , τ s 0 but Theorem 1 does not need them. This leads to the conclusion that Theorem 1 improves many previous results in the literature, even without the usual restrictive suppositions on the diverging argument.
Remark 3.
We make a simulation experiment for Example 1, by considering the ODE of (42)
x s + p 0 x s 4 s + q 0 δ * 3 s 4 x s = 0 ,
where p 0 = 0.5   , q 0 = q 0 δ * 3 s 4 and δ * = 1 . By using Theorem 1, Equation (51) is oscillatory if
lim inf s s 3 s q 0 s 4 > max c ( 1 c ) ( 2 c ) ( 3 c ) 1 p 0 ,
which is
q 0 > max c ( 1 c ) ( 2 c ) ( 3 c ) 1 p 0 , where 0 < c < 1 .

3. Conclusions

The oscillatory behavior of the solutions of the equation of the neutral type was studied, where the positive solutions of the equation were classified as ρ 1 and ρ 2 , and then we studied the monotonic properties of these positive solutions by providing a series of lemmas for each case to iteratively improve the monotonicity features of non-oscillatory solutions. We observe that the main difference between ρ 1 and ρ 2 is the change of the second derivative of the function z, but this simple change constitutes a change in many of the monotonic and asymptotic properties. For example, in case ρ 1 , we can get new monotonic properties of the first, second, and third derivatives of the function z, while in case ρ 2 we can only get an ordinal behavior of z; there is no change except for the function z only. By using the iterative method based on these characteristics, we have provided a single criterion to eliminate the positive solutions to our equation that ensures the oscillatory nature of the solutions. The example improved upon previous results and showed the importance of the new properties. The most important thing that distinguishes our results is that they can be applied without some of the restrictions that other results require.

Author Contributions

Conceptualization, H.S., O.M., S.S.A., A.M.A. and E.M.E.; Methodology, H.S., O.M., S.S.A., A.M.A. and E.M.E.; Investigation, H.S., O.M., S.S.A., A.M.A. and E.M.E.; Writing—original draft, O.M., S.S.A. and A.M.A.; Writing—review & editing, H.S. and E.M.E. All authors have read and agreed to the published version of the manuscript.

Funding

This project is funded by King Saud University, Riyadh, Saudi Arabia.

Acknowledgments

The authors present their appreciation to King Saud University for funding the publication of this research through the Researchers Supporting Program (RSPD2023R533), King Saud University, Riyadh, Saudi Arabia.

Conflicts of Interest

This work does not have any conflict of interest.

References

  1. Rihan, F.A. Delay Differential Equations and Applications to Biology; Springer Nature Singapore Pte Ltd.: Singapore, 2021. [Google Scholar] [CrossRef]
  2. Hale, J.K. Functional Differential Equations. In Oxford Applied Mathematical Sciences; Springer: New York, NY, USA, 1971. [Google Scholar]
  3. Gyori, I.; Ladas, G. Oscillation Theory of Delay Differential Equations with Applications; Clarendon Press: Oxford, UK, 1991. [Google Scholar]
  4. Erbe, L.H.; Kong, Q.; Zhong, B.G. Oscillation Theory for Functional Differential Equations; Marcel Dekker: New York, NY, USA, 1995. [Google Scholar]
  5. Ladde, G.S.; Lakshmikantham, V.; Zhang, B.G. Oscillation Theory of Differential Equations with Deviating Arguments; Marcel Dekker: New York, NY, USA, 1987. [Google Scholar]
  6. Grace, S.R.; Jadlovska, I.; Zafer, A. On oscillation of second order delay differential equations with a sublinear neutral term. Mediterr. J. Math. 2020, 17, 116. [Google Scholar] [CrossRef]
  7. Jadlovská, I. New criteria for sharp oscillation of second-order neutral delay differential equations. Mathematics 2021, 9, 2089. [Google Scholar] [CrossRef]
  8. Agarwal, R.P.; Grace, S.R.; O’Regan, D. Oscillation criteria for certain nth order differential equations with deviating arguments. J. Math. Anal. Appl. 2001, 262, 601–622. [Google Scholar] [CrossRef]
  9. Grace, S.R.; Agarwal, R.P.; Pavani, R.; Thandapani, E. On the oscillation of certain third order nonlinear functional differential equations. Appl. Math. Comput. 2008, 202, 102–112. [Google Scholar] [CrossRef]
  10. Agarwal, R.P.; Zhang, C.; Li, T. Some remarks on oscillation of second order neutral differential equations. Appl. Math. Comput. 2016, 274, 178–181. [Google Scholar] [CrossRef]
  11. Li, T.; Baculìková, B.; Džurina, J.; Zhang, C. Oscillation of fourth-order neutral differential equations with p-Laplacian like operators. Bound. Value Probl. 2014, 2014, 56. [Google Scholar] [CrossRef]
  12. Grace, S.R.; Džurina, J.; Jadlovská, I.; Li, T. An improved approach for studying oscillation of second-order neutral delay differential equations. J. Inequalities Appl. 2018, 2018, 193. [Google Scholar] [CrossRef] [PubMed]
  13. Agarwal, R.; Grace, S.; Manojlovic, J. Oscillation criteria for certain fourth order nonlinear functional differential equations. Math. Comput. Model. 2006, 44, 163–187. [Google Scholar] [CrossRef]
  14. Cesarano, C.; Bazighifan, O. Oscillation of fourth-order functional differential equations with distributed delay. Axioms 2019, 8, 61. [Google Scholar] [CrossRef]
  15. Hale, J.K. Theory of Functional Differential Equations; Springer: New York, NY, USA, 1977. [Google Scholar] [CrossRef]
  16. Bartušek, M.; Cecchi, M.; Došlá, Z.; Marini, M. Fourth-order differential equation with deviating argument. Abstr. Appl. Anal. 2012, 2012, 185242. [Google Scholar] [CrossRef]
  17. Chatzarakis, G.E.; Elabbasy, E.M.; Bazighifan, O. An oscillation criterion in 4th-order neutral differential equations with a continuously distributed delay. Adv. Differ. Equ. 2019, 2019, 336. [Google Scholar] [CrossRef]
  18. Bazighifan, O.; Ruggieri, M.; Scapellato, A. An improved criterion for the oscillation of fourth-order differential equations. Mathematics 2020, 8, 610. [Google Scholar] [CrossRef]
  19. Kiguradze, I.; Chanturia, T.A. Asymptotic Properties of Solutions of Nonautonomous Ordinary Differential Equations; Kluwer Academic, Publishers: Dordrecht, The Netherlands, 1993. [Google Scholar] [CrossRef]
  20. Zafer, A. Oscillation criteria for even order neutral differential equations. Appl. Math. Lett. 1998, 11, 21–25. [Google Scholar] [CrossRef]
  21. Agarwal, R.P.; Bohner, M.; Li, T.; Zhang, C. A new approach in the study of oscillatory behavior of even-order neutral delay differential equations. Appl. Math. Comput. 2013, 225, 787–794. [Google Scholar] [CrossRef]
  22. Cesarano, C.; Moaaz, O.; Qaraad, B.; Muhib, A. Oscillatory and asymptotic properties of higher-order quasilinear neutral differential equations. AIMS Math. 2021, 6, 11124–11138. [Google Scholar] [CrossRef]
  23. Moaaz, O.; Alhgilan, A.E. A Study of the Monotonic Properties of Solutions of Neutral Differential Equations and Their Applications. Axioms 2023, 12, 346. [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

Salah, H.; Moaaz, O.; Askar, S.S.; Alshamrani, A.M.; Elabbasy, E.M. Optimizing the Monotonic Properties of Fourth-Order Neutral Differential Equations and Their Applications. Symmetry 2023, 15, 1744. https://doi.org/10.3390/sym15091744

AMA Style

Salah H, Moaaz O, Askar SS, Alshamrani AM, Elabbasy EM. Optimizing the Monotonic Properties of Fourth-Order Neutral Differential Equations and Their Applications. Symmetry. 2023; 15(9):1744. https://doi.org/10.3390/sym15091744

Chicago/Turabian Style

Salah, Hend, Osama Moaaz, Sameh S. Askar, Ahmad M. Alshamrani, and Elmetwally M. Elabbasy. 2023. "Optimizing the Monotonic Properties of Fourth-Order Neutral Differential Equations and Their Applications" Symmetry 15, no. 9: 1744. https://doi.org/10.3390/sym15091744

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