This paper investigates the oscillatory behavior of nonlinear third-order dynamic equations on time scales. Our main approach is to transform the equation from its semi-canonical form into a more tractable canonical form. This transformation simplifies the analysis of oscillation behavior and allows us to derive new oscillation criteria. These criteria guarantee that all solutions to the equation oscillate. Our results extend and improve upon existing findings in the literature, particularly for the special cases where T=R and T=Z. Additionally, we provide illustrative examples to demonstrate the practical application of the developed criteria.
Citation: Ahmed M. Hassan, Clemente Cesarano, Sameh S. Askar, Ahmad M. Alshamrani. Oscillatory behavior of solutions of third order semi-canonical dynamic equations on time scale[J]. AIMS Mathematics, 2024, 9(9): 24213-24228. doi: 10.3934/math.20241178
[1] | A. A. El-Gaber, M. M. A. El-Sheikh, M. Zakarya, Amirah Ayidh I Al-Thaqfan, H. M. Rezk . On the oscillation of solutions of third-order differential equations with non-positive neutral coefficients. AIMS Mathematics, 2024, 9(11): 32257-32271. doi: 10.3934/math.20241548 |
[2] | Maged Alkilayh . On The oscillatory behavior of solutions to a class of second-order nonlinear differential equations. AIMS Mathematics, 2024, 9(12): 36191-36201. doi: 10.3934/math.20241718 |
[3] | Samy E. Affan, Elmetwally M. Elabbasy, Bassant M. El-Matary, Taher S. Hassan, Ahmed M. Hassan . Second-order advanced dynamic equations on time scales: Oscillation analysis via monotonicity properties. AIMS Mathematics, 2025, 10(2): 4473-4491. doi: 10.3934/math.2025206 |
[4] | S. Sangeetha, S. K. Thamilvanan, S. S. Santra, S. Noeiaghdam, M. Abdollahzadeh . Property $ \bar{A} $ of third-order noncanonical functional differential equations with positive and negative terms. AIMS Mathematics, 2023, 8(6): 14167-14179. doi: 10.3934/math.2023724 |
[5] | Mohammed Ahmed Alomair, Ali Muhib . On the oscillation of fourth-order canonical differential equation with several delays. AIMS Mathematics, 2024, 9(8): 19997-20013. doi: 10.3934/math.2024975 |
[6] | Abdulaziz khalid Alsharidi, Ali Muhib . Oscillation criteria for mixed neutral differential equations. AIMS Mathematics, 2024, 9(6): 14473-14486. doi: 10.3934/math.2024703 |
[7] | Yuanbin Ma, Zhi Li . Perturbed uncertain differential equations and perturbed reflected canonical process. AIMS Mathematics, 2021, 6(9): 9647-9659. doi: 10.3934/math.2021562 |
[8] | Osama Moaaz, Wedad Albalawi . Differential equations of the neutral delay type: More efficient conditions for oscillation. AIMS Mathematics, 2023, 8(6): 12729-12750. doi: 10.3934/math.2023641 |
[9] | Dmitry Sinelshchikov . On an integrability criterion for a family of cubic oscillators. AIMS Mathematics, 2021, 6(11): 12902-12910. doi: 10.3934/math.2021745 |
[10] | Abdulaziz khalid Alsharidi, Ali Muhib . Functional differential equations in the non-canonical case: New conditions for oscillation. AIMS Mathematics, 2025, 10(3): 7256-7268. doi: 10.3934/math.2025332 |
This paper investigates the oscillatory behavior of nonlinear third-order dynamic equations on time scales. Our main approach is to transform the equation from its semi-canonical form into a more tractable canonical form. This transformation simplifies the analysis of oscillation behavior and allows us to derive new oscillation criteria. These criteria guarantee that all solutions to the equation oscillate. Our results extend and improve upon existing findings in the literature, particularly for the special cases where T=R and T=Z. Additionally, we provide illustrative examples to demonstrate the practical application of the developed criteria.
A recently developed theory gaining significant interest is time scale T, introduced by Stefan Hilger to bridge the gap between continuous and discrete analysis [1]. In simpler terms, it aims to unify the study of differential equations (governing continuous change) and difference equations (modeling discrete jumps) [2].
The core concept of time scale T involves defining a time domain as any non-empty, closed set of real numbers. The familiar differential and difference equations emerge as special cases when the time scale is the set of all real numbers or integers, respectively.
To gain a comprehensive understanding, it is necessary to review some basic concepts of time scale theory. The forward and backward jump operators σ, ρ:T→T are defined by
σ(ℓ)=inf{s∈T∣s>ℓ} and ρ(ℓ)=sup{s∈T∣s<ℓ}, |
(supplemented by inf∅=supT and sup∅=infT). A point ℓ∈T is called right-scattered, right-dense, left-scattered or left-dense if σ(ℓ)>ℓ, σ(ℓ)=ℓ, ρ(ℓ)<ℓ, ρ(ℓ)=ℓ holds, respectively. The set Tκ is defined to be T if T does not have a left-scattered maximum; otherwise, it is T without this left-scattered maximum. The graininess function μ:T→[0,∞) is defined by μ(ℓ)=σ(ℓ)−ℓ. Hence, the graininess function is constant 0 if T=R, while it is constant ℓ for T=Z. However, a time scale T could have nonconstant graininess. A function h:T→R is said to be rd-continuous and is written h∈Crd(T,R), provided that h is continuous at right dense points and at left dense points in T, left hand limits exist, and are finite. We say that h:T→R is differentiable at ℓ∈T whenever
hΔ:=lims→ℓh(ℓ)−h(s)ℓ−s |
exists when σ(ℓ)=ℓ (here by s→ℓ it is understood that s approaches ℓ in the time scale), and when h is continuous at ℓ and σ(ℓ)>ℓ it is
hΔ:=lims→ℓh(σ(ℓ))−h(ℓ)μ(ℓ). |
The product and quotient rules [3, Theorem 1.20] for the derivative of the product hk and the quotient h/k of two differentiable functions h and k are as follows:
(hk)Δ(ℓ)=hΔ(ℓ)k(ℓ)+h(σ(ℓ))kΔ=h(ℓ)kΔ(ℓ)+hΔ(ℓ)k(σ(ℓ)), | (1.1) |
(hk)Δ(ℓ)=hΔ(ℓ)k(ℓ)−h(ℓ)kΔ(ℓ)k(ℓ)k(σ(ℓ)). | (1.2) |
The chain rule [3, Theorem 1.90] for the derivative of the composite function h∘k of a continuously differentiable function h:R→R and a (delta) differentiable function k:T→R results in
(h∘k)Δ={∫10h′(k+sμkΔ)ds}gΔ. | (1.3) |
A function h:T→R is called rd-continuous provided it is continuous at right-dense points in T and its left-sided limits exist (finite) at left-dense points in T. The set of rd-continuous functions f:T→R is denoted by
Crd=Crd(T)=Crd(T,R). |
The set of functions h:T→R that are differentiable and whose derivative is rd-continuous is denoted by
C1rd=C1rd(T)=C1rd(T,R). |
Finally, if h:T→R is a function, then we define the function hσ:T→R by hσ(ℓ)=h(σ(ℓ)) for all ℓ∈T.
Bohner and Peterson's book [3] provides a comprehensive overview and organization of this new calculus. Beyond these basic cases, numerous other time scales can be defined, leading to a wealth of applications. One such application is the study of population dynamic models, as explored in [4]. To delve deeper into the theory, readers can consult the referenced papers [5,6] and monographs [3,7].
Recent years have seen a surge in research on the oscillation and non-oscillation of solutions to dynamic equations on time scales. For further exploration, readers can refer to the references provided [8,9,13,14,15,16,17,18,19].
The present paper investigates the asymptotic behavior of solutions to the semi-canonical third dynamic equation
(a2(ℓ)(a1(ℓ)xΔ(ℓ))Δ)Δ+p(ℓ)xγ(δ(ℓ))=0, ℓ∈[ℓ0,∞)T, | (1.4) |
where γ is the ratio of positive odd integers.
In this paper, we consider the following conditions:
(ⅰ) a1(ℓ)∈C2rd([ℓ0,∞)T,(0,∞)), a2(ℓ)∈C1rd([ℓ0,∞)T,(0,∞)), p(ℓ)∈C([ℓ0,∞)T,(0,∞)) and Eq (1.4) is in semi-canonical form, i.e.,
∫∞ℓ0Δsa2(s)=∞ and ∫∞ℓ0Δsa1(s)<∞; | (1.5) |
(ⅱ) δ∈C1rd([ℓ0,∞)T), δΔ(ℓ)≥0, and limℓ→∞δ(ℓ)=∞.
Let us recall that a solution of Eq (1.4) is a nontrivial real-valued function x satisfying the equation for ℓ≥ℓx for some ℓx≥ℓx0 such that x∈C1([ℓx,∞)T,R), a1xΔ∈C1([ℓx,∞)T,R), and a2(a1(xΔ))Δ∈C1([ℓx,1)T,R). We exclude solutions that vanish identically in some neighborhood of infinity, assuming that such solutions exist for Eq (1.4). A solution x(ℓ) of Eq (1.4) is termed oscillatory if it exhibits arbitrarily large zeros on [ℓx,∞)T; otherwise, it is classified as non-oscillatory.
The study of oscillatory behavior in Eq (1.4) often hinges on its form. Equation (1.4) is in canonical form if
∫∞ℓ0Δsa1(s)=∫∞ℓ0Δsa2(s)=∞, |
and it is in non-canonical form if
∫∞ℓ0Δsa1(s)<∞ and ∫∞ℓ0Δsa2(s)<∞. |
If either
∫∞ℓ0Δsa1(s)<∞ and ∫∞ℓ0Δsa2(s)=∞, | (S1) |
or
∫∞ℓ0Δsa1(s)=∞ and ∫∞ℓ0Δsa2(s)<∞, | (S2) |
then we will say that (1.4) is in semicanonical form.
The groundwork for studying third-order dynamic equations on general time scales was laid by Erbe et al. [13], who focused on equations of the form
(a2(ℓ)((a1(ℓ)xΔ(ℓ))Δ))Δ+p(ℓ)f(x(ℓ))=0, ℓ∈[ℓ0,∞)T, | (1.6) |
where a1,a2,p∈Crd(ℓ0,∞)T, f∈C(R,R) R is continuous and satisfies uf(u)>0 for u≠0. Additionally, for each k>0, there exists M=Mk>0 such that f(u)/u≥M, |u|≥k. Using the Riccati transformation technique, they established sufficient conditions that guarantee every solution to this equation either oscillates or converges to zero.
Building on Erbe et al.'s work [13], Hassan [14] investigated a more general form of the third-order equation
(a2(ℓ)((a1(ℓ)xΔ(ℓ))Δ)α)Δ+f(ℓ,x(δ(ℓ)))=0, ℓ∈[ℓ0,∞)T, | (1.7) |
where α≥1 and δ(ℓ)≤ℓ, in the canonical form.
In the particular case of T=R and γ=1, Chatzarakis et al. [11] established new oscillation criteria for the differential equation
(a2(ℓ)(a1(ℓ)x′(ℓ))′)′+p(ℓ)x(δ(ℓ))=0, |
in the canonical form. Recently, techniques have been developed to study the oscillatory behavior of solutions to third-order equations. Moaaz et al. [21,22] extended the improved methods used in studying second-order equations [23,24]. The development of oscillation criteria for delay differential equations of odd orders can also be observed through the works [25,26].
Our literature review indicates a scarcity of research on the oscillatory behavior of solutions to Eq (1.4) when it takes the semi-canonical form (S1). This paper tackles Eq (1.4) in its less-studied semi-canonical form. We begin by transforming it into the more common canonical form. This transformation allows us to then establish new criteria for determining when solutions to Eq (1.4) oscillate.
To enhance readability, we'll use the following symbols:
A(ℓ):=∫∞ℓΔsa1(s), a(ℓ):=a1(ℓ)A(ℓ)Aσ(ℓ), r(ℓ):=a2(ℓ)Aσ(ℓ), |
P(ℓ):=p(ℓ)Aγ(ℓ), ϕ(ℓ):=∫ℓℓ1Δsr(s), ψ(ℓ):=∫ℓℓ1ϕ(s)Δsa(s), and z(ℓ):=x(ℓ)A(ℓ). |
Lemma 2.1. [27] Assume that x is an eventually positive solution of (1.4) satisfying (1.5). Then there exists ℓ1∈[ℓ0,∞)T such that x satisfies one of the following three cases:
(I) xΔ>0, (a1(ℓ)(xΔ(ℓ)))Δ>0, (a2(ℓ)((a1(ℓ)(xΔ(ℓ)))Δ)Δ<0;
(II) xΔ<0, (a1(ℓ)(xΔ(ℓ)))Δ>0, (a2(ℓ)((a1(ℓ)(xΔ(ℓ)))Δ)Δ<0;
(III) xΔ<0, (a1(ℓ)(xΔ(ℓ)))Δ<0, (a2(ℓ)((a1(ℓ)(xΔ(ℓ)))Δ)Δ<0.
Theorem 2.1. Assume that
∫∞ℓΔsr(s)=∞. | (2.1) |
Then the semi-canonical dynamic Eq (1.4) has a solution x(ℓ) if and only if the corresponding canonical equation
(r(ℓ)(a(ℓ)zΔ(ℓ))Δ)Δ+P(ℓ)zγ(δ(ℓ))=0, | (2.2) |
admits the solution z(ℓ)=x(ℓ)A(ℓ).
Proof. Referring back to σ(ℓ) as the forward jump operator and performing differentiation yields
a2(ℓ)Aσ(ℓ)(a1(ℓ)A(ℓ)Aσ(ℓ)(x(ℓ)A(ℓ))Δ)Δ=a2(ℓ)Aσ(ℓ){a1(ℓ)A(ℓ)Aσ(ℓ)[xΔ(ℓ)A(t)−x(t)AΔ(t)A(t)Aσ(ℓ)]}Δ=a2(ℓ)Aσ(ℓ){a1(ℓ)xΔ(ℓ)A(ℓ)−a1(ℓ)x(ℓ)AΔ(ℓ)}Δ=a2(ℓ)Aσ(ℓ){(a1(ℓ)xΔ(ℓ))ΔAσ(ℓ)+a1(ℓ)xΔ(ℓ)AΔ(ℓ)+xΔ(ℓ)}=a2(ℓ)(a1(ℓ)xΔ(ℓ))Δ. | (2.3) |
From (2.1), we have
∫∞ℓAσ(s)a2(s)Δs=∞, | (2.4) |
and
∫∞ℓΔsa1(s)A(s)Aσ(s)=∫∞ℓ(1A(s))ΔΔs=limℓ→∞(1A(ℓ)−1A(ℓ0))=∞. | (2.5) |
Combining (2.3) with (1.4), we obtain
(a2(ℓ)(a1(ℓ)xΔ(ℓ))Δ)Δ+p(ℓ)xγ(δ(ℓ))=0(a2(ℓ)Aσ(ℓ)(a1(ℓ)A(ℓ)Aσ(ℓ)(x(ℓ)A(ℓ))Δ)Δ)Δ+p(ℓ)Aγ(ℓ)xγ(δ(ℓ))Aγ(ℓ)=0(r(ℓ)(a(ℓ)zΔ(ℓ))Δ)Δ+P(ℓ)zγ(δ(ℓ))=0. | (2.6) |
It is clear that x(ℓ)A(ℓ) is a solution of (2.6). Moreover, considering (2.4) and (2.5), it is apparent that Eq (2.6) is in canonical form and from [28] this canonical from is unique.
Theorem 2.1 significantly streamlines the analysis of Eq (1.4) by reducing it to the scope of (2.2), thereby directing our focus towards only two classifications of solutions that ultimately exhibit positivity, i.e., either
z(ℓ)>0, a(ℓ)zΔ(ℓ)<0, r(ℓ)(a(ℓ)zΔ(ℓ)>0, (r(ℓ)(a(ℓ)zΔ(ℓ))Δ)Δ<0, |
and in this case, we denote z∈ℵ0 or
z(ℓ)>0, a(ℓ)zΔ(ℓ)>0, r(ℓ)(a(ℓ)zΔ(ℓ))Δ>0, (r(ℓ)(a(ℓ)zΔ(ℓ))Δ)Δ<0, |
and for this characteristic, we indicate that z∈ℵ2.
Theorem 2.2. Let γ≥1 and (2.1) hold. Suppose that
lim supℓ→∞{1ϕγ(δ(ℓ))∫δ(ℓ)ℓ1ϕ(σ(s))P(s)ψγ(s)Δs+1ϕγ−1(δ(ℓ))∫ℓδ(ℓ)P(s)ψγ(s)ϕγ(δ(s))Δs +ϕ(δ(ℓ))∫∞ℓP(s)ψγ(s)Δs}={∞,γ>1,1,γ=1, | (2.7) |
and
∫∞ℓ01a(u)∫∞u1r(v)∫∞vP(s)ΔsΔvΔu=∞. | (2.8) |
Then every non-oscillatory solution z(ℓ) of (1.4) satisfies limℓ→∞x(ℓ)A(ℓ)=0.
Proof. Let x(ℓ) be a non-oscillatory solution of Eq (1.4), where x(ℓ)>0, and x(δ(ℓ))>0 for ℓ≥ℓ1 for some ℓ1≥ℓ0. According to Theorem 2.1, the corresponding function z(ℓ)=x(ℓ)A(ℓ) is a positive solution of (2.2), implying that either z∈ℵ0 or z∈ℵ2 for ℓ≥ℓ1.
Let us examine the case where z∈ℵ2. In this case, we observe that
a(ℓ)zΔ(ℓ)≥∫ℓℓ1r−1(s)r(s)(a(s)zΔ(s))ΔΔs≥r(ℓ)(a(ℓ)zΔ(ℓ))Δ∫ℓℓ1Δsr(s)≥r(ℓ)(a(ℓ)zΔ(ℓ))Δϕ(ℓ). |
Hence,
(a(ℓ)zΔ(ℓ)ϕ(ℓ))Δ=ϕ(ℓ)(a(ℓ)zΔ(ℓ))Δ−(a(ℓ)zΔ(ℓ))ϕΔ(ℓ)ϕ(ℓ)ϕσ(ℓ)=ϕ(ℓ)r(ℓ)(a(ℓ)zΔ(ℓ))Δ−(a(ℓ)zΔ(ℓ))r(ℓ)ϕ(ℓ)ϕσ(ℓ)≤0. | (2.9) |
Consequently, it can be inferred from (2.9) that
z(ℓ)≥∫ℓℓ1zΔ(s)Δs=∫ℓℓ1a(s)zΔ(s)ϕ(s)ϕ(s)a(s)Δs≥a(ℓ)zΔ(ℓ)ϕ(ℓ)ψ(ℓ). | (2.10) |
Combining (2.10) with (2.2), we see that a(ℓ)zΔ(ℓ)ϕ(ℓ)ψ(ℓ) is a positive solution to the dynamic inequality
(r(ℓ)χΔ(ℓ))Δ+P(ℓ)ψγ(ℓ)ϕγ(δ(ℓ))χγ(δ(ℓ))≤0, | (2.11) |
where χ(ℓ):=a(ℓ)zΔ(ℓ). Integration (2.11) from ℓ to ∞ and considering the nonincreasing nature of χ(ℓ)/ϕ(ℓ), we obtain
χΔ(ℓ)≥1r(ℓ)∫∞ℓP(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs. |
Therefore,
χ(ℓ)≥∫ℓℓ11r(s)∫∞sP(u)ψγ(u)ϕγ(δ(u))χγ(δ(u))ΔuΔs=∫ℓℓ11r(s)∫ℓsP(u)ψγ(u)ϕγ(δ(u))χγ(δ(u))ΔuΔs+∫ℓℓ11r(s)∫∞ℓP(u)ψγ(u)ϕγ(δ(u))χγ(δ(u))ΔuΔs. | (2.12) |
Integrating by parts, we obtain
χ(ℓ)≥∫ℓℓ1ϕ(σ(s))P(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs+ϕ(ℓ)∫∞ℓP(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs. | (2.13) |
It follows that
χ(δ(ℓ))≥∫δ(ℓ)ℓ1ϕ(σ(s))P(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs+ϕ(δ(ℓ))∫∞δ(ℓ)P(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs=∫δ(ℓ)ℓ1ϕ(σ(s))P(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs+ϕ(δ(ℓ))∫ℓδ(ℓ)P(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs +ϕ(δ(ℓ))∫∞ℓP(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs. | (2.14) |
Utilizing the monotonicity characteristics of χ(ℓ) and χ(ℓ)/ϕ(ℓ), we have χ(δ(ℓ))≤χ(δ(s)) and χ(δ(s))ϕ(δ(s))≥χ(δ(ℓ))ϕ(δ(ℓ)) for s≥ℓ, hence (2.13) takes the form
χ(δ(ℓ))≥χγ(δ(ℓ))ϕγ(δ(ℓ))∫δ(ℓ)ℓ1ϕ(σ(s))P(s)ψγ(s)Δs+χγ(δ(ℓ))ϕγ−1(δ(ℓ))∫ℓδ(ℓ)P(s)ψγ(s)ϕγ(δ(s))Δs +ϕ(δ(ℓ))χγ(δ(ℓ))∫∞ℓP(s)ψγ(s)Δs, | (2.15) |
χ1−γ(δ(ℓ))≥1ϕγ(δ(ℓ))∫δ(ℓ)ℓ1ϕ(σ(s))P(s)ψγ(s)Δs+1ϕγ−1(δ(ℓ))∫ℓδ(ℓ)P(s)ψγ(s)ϕγ(δ(s))Δs +ϕ(δ(ℓ))∫∞ℓP(s)ψγ(s)Δs. | (2.16) |
This contradicts (2.7). Subsequently, let us assume that z∈ℵ0. Then limℓ→∞z(ℓ)=k≥0, and we propose that k=0. If not, it would imply z(ℓ)≥k>0. Integrating (2.2) from ℓ to ∞ yields
r(ℓ)(a(ℓ)zΔ(ℓ))Δ≥∫∞ℓP(s)zγ(δ(s))Δs≥kγ∫∞ℓP(s)Δs. |
Therefore,
−a(ℓ)zΔ(ℓ)≥kγ∫∞ℓ1r(u)∫∞uP(s)ΔsΔu, |
and
z(ℓ1)≥kγ∫∞ℓ11a(u)∫∞u1r(v)∫∞vP(s)ΔsΔvΔu. |
This leads to a contradiction to (2.8). Thus, we conclude: limℓ→∞z(ℓ)=limℓ→∞x(ℓ)A(ℓ)=0, and, the proof of the theorem is complete.
Theorem 2.3. Let 0<γ<1 and (2.1) hold. If (2.8) and
lim supℓ→∞{1ϕ(δ(ℓ))∫δ(ℓ)ℓ1ϕ(σ(s))P(s)ψγ(s)Δs+∫ℓδ(ℓ)P(s)ψγ(s)ϕγ(δ(s))Δs +ϕγ(δ(ℓ))∫∞ℓP(s)ψγ(s)Δs}=∞ | (2.17) |
hold, then every non-oscillatory solutionz(ℓ) of (1.4) satisfies limℓ→∞x(ℓ)A(ℓ)=0.
Proof. Let x(ℓ) be a non-oscillatory solution of Eq (1.4), where x(ℓ)>0, and x(δ(ℓ))>0 for ℓ≥ℓ1 for some ℓ1≥ℓ0. According to Theorem 2.1, the corresponding function z(ℓ)=x(ℓ)A(ℓ) is a positive solution of (2.2), implying that either z∈ℵ0 or z∈ℵ2 for ℓ≥ℓ1.
First, let us assume that z∈ℵ2. Proceeding similarly to the proof of Theorem 2.7, we arrive at (2.15). Dividing (2.16) by ϕ1−γ(δ(ℓ), we obtain
(χ(δ(ℓ))ϕ(δ(ℓ)))1−γ≥1ϕ(δ(ℓ))∫δ(ℓ)ℓ1ϕ(σ(s))P(s)ψγ(s)Δs+∫ℓδ(ℓ)P(s)ψγ(s)ϕγ(δ(s))Δs +ϕγ(δ(ℓ))∫∞ℓP(s)ψγ(s)Δs. | (2.18) |
In view of the decreasing nature of χ(δ(ℓ))/ϕ(δ(ℓ)) and the fact that 0<γ<1, there exists a constant C>0 such that
(χ(δ(ℓ))ϕ(δ(ℓ)))1−γ≤C. |
Taking the lim sup as ℓ→∞, we establish a contradiction to (2.18), and consequently, z∉ℵ2.
Subsequently, let us assume that z∈ℵ0. Proceeding similarly to the proof of Theorem 2.7, it becomes evident that condition (2.8) once more leads to the conclusion that limℓ→∞x(ℓ)A(ℓ)=0. This completes the proof.
Theorem 2.4. Suppose that conditions (i), (ii), and δΔ(ℓ)>0 are satisfied on [ℓ0,∞)T, γ≤1, and there exists a function ξ(ℓ) such that
ξΔ(ℓ)≥0, ξ(ℓ)>ℓ ,and θ(ℓ)=δ(ξ(ξ(ℓ)))<ℓ. | (2.19) |
If
lim infℓ→∞∫ℓδ(ℓ)P(s)ψγ(δ(s))Δs{=∞,γ<1,>1/e,γ=1, | (2.20) |
and
lim infℓ→∞∫ℓθ(ℓ)(1a(s)∫ξ(ℓ)s1r(u)∫ξ(u)uP(v)ΔvΔu)Δs{=∞,γ<1,>1/e,γ=1, | (2.21) |
for all ℓ1≥ℓ0, then Eq (1.4) is oscillatory.
Proof. Let x(ℓ) be a non-oscillatory solution of Eq (1.4), where x(ℓ)>0, and x(δ(ℓ))>0 for ℓ≥ℓ1 for some ℓ1≥ℓ0. According to Theorem 2.1, the corresponding function z(ℓ)=x(ℓ)A(ℓ) is a positive solution of (2.2), implying that either z∈ℵ0 or z∈ℵ2 for ℓ≥ℓ1. Assuming that z(ℓ)∈ℵ2, we have
a(ℓ)zΔ(ℓ)≥∫ℓℓ1r−1(s)r(s)(a(s)zΔ(s))ΔΔs≥r(ℓ)(a(ℓ)zΔ(ℓ))Δϕ(ℓ). |
It follows that
zΔ(ℓ)≥r(ℓ)(a(ℓ)zΔ(ℓ))Δϕ(ℓ)a(ℓ). | (2.22) |
Integrating the above inequality from ℓ2 to ℓ, we obtain
z(ℓ)≥∫ℓℓ2r(s)(a(s)zΔ(s))Δϕ(s)a(s)Δs≥r(ℓ)(a(ℓ)zΔ(ℓ))Δ∫ℓℓ2ϕ(s)a(s)Δs=r(ℓ)(a(ℓ)zΔ(ℓ))Δψ(ℓ). | (2.23) |
There exists ℓ3≥ℓ2 such that δ(ℓ)≥ℓ2 for all ℓ≥ℓ3. Then, we have
z(δ(ℓ))≥r(δ(ℓ))(a(δ(ℓ))zΔ(δ(ℓ)))Δψ(δ(ℓ)), for all ℓ≥ℓ3. |
Combining this with (2.2) yields
YΔ(ℓ)+P(ℓ)ψγ(δ(ℓ))Yγ(δ(ℓ))≤0, for ℓ≥ℓ3, | (2.24) |
where Y(ℓ):=r(ℓ)(a(ℓ)zΔ(ℓ))Δ. Integrating (2.24) from δ(ℓ) to ℓ, we have
Y(δ(ℓ))≥Y(δ(ℓ))−Y(ℓ)≥Yγ(δ(ℓ))∫ℓδ(ℓ)P(s)ψγ(δ(s))Δs. | (2.25) |
Hence,
Y1−γ(δ(ℓ))≥∫ℓδ(ℓ)P(s)ψγ(δ(s))Δs for ℓ≥ℓ3. |
According to [29, Theorem 1], we reach the intended contradiction.
Now, consider z∈ℵ0. Integrating (2.2) from ℓ to ξ(ℓ), we obtain
r(ℓ)(a(ℓ)zΔ(ℓ))Δ≥∫ξ(ℓ)ℓP(s)zγ(δ(s))Δs≥zγ(δ(ξ(ℓ)))∫ξ(ℓ)ℓP(s)Δs, |
where θ(ℓ):=δ(ξ(ξ(ℓ))). Consequently,
(a(ℓ)zΔ(ℓ))Δ≥zγ(δ(ξ(ℓ)))r(ℓ)∫ξ(ℓ)ℓP(s)Δs. | (2.26) |
Integrating (2.26) from ℓ to ξ(ℓ), we have
−a(ℓ)zΔ(ℓ)∫ξ(ℓ)ℓzγ(δ(ξ(s)))r(s)∫ξ(s)sP(u)ΔuΔs≥zγ(δ(ξ(ξ(ℓ))))∫ξ(ℓ)ℓ1r(s)∫ξ(s)sP(u)ΔuΔs=zγ(θ(ℓ))∫ξ(ℓ)ℓ1r(s)∫ξ(s)sP(u)ΔuΔs. | (2.27) |
It follows that
zΔ(ℓ)+(1a(ℓ)∫ξ(ℓ)ℓ1r(s)∫ξ(s)sP(u)ΔuΔs)zγ(θ(ℓ))≤0. | (2.28) |
The remainder of the proof follows a similar pattern to the one described above and is therefore omitted.
Theorem 2.5. Let (2.1) hold. Assume that there exists a function ρ(ℓ)∈C1rd(T,R+), such that
lim supℓ→∞∫ℓℓ0(P(s)ρ(s)ψ(δ(s))ϕ(s)λγ−1−ρΔ(s)r(s)4ρ(s))Δs=∞, | (2.29) |
and (2.8) hold. Then every solution z(ℓ) of (1.4) is oscillatory or satisfies limℓ→∞x(ℓ)A(ℓ)=0.
Proof. Let x(ℓ) be a non-oscillatory solution of Eq (1.4), where x(ℓ)>0, and x(δ(ℓ))>0 for ℓ≥ℓ1 for some ℓ1≥ℓ0. According to Theorem 2.1, the corresponding function z(ℓ)=x(ℓ)A(ℓ) is a positive solution of (2.2), implying that either z∈ℵ0 or z∈ℵ2 for ℓ≥ℓ1.
Firstly, let us consider z∈ℵ2; then we have r(ℓ)(a(ℓ)zΔ(ℓ))Δ is decreasing, and moreover,
r(ℓ)(a(ℓ)zΔ(ℓ))Δ≥∫∞ℓP(s)zγ(δ(s))Δs≥zγ(ℓ(s))∫∞ℓP(s)Δs. | (2.30) |
Let us define the generalized Riccati substitution
ω(ℓ)=ρ(ℓ)r(ℓ)(a(ℓ)zΔ(ℓ))Δa(ℓ)zΔ(ℓ). | (2.31) |
Applying both the product rule and the quotient rule, we obtain
ωΔ(ℓ)=(r(ℓ)(a(ℓ)zΔ(ℓ))Δ)Δ(ρ(ℓ)a(ℓ)zΔ(ℓ))+(r(ℓ)(a(ℓ)zΔ(ℓ))Δ)σ(ρ(ℓ)a(ℓ)zΔ(ℓ))Δ=(r(ℓ)(a(ℓ)zΔ(ℓ))Δ)Δ(ρ(ℓ)a(ℓ)zΔ(ℓ)) +(r(ℓ)(a(ℓ)zΔ(ℓ))Δ)σ((a(ℓ)zΔ(ℓ))ρΔ(ℓ)−ρ(ℓ)(a(ℓ)zΔ(ℓ))Δ(a(ℓ)zΔ(ℓ))(a(ℓ)zΔ(ℓ))σ)≤−P(ℓ)ρ(ℓ)(zγ(δ(ℓ))a(ℓ)zΔ(ℓ))+ρΔ+(ℓ)ρ(σ(ℓ))ω(σ(ℓ)) −ρ(ℓ)(r(ℓ)(a(ℓ)zΔ(ℓ))Δ)σ(a(ℓ)zΔ(ℓ))Δ(a(ℓ)zΔ(ℓ))(a(ℓ)zΔ(ℓ))σ. | (2.32) |
Using the monotonicity of r(ℓ)(a(ℓ)rΔ(ℓ))Δ and (a(ℓ)rΔ(ℓ)), we have
(a(ℓ)zΔ(ℓ))Δ≥r(σ(ℓ))r(ℓ)(a(σ(ℓ))zΔ(σ(ℓ)))Δ, | (2.33) |
and
1a(ℓ)zΔ(ℓ)≥1a(σ(ℓ))zΔ(σ(ℓ)). | (2.34) |
Combining (2.33) and (2.34) with (2.32), we obtain
ωΔ(ℓ)≤−P(ℓ)ρ(ℓ)(zγ(δ(ℓ))a(ℓ)zΔ(ℓ))+ρΔ+(ℓ)ρ(σ(ℓ))ω(σ(ℓ))−ρ(ℓ)r(ℓ)ρ2(σ((ℓ))ω2(σ(ℓ)). | (2.35) |
From (2.9), (2.10), and the fact that δ(ℓ)≤ℓ, we have
z(δ(ℓ))a(ℓ)zΔ(ℓ)≥ψ(δ(ℓ))ϕ(ℓ), for ℓ≥ℓ3. | (2.36) |
This, together with (2.33), leads to
ωΔ(ℓ)≤−P(ℓ)ρ(ℓ)ψ(δ(ℓ))ϕ(ℓ)zγ−1(δ(ℓ))+ρΔ+(ℓ)ρ(σ(ℓ))ω(σ(ℓ))−ρ(ℓ)r(ℓ)ρ2(σ((ℓ))ω2(σ(ℓ)). | (2.37) |
Since zΔ(ℓ)>0, then there exists a constant λ>0 such that z(ℓ)≥λ for ℓ≥ℓ3. Consequently, (2.37) can be expressed as
ωΔ(ℓ)≤−P(ℓ)ρ(ℓ)ψ(δ(ℓ))ϕ(ℓ)λγ−1+ρΔ+(ℓ)ρ(σ(ℓ))ω(σ(ℓ))−ρ(ℓ)r(ℓ)ρ2(σ((ℓ))ω2(σ(ℓ))≤−P(ℓ)ρ(ℓ)ψ(δ(ℓ))ϕ(ℓ)λγ−1+ρΔ(ℓ)r(ℓ)4ρ(ℓ). | (2.38) |
Integrating both sides of (2.38) from ℓ4>ℓ3 to ℓ, we obtain
∫ℓℓ4(P(s)ρ(s)ψ(δ(s))ϕ(s)λγ−1−ρΔ(s)r(s)4ρ(s))Δs≤ω(ℓ4), | (2.39) |
which contradicts (2.29). Now, assume that z∈ℵ0. Proceeding similarly to the proof of Theorem 2.7, it becomes evident that condition (2.8) once more leads to the conclusion that limℓ→∞x(ℓ)A(ℓ)=0. This completes the proof.
Theorem 2.6. Let γ=1,
lim supℓ→∞∫ℓδ(ℓ)(1a(v)∫ℓv1r(u)∫ℓuP(s)ΔsΔu)Δv>1, | (2.40) |
and assume that there exists a function ρ(ℓ)∈C1rd(T,R+), such that
lim supℓ→∞∫ℓℓ0(P(s)ρ(s)ψ(δ(s))ϕ(s)−ρΔ(s)r(s)4ρ(s))Δs=∞. | (2.41) |
Then every solution of (1.4) is oscillatory.
Proof. Let x(ℓ) be a non-oscillatory solution of Eq (1.4), where x(ℓ)>0, and x(δ(ℓ))>0 for ℓ≥ℓ1 for some ℓ1≥ℓ0. According to Theorem 2.1, the corresponding function z(ℓ)=x(ℓ)A(ℓ) is a positive solution of (2.2), implying that either z∈ℵ0 or z∈ℵ2 for ℓ≥ℓ1. Assume z(ℓ)∈ℵ0. Integrating (2.2) from ν to ℓ yields
r(ν)(a(ν)zΔ(ν))Δ≥∫ℓνP(s)z(δ(s))Δs≥z(δ(ℓ))∫ℓνP(s)Δs. |
Integrating again twice from ν to ℓ, we obtain
z(ν)≥z(δ(ℓ))∫ℓν(1a(v)∫ℓv1r(u)∫ℓuP(s)ΔsΔu)Δv. |
Replacing ν with δ(ℓ) leads to contradiction to (2.40). Hence, every positive solution z(ℓ) does not satisfy ℵ0. Therefore, if (2.40) holds, then z(ℓ)∈ℵ2. Proceeding as in Theorem (2.29) with γ=1, completes the proof.
Example 2.1. Consider the third order linear differential equation
(1ℓ(ℓ2(x′(ℓ)))′)′+p0√ℓx(αℓ)=0, ℓ≥1, | (2.42) |
where p0 is a constant and α∈(0,1). Here a2(ℓ)=1ℓ, a1(ℓ)=ℓ2, p(ℓ)=p0√ℓ and δ(ℓ)=αℓ. It is clear that (2.42) is semi-canonical. Since A(ℓ)=1ℓ, a(ℓ)=r(ℓ)=1, and P(ℓ)=p0ℓ3/2, the corresponding canonical equation is
z′′′(ℓ)+p0ℓ3/2x(αℓ)=0. | (2.43) |
It is clear that (2.1) holds. Applying Theorem 2.6, we have
lim supℓ→∞∫ℓδ(ℓ)(1a(v)∫ℓv1r(u)∫ℓuP(s)ΔsΔu)Δv=lim supℓ→∞∫ℓαℓ(∫ℓv∫ℓup0s3/2dsdu)dv=limℓ→∞(−α2+(1−2α)√α+4)ℓ3/2>1 |
and by choosing ρ(ℓ)=ℓ
lim supℓ→∞∫ℓℓ0(P(s)ρ(s)ψ(δ(s))ϕ(s)−ρΔ(s)r(s)4ρ(s))Δs=lim supℓ→∞∫ℓℓ0(p0s3/2α2s3s−14s)ds=∞. |
It follows that (2.42) is oscillatory. Also, by Theorems 2.2 and 2.4, Eq (2.42) is oscillatory or limℓ→∞x(ℓ)A(ℓ)=0.
Remark 2.1. It is worth noting that the existing results in [10,30,31,32] cannot be directly applied to Eq (2.42) due to the fact that a1(ℓ)≠1.
Example 2.2. Consider the second order difference equation
Δ(1ℓ+1Δ(ℓ(ℓ+1))Δx(ℓ))+p0x1/2(ℓ−2)=0, ℓ≥1, | (2.44) |
where p0 is a constant. Here a2(ℓ)=1ℓ+1, a1(ℓ)=ℓ(ℓ+1), p(ℓ)=p0 and δ(ℓ)=ℓ−2. It is clear that (2.44) is semi-canonical. Since A(ℓ)=1ℓ, a(ℓ)=r(ℓ)=1, and P(ℓ)=p0ℓ1/2, the corresponding canonical equation is
Δ(Δ(Δ(z(ℓ))))+p0ℓ1/2x(ℓ−2)=0. | (2.45) |
It is clear that (2.1) and (2.8) hold. Further, (2.17) becomes
lim supℓ→∞{1ℓ−3ℓ−2∑1(s+1)p0s1/2s2+ℓ∑ℓ−2sp0√2s(s−2) +√(ℓ−2)∞∑ℓp0s√2s}=∞. |
Hence, by Theorem 2.3, every solution is oscillatory or limℓ→∞x(ℓ)A(ℓ)=0.
The results of this study are presented in a novel and generalizable framework, highlighting their broad applicability. Our approach involves a unique transformation that converts the equation from the semi-canonical form to the more tractable canonical form. This transformation facilitates the derivation of new oscillation criteria with fewer restrictions compared to the existing literature. Theorems 2.19 and 2.40 illustrate our criteria, ensuring that all solutions oscillate. The results obtained are consistent with the results in [11,13,14] and can be extended to non linear difference equations. Our approach has the potential to be extended to both non-canonical and semi-canonical forms (as defined in (S2)), potentially leading to new oscillation conditions.
Ahmed M. Hassan: Writing-original draft, Writing-review and editing, Making major revisions; Clemente Cesarano: Supervision, Writing-review and editing; Sameh S. Askar: Formal analysis, Writing-original draft; Ahmad M. Alshamrani: Writing-original draft, Making major revisions. All authors have read and approved the final version of the manuscript for publication.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors present their appreciation to King Saud University for funding this research through Researchers Supporting Project number (RSPD2024R533), King Saud University, Riyadh, Saudi Arabia.
The authors declare that there is no conflict of interest regarding the publication of this paper.
[1] | S. Hilger, Ein maßkettenkalkül mit anwendung auf zentrumsmannigfaltigkeiten, PhD thesis, 1988. |
[2] |
S. Hilger, Analysis on measure chains-a unified approach to continuous and discrete calculus, Results Math., 18 (1990), 18–56. https://doi.org/10.1007/BF03323153 doi: 10.1007/BF03323153
![]() |
[3] | M. Bohner, A. Peterson, Dynamic equations on time scales: An introduction with applications, Springer Science & Business Media, 2001. https://doi.org/10.1007/978-1-4612-0201-1 |
[4] |
G. S. Guseinov, B. Kaymakçalan, On a disconjugacy criterion for second order dynamic equations on time scales, J. Comput. Appl. Math., 141 (2002), 187–196. https://doi.org/10.1016/S0377-0427(01)00445-9 doi: 10.1016/S0377-0427(01)00445-9
![]() |
[5] |
R. Agarwal, M. Bohner, D. O'Regan, A. Peterson, Advances in dynamic equations on time scales, J. Comput. Appl. Math., 141 (2002), 1–26. https://doi.org/10.1016/S0377-0427(01)00432-0 doi: 10.1016/S0377-0427(01)00432-0
![]() |
[6] |
R. P. Agarwal, M. Bohner, Basic calculus on time scales and some of its applications, Results Math., 35 (1999), 3–22. https://doi.org/10.1007/BF03322019 doi: 10.1007/BF03322019
![]() |
[7] | M. Bohner, A. Peterson, Advances in dynamic equations on time scales, Boston: Birkhäuser, 2003. https://doi.org/10.1007/978-0-8176-8230-9 |
[8] |
R. P. Agarwal, M. Bohner, T. Li, C. Zhang, Hille and Nehari type criteria for third-order delay dynamic equations, J. Differ. Equ. Appl., 19 (2013), 1563–1579. https://doi.org/10.1080/10236198.2013.766729 doi: 10.1080/10236198.2013.766729
![]() |
[9] |
T. S. Hassan, R. Agarwal, W. Mohammed, Oscillation criteria for third-order functional half-linear dynamic equations, Adv. Differ. Equ., 2017 (2017), 1–28. https://doi.org/10.1186/s13662-016-1057-2 doi: 10.1186/s13662-016-1057-2
![]() |
[10] |
G. Chatzarakis, J. Džurina, I. Jadlovská, Oscillatory and asymptotic properties of third-order quasilinear delay differential equations, J. Inequal. Appl., 2019 (2019), 1–17. https://doi.org/10.1186/s13660-019-1967-0 doi: 10.1186/s13660-019-1967-0
![]() |
[11] |
G. E. Chatzarakis, S. R. Grace, I. Jadlovská, Oscillation criteria for third-order delay differential equations, Adv. Differ. Equ., 2017 (2017), 1–11. https://doi.org/10.1186/s13662-016-1057-2 doi: 10.1186/s13662-016-1057-2
![]() |
[12] | L. Erbe, B. Karpuz, A. Peterson, Kamenev-type oscillation criteria for higher-order neutral delay dynamic equations, Int. J. Differ. Equ., 6 (2011), 1–16. |
[13] |
L. Erbe, A. Peterson, S. H. Saker, Asymptotic behavior of solutions of a third-order nonlinear dynamic equation on time scales, J. Comput. Appl. Math., 181 (2005), 92–102. https://doi.org/10.1016/j.cam.2004.11.021 doi: 10.1016/j.cam.2004.11.021
![]() |
[14] |
T. S. Hassan, Oscillation of third order nonlinear delay dynamic equations on time scales, Math. Comput. Model., 49 (2009), 1573–1586. https://doi.org/10.1016/j.mcm.2008.12.011 doi: 10.1016/j.mcm.2008.12.011
![]() |
[15] |
Y. Wang, Z. Xu, Asymptotic properties of solutions of certain third-order dynamic equations, J. Comput. Appl. Math., 236 (2012), 2354–2366. https://doi.org/10.1016/j.cam.2011.11.021 doi: 10.1016/j.cam.2011.11.021
![]() |
[16] |
S. Salem, A. M. Hassan, Oscillatory behavior of solutions of third-order nonlinear neutral delay dynamic equations on time scales, Mediterr. J. Math., 20 (2023), 308. https://doi.org/10.1007/s00009-023-02506-y doi: 10.1007/s00009-023-02506-y
![]() |
[17] |
F. Masood, O. Moaaz, G. AlNemer, H. El-Metwally, More effective criteria for testing the asymptotic and oscillatory behavior of solutions of a class of third-order functional differential equations, Axioms, 12 (2023), 1112. https://doi.org/10.3390/axioms12121112 doi: 10.3390/axioms12121112
![]() |
[18] |
L. Gao, S. Liu, X. Zheng, New oscillatory theorems for third-order nonlinear delay dynamic equations on time scales, J. Appl. Math. Phys., 6 (2018), 232. https://doi.org/10.4236/jamp.2018.61023 doi: 10.4236/jamp.2018.61023
![]() |
[19] | H. Wu, L. Erbe, A. Peterson, Oscillation of solution to second-order half-linear delay dynamic equations on time scales, 2016. |
[20] |
Y. Yu, Q. Wang, Q. Bi, C. Lim, Multiple-s-shaped critical manifold and jump phenomena in low frequency forced vibration with amplitude modulation, Int. J. Bifurcat. Chaos, 29 (2019), 1930012. https://doi.org/10.1142/S021812741930012X doi: 10.1142/S021812741930012X
![]() |
[21] |
O. Moaaz, B. Qaraad, R. A. El-Nabulsi, O. Bazighifan, New results for kneser solutions of third-order nonlinear neutral differential equations, Mathematics, 8 (2020), 686. https://doi.org/10.3390/math8050686 doi: 10.3390/math8050686
![]() |
[22] |
O. Moaaz, D. Chalishajar, O. Bazighifan, Asymptotic behavior of solutions of the third order nonlinear mixed type neutral differential equations, Mathematics, 8 (2020), 485. https://doi.org/10.3390/math8040485 doi: 10.3390/math8040485
![]() |
[23] |
O. Moaaz, R. A. El-Nabulsi, W. Muhsin, O. Bazighifan, Improved oscillation criteria for 2nd-order neutral differential equations with distributed deviating arguments, Mathematics, 8 (2020), 849. https://doi.org/10.3390/math8050849 doi: 10.3390/math8050849
![]() |
[24] |
S. S. Santra, A. K. Sethi, O. Moaaz, K. M. Khedher, S. W. Yao, New oscillation theorems for second-order differential equations with canonical and non-canonical operator via riccati transformation, Mathematics, 9 (2021), 1111. https://doi.org/10.3390/math9101111 doi: 10.3390/math9101111
![]() |
[25] |
A. Muhib, T. Abdeljawad, O. Moaaz, E. M. Elabbasy, Oscillatory properties of odd-order delay differential equations with distribution deviating arguments, Appl. Sci., 10 (2020), 5952. https://doi.org/10.3390/app10175952 doi: 10.3390/app10175952
![]() |
[26] |
O. Moaaz, J. Awrejcewicz, A. Muhib, Establishing new criteria for oscillation of odd-order nonlinear differential equations, Mathematics, 8 (2020), 937. https://doi.org/10.3390/math8060937 doi: 10.3390/math8060937
![]() |
[27] |
J. R. Graef, Canonical, noncanonical, and semicanonical third order dynamic equations on time scales, Results Nonlinear Anal., 5 (2022), 273–278. https://doi.org/10.53006/rna.1075859 doi: 10.53006/rna.1075859
![]() |
[28] |
W. F. Trench, Canonical forms and principal systems for general disconjugate equations, T. Am. Math. Soc., 189 (1974), 319–327. https://doi.org/10.1090/S0002-9947-1974-0330632-X doi: 10.1090/S0002-9947-1974-0330632-X
![]() |
[29] | R. Agarwal, M. Bohner, An oscillation criterion for first order delay dynamic equations, Funct. Differ. Equ., 16 (2004), 11–17. |
[30] | R. Agarwal, M. Bohner, T. Li, C. Zhang, Oscillation of third-order nonlinear delay differential equations, Taiwan. J. Math., 17 (2013). https://doi.org/10.11650/tjm.17.2013.2095 |
[31] | B. Baculíková, J. Džurina, Oscillation of third-order neutral differential equations, Math. Comput. Model., 35 (2010), 215–226. |
[32] | S. R. Grace, R. P. Agarwal, M. F. Aktas, On the oscillation of third order functional differential equations, Indian J. Pure Appl. Math., 39 (2008), 491–507. |