407 Β© 2025 The Author(s). Published by College of Education for Pure Science (Ibn Al-Haitham), University of Baghdad. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License Oscillation Criteria for Solutions of Neutral Differential Equations of the Second-Order Emden-Fowler Type with Forcing Term Jihan Saad1* and Hussain Ali Mohamad2 1,2Department of Mathematics, College of Science for Women, University of Baghdad, Baghdad, Iraq. *Corresponding Author. Received:22 May 2023 Accepted:14 August 2023 Published:20 January 2025 doi.org/10.30526/38.1.3515 Abstract In this paper, the oscillation property and the asymptotic behavior of solutions of neutral second- order differential Equations of the Emden-Fowler type were studied under the influence of the coefficients of forces. It has been shown through this research that the coefficients of forces in addition to the Emden-Fowler type have a major role on the oscillation of solutions of neutral Equations. As well as its effect on the convergence and divergence of nonoscillatory solutions. For this purpose, some conditions are obtained to ensure that all solutions of the neutral Equations Emden-Fowler type oscillating or nonoscillating go to ∞, as t β†’ ∞. Some of these conditions are the development of conditions similar to them in some of the well-known results included in the references, for example, condition (8) in this research with condition (4) in (9). The obtained results included some illustrative examples showing that the resulting conditions are easy to apply and guarantee oscillation. Keywords: Oscillation Criteria, Asymptotic Behavior, Emden-Fowler Type, Neutral Second Order with Forcing Term. 1. Introduction This paper aims to obtain sufficient conditions to ensure that every solution of the neutral force Equation of the second-order type Emden Fowler oscillates. Consider the Equation: (𝝃(𝒕)(πŽβ€²(𝒕)) 𝜸 ) β€² + βˆ‘ π’’π’Š(𝒕)π’™πœΈ(πœΉπ’Š(𝒕)) 𝒏 π’Š=𝟏 𝐬𝐠𝐧(𝒙) = βˆ‘ 𝒓𝒋(𝒕) π’Œ 𝒋=𝟏 . (𝟏) πœ”(𝑑) = π‘₯(𝑑) + 𝑝(𝑑)π‘₯(𝜏(𝑑)). (2) https://creativecommons.org/licenses/by/4.0/ https://creativecommons.org/licenses/by/4.0/ https://doi.org/10.30526/38.1.3501 https://orcid.org/0009-0000-2264-7057 mailto:jihan.saad2103m@csw.uobaghdad.edu.iq https://orcid.org/0000-0002-4684-786x mailto:hussainam_math@csw.uobaghdad.edu.iq IHJPAS. 2024, 38 (1) 408 The number 𝛾 is a quotient of odd positive integers. 𝜏, 𝛿𝑖 ∈ ∁( [𝑑0, ∞)𝕋 , 𝑅), 𝑖 = 1,2, … , 𝑛, lim tβ†’βˆž 𝜏(𝑑) = ∞ , lim tβ†’βˆž 𝛿𝑖(𝑑) = ∞ , sgn(π‘₯) = Β±1 if π‘₯ β‰· 0, 𝑠𝑔𝑛(0) = 0, 𝑝 ∈ ∁([𝑑0, ∞) , 𝑅+) and π‘žπ‘–, π‘Ÿπ‘— ∈ ∁( [𝑑0, ∞) , 𝑅), 𝑖 = 1,2, … , 𝑛, 𝑗 = 1,2, … , π‘˜. During this research, the following assumptions will be used as needed: (M1) lim sup π‘‘β†’βˆž ∫ ( 1 πœ‰(𝑠) ) 1 𝛾 𝑑𝑠 = ∞; 𝑑 𝑑0 (M2) 0 < 𝑝(𝑑) ≀ π‘Ž, πœ‰(𝑑) > 0 𝑄𝑖(𝑑) = min 𝑑β‰₯𝑑0 { π‘žπ‘–(𝑑), π‘žπ‘–(𝜏(𝑑)), 𝑄(𝑑) = min{𝑄𝑖(𝑑), 𝑖 = 1,2, … , 𝑛} , and 𝐺𝑖(𝑑) = max 𝑑β‰₯𝑑0 { π‘Ÿπ‘–(𝑑), π‘Ÿπ‘–(𝜏(𝑑)), 𝑖 = 1,2, … , π‘˜} , 𝐺(𝑑) = max{𝐺𝑖(𝑑), 𝑖 = 1,2, … , π‘˜}. (M3) ∫ |𝐺(𝑑)|𝑑𝑑 ∞ 𝑇 < ∞, 𝑇 β‰₯ 𝑑0. (M4) ∫ |𝑄(𝑑)|𝑑𝑑 ∞ 𝑇 = ∞, 𝑇 β‰₯ 𝑑0 . The Emden-Fowler Equation has emerged in recent decades as a focus of interest for many researchers, specifically research in oscillation and the asymptotic behavior of the solutions of these Equations, Emden-Fowler Equation has been classified as unconventional Equations and its importance has emerged for its use in many applications, and for this reason many researches have appeared that produce a lot of conditions to ensure that each solution of these Equations oscillates, or that their non-oscillating solutions are convergent. Ahmed et al. (1) studied the second-order neutral dynamic linear Equation and established some conditions for the oscillation of every solution of this Equation. Yingzhu et al. (2) obtained oscillation conditions of each solution of second-order neutral nonlinear differential Equations. The obtained results in (3) are based on comparison theorems which enable to address the problem of second order Equation oscillation to first order Equation oscillation. Mehta et al. (4) investigated the Emden-Fowler Equation of the form 𝑑 𝑑𝑑 (𝑑𝛼 𝑑𝑀 𝑑𝑑 ) = π‘‘πœŽπ‘€π›Ό, where 𝑀(𝑑) = π‘₯(𝑑) Β± π‘ž(𝑑)π‘₯(𝜏(𝑑)), and established some conditions for all solutions to oscillate. Mohamad et al. (5, 6) discussed the oscillation property of third order neutral half-linear Equations and established some conditions to insure the oscillation of every solution of these Equations. Moaaz et al. (7) and Xu et al. (8) obtained oscillation conditions of each solution of second order neutral Emden-Fowler type (π‘Ž(𝑑)[(π‘₯(𝑑) βˆ’ 𝑝(𝑑)π‘₯(𝜏(𝑑)))β€²]𝛾)β€² + π‘ž(𝑑)π‘₯𝛾(𝜎(𝑑)) = 0, t β‰₯ 𝑑0. Thandapani et al. (9) studied asymptotic properties of the third order quasi-linear neutral functional differential Equation (π‘Ž(𝑑)[(π‘₯(𝑑) βˆ’ 𝑝(𝑑)π‘₯(𝜏(𝑑)))β€²β€²]𝛾)β€² + π‘ž(𝑑)π‘₯𝛾(𝜎(𝑑)) = 0, by using the Riccati transformation, and establishing some conditions which ensure that every solution of that Equation is either oscillatory or converges to zero. Hassan et al. (10) he dealt with new standards for the oscillation of half-linear differential Equations developed of the second order, where he concluded that the results obtained work to expand and develop modern standards for the same Equations that have been developed by many authors. Tripathy et al. (11) find the necessary and sufficient conditions for volatility one of the impulsive neutral differential system solutions of the second order under certain conditions that ensure the occurrence of oscillation . See Mehta et al. (12), and Vidhyaa et al. (13) they studied the differential Equations and obtained the oscillation criteria for all the IHJPAS. 2024, 38 (1) 409 solutions of the neutral differential Equations of the second degree, half-linear (π‘˜(𝑑)((β„Ž(𝑑)𝑧′ (𝑑))β€²)πœ€)β€² + π‘˜(𝑑)π‘₯πœ€(𝑑) = 0, 𝑑 β‰₯ 𝑑0, where 𝑧(𝑑) = π‘₯(𝑑) + 𝑝(𝑑)π‘₯(𝜏(𝑑)). Dassios et al. (14) he studied the delayed and neutral differential Equations, where he focused on the stability of the important joins because these joins contain delays in each of the state variables and their time derivatives, where the proposed approach consists of model transformation that builds an equivalent set of algebraic differential Equations. Li et al. (15) and Marappan et al. (16) the oscillatory behavior of solutions of mixed nonlinear neutral differential Equations of the Emden-Fowler type was studied by applying the integral conditions and the integral average method. Baty (17) he studied second- order Lane-Emden-Fowler differential Equations, third-order Emden-Fowler Equations, and fourth- order Lane-Emden-Fowler Equations. He presented numerical methods using neural networks based on physics with the aim of solving higher-order differential Equations. Naeif et al. (18) he studied the oscillation and asymptotic behavior of a half-linear three-dimensional neutral system of second order and gave sufficient conditions to ensure oscillation or not. See (19-22) they studied the optimal decomposition method for solving third-order nonlinear Emden-Fowler differential Equations, to avoid the singularity at x=0, by transforming the Emden-Fowler Equation into an integral Volterra Equation. Our paper was based on article (9), where a more general Equation with a forcing term is used and condition (2) in (9) has been developed into a more general case. A solution 𝒙(𝒕) is said to be oscillatory if it has arbitrarily large zeros on (π’•πŸŽ, ∞), otherwise it is said to be nonoscillatory that is either eventually positive or eventually negative (6). 2. Main Results In this section some results .established for oscillation for every solution of Equation (1). In the beginning, it is shown that every non-oscillatory solution is achieves the following cases. Lemma 2.1: Assume that π‘žπ‘–(𝑑) β‰₯ 0, βˆ‘ π‘Ÿπ‘—(𝑑)π‘˜ 𝑗=1 ≀ 0, and (M1) holds. Let π‘₯(𝑑) be a non- oscillatory solution of Equation (1). Then πœ”β€²(𝑑) > 0, and either lim π‘‘β†’βˆž π‘₯(𝑑) = ∞, or lim π‘‘β†’βˆž πœ‰(𝑑)(πœ”β€²(𝑑)) 𝛾 = 0. Proof: Assume that π‘₯(𝑑) be eventually positive solution of Equation (1). From (1) it follows that (πœ‰(𝑑)(πœ”β€²(𝑑))𝛾)β€² ≀ 0, that is πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 is non-increasing for 𝑑 β‰₯ 𝑑0. we claim that πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 is eventually positive, otherwise if πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 is eventually negative then there is πœ‡ < 0 and 𝑑1 β‰₯ 𝑑0 such that πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 ≀ πœ‡ < 0, 𝑑 β‰₯ 𝑑1 , so it follows πœ”β€²(𝑑) ≀ ( πœ‡ πœ‰(𝑑) ) 1 𝛾 , 𝑑 β‰₯ 𝑑1 . (3) Integrating (3) from 𝑑1 to 𝑑 we get πœ”(𝑑) βˆ’ πœ”(𝑑1) ≀ πœ‡ 1 𝛾 ∫ ( 1 πœ‰(𝑠) ) 1 𝛾𝑑 𝑑1 𝑑𝑠. (4) Letting 𝑑 β†’ ∞, then from inequality (4) yields lim π‘‘β†’βˆž πœ”(𝑑) = βˆ’βˆž, a contradiction. Hence our claim verified and πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 > 0, that is πœ”β€²(𝑑) > 0, 𝑑 β‰₯ 𝑑1 β‰₯ 𝑑0, then lim π‘‘β†’βˆž πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 = 𝑙 β‰₯ 0, thus πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 β‰₯ 𝑙, 𝑑 β‰₯ 𝑑1 or πœ”β€²(𝑑) β‰₯ ( 𝑙 πœ‰(𝑑) ) 1 𝛾 , 𝑑 β‰₯ 𝑑1 . (5) IHJPAS. 2024, 38 (1) 410 Integrating (5) from 𝑑1 to 𝑑, it follows πœ”(𝑑) βˆ’ πœ”(𝑑1) β‰₯ 𝑙 1 𝛾 ∫ ( 1 πœ‰(𝑠) ) 1 𝛾𝑑 𝑑1 𝑑𝑠. If 𝑙 > 0 then lim π‘‘β†’βˆž πœ”(𝑑) = ∞, implies that lim π‘‘β†’βˆž π‘₯(𝑑) = ∞. If 𝑙 = 0 then lim π‘‘β†’βˆž πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 = 0. Lemma 2.2: Assume that π‘žπ‘–(𝑑) ≀ 0, βˆ‘ π‘Ÿπ‘—(𝑑)π‘˜ 𝑗=1 β‰₯ 0, and (M1) holds. Let π‘₯(𝑑) be a non- oscillatory solution of Equation (1). Then the following statements hold: (a) πœ”β€²(𝑑) > 0, and lim π‘‘β†’βˆž π‘₯(𝑑) = ∞. (b) πœ”β€²(𝑑) < 0, and lim π‘‘β†’βˆž πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 = 0. Proof: Assume that π‘₯(𝑑) be eventually positive solution of Equation (1). From Equation (1) we get (πœ‰(𝑑)(πœ”β€²(𝑑))𝛾)β€² β‰₯ 0 that is πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 is non-decreasing, we have two cases to consider: 1. πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 > 0; 2. πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 < 0, 𝑑 β‰₯ 𝑑1 β‰₯ 𝑑0. Case 1: πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 > 0 that is πœ”β€²(𝑑) > 0, 𝑑 β‰₯ 𝑑1 then there exist πœ‡ > 0 such that πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 β‰₯ πœ‡, 𝑑 β‰₯ 𝑑2 β‰₯ 𝑑1 πœ”β€²(𝑑) β‰₯ ( πœ‡ πœ‰(𝑑) ) 1 𝛾 , 𝑑 β‰₯ 𝑑2. (6) By integrating (6) from 𝑑2 to 𝑑 we get πœ”(𝑑) βˆ’ πœ”(𝑑2) β‰₯ πœ‡ 1 𝛾 ∫ ( 1 πœ‰(𝑠) ) 1 𝛾𝑑 𝑑2 𝑑𝑠. As 𝑑 β†’ ∞ it follows lim π‘‘β†’βˆž πœ”(𝑑) = ∞, which implies that lim π‘‘β†’βˆž π‘₯(𝑑) = ∞. Case 2: πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 < 0 that is πœ”β€²(𝑑) < 0, 𝑑 β‰₯ 𝑑1 and lim π‘‘β†’βˆž πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 = 𝑙 ≀ 0. we claim that 𝑙 = 0 otherwise 𝑙 < 0 thus πœ”β€²(𝑑) ≀ ( 𝑙 πœ‰(𝑑) ) 1 𝛾 , 𝑑 β‰₯ 𝑑2 β‰₯ 𝑑1 . (7) By integrating Equation (7) from 𝑑2 to t we get πœ”(𝑑) βˆ’ πœ”(𝑑2) ≀ 𝑙 1 𝛾 ∫ ( 1 πœ‰(𝑠) ) 1 𝛾𝑑 𝑑2 𝑑𝑠. As 𝑑 β†’ ∞ it follows that lim π‘‘β†’βˆž πœ”(𝑑) = βˆ’βˆž, a contradiction. Then 𝑙 = 0. Theorem 2.1. Assume that π‘žπ‘–(𝑑) β‰₯ 0, βˆ‘ π‘Ÿπ‘—(𝑑)π‘˜ 𝑗=1 ≀ 0, (𝑀1) βˆ’ (𝑀4) hold, and for any continuous functions 𝑒(𝑑), 𝑣(𝑑), 𝑒𝑣 > 0, there exists πœ† > 0, such that 𝑒𝛾(𝑑) + 𝑣𝛾(𝑑) β‰₯ πœ†(𝑒(𝑑) + 𝑣(𝑑)) 𝛾 . (8) Then every solution of Equation (1) oscillates. Proof. Assume that Equation (1) has a non-oscillatory solution π‘₯(𝑑). For lack of prolongation and repetition, it can be assumed that π‘₯(𝑑) > 0, π‘₯(𝜏(𝑑)) > 0, π‘₯(𝛿𝑖(𝑑)) > 0, 𝑖 = 1,2, … , 𝑛, for 𝑑 β‰₯ 𝑑0 . Let 𝑒(𝑑) = π‘₯(𝑑), 𝑣(𝑑) = 𝑝(𝑑) π‘₯(𝜏(𝑑)), πœ‚(𝑑) = max {π‘₯(𝑑), π‘₯(𝜏(𝑑))}, then (π‘₯(𝑑) + 𝑝(𝑑)π‘₯(𝜏(𝑑)))𝛾 ≀ (π‘₯(𝑑) + π‘Ž π‘₯(𝜏(𝑑)))𝛾 ≀ (πœ‚(𝑑) + π‘Ž πœ‚(𝑑))𝛾 = πœ‚π›Ύ(𝑑)(1 + π‘Ž)𝛾. Since πœ‚(𝑑) = max {π‘₯(𝑑), π‘₯(𝜏(𝑑))} so there exists νœ€ > 0, such that π‘₯𝛾(𝑑) + 𝑝𝛾(𝑑)π‘₯𝛾(𝜏(𝑑)) β‰₯ νœ€πœ‚π›Ύ(𝑑) = πœ€(1+π‘Ž)𝛾 (1+π‘Ž)𝛾 πœ‚π›Ύ(𝑑) = πœ€ (1+π‘Ž)𝛾 (πœ‚(𝑑) + π‘Ž πœ‚(𝑑))𝛾 β‰₯ πœ€ (1+π‘Ž)𝛾 (π‘₯(𝑑) + 𝑝(𝑑)π‘₯(𝜏(𝑑)))𝛾, Choose πœ† = πœ€ (1+𝒢)𝛾 > 0. Hence, IHJPAS. 2024, 38 (1) 411 π‘₯𝛾(𝑑) + 𝑝𝛾(𝑑)π‘₯𝛾(𝜏(𝑑)) β‰₯ πœ†(π‘₯(𝑑) + 𝑝(𝑑)π‘₯(𝜏(𝑑)))𝛾 = πœ†πœ”π›Ύ(𝑑). (9) From Equation (1) it follows that: (πœ‰(𝑑)(πœ”β€²(𝑑))𝛾)β€² + βˆ‘ π‘žπ‘–(𝑑)(π‘₯(𝛿𝑖(𝑑)) 𝛾𝑛 𝑖=1 + π‘Žπ›Ύ[πœ‰(𝜏(𝑑))(πœ”β€²(𝜏(𝑑)))𝛾]β€² + π‘Žπ›Ύ βˆ‘ π‘žπ‘–(𝜏(𝑑))(π‘₯ (𝛿𝑖(𝜏(𝑑))) 𝛾𝑛 𝑖=1 = βˆ‘ (π‘Ÿπ‘—(𝑑) π‘˜ 𝑗=1 + π‘Žπ›Ύπ‘Ÿπ‘—(𝜏(𝑑)) (10) (πœ‰(𝑑)(πœ”β€²(𝑑))𝛾)β€² + π‘Žπ›Ύ[πœ‰(𝜏(𝑑))(πœ”β€²(𝜏(𝑑)))𝛾]β€² + 𝑄(𝑑) βˆ‘ [π‘₯𝛾(𝛿𝑖(𝑑) 𝑛 𝑖=1 + π‘Žπ›Ύπ‘₯𝛾(𝛿𝑖(𝜏(𝑑))] βˆ’ βˆ‘ (π‘Ÿπ‘—(𝑑) π‘˜ 𝑗=1 + π‘Žπ›Ύπ‘Ÿπ‘—(𝜏(𝑑)) ≀ 0 (πœ‰(𝑑)(πœ”β€²(𝑑))𝛾)β€² + π‘Žπ›Ύ[πœ‰(𝜏(𝑑))(πœ”β€²(𝜏(𝑑)))𝛾]β€² +𝑄(𝑑) βˆ‘ [π‘₯𝛾(𝛿𝑖(𝑑)) + 𝑝𝛾(𝛿𝑖(𝑑))π‘₯𝛾(𝛿𝑖(𝜏(𝑑))] 𝑛 𝑖=1 βˆ’ 𝐺(𝑑) βˆ‘ (1 + π‘˜ 𝑗=1 π‘Žπ›Ύ) ≀ 0 by using(9) the last inequality yields: (πœ‰(𝑑)(πœ”β€²(𝑑))𝛾)β€² + π‘Žπ›Ύ[πœ‰(𝜏(𝑑))(πœ”β€²(𝜏(𝑑)))𝛾]β€² + πœ†π‘„(𝑑) βˆ‘ πœ”π›Ύ(𝛿𝑖(𝑑)) 𝑛 𝑖=1 βˆ’ π‘˜πΊ(𝑑)(1 + π‘Žπ›Ύ) ≀ 0 Let 𝛿(𝑑) = min 𝑑β‰₯𝑑1 {𝛿𝑖(𝑑), 𝑖 = 1,2 … , 𝑛 }, by lemma 2.1, πœ”(𝑑) is positive and increasing so there exist a constant 𝑏 > 0, and 𝑑2 β‰₯ 𝑑1 such that πœ”(𝑑) β‰₯ 𝑏, 𝑑 β‰₯ 𝑑2. Hence the last inequality leads to: (πœ‰(𝑑)(πœ”β€²(𝑑))𝛾)β€² + π‘Žπ›Ύ[πœ‰(𝜏(𝑑))(πœ”β€²(𝜏(𝑑)))𝛾]β€² + π‘›πœ†π‘„(𝑑)πœ”π›Ύ(𝛿(𝑑)) βˆ’ π‘˜πΊ(𝑑)(1 + π‘Žπ›Ύ) ≀ 0, (11) π‘›πœ†π‘„(𝑑)𝑏𝛾 ≀ βˆ’(πœ‰(𝑑)(πœ”β€²(𝑑))𝛾)β€² βˆ’ π‘Žπ›Ύ[πœ‰(𝜏(𝑑))(πœ”β€²(𝜏(𝑑)))𝛾 ] β€² + π‘˜πΊ(𝑑)(1 + π‘Žπ›Ύ), (12) Consequently, by integrating (12) from 𝑑2 to 𝑑 yields π‘›πœ†π‘π›Ύ ∫ 𝑄(𝑠) 𝑑 𝑑2 𝑑𝑠 ≀ βˆ’ ∫ (πœ‰(𝑠)(πœ”β€²(𝑠))𝛾)β€² 𝑑 𝑑2 𝑑𝑠 βˆ’ π‘Žπ›Ύ ∫ [πœ‰(𝜏(𝑠))(πœ”β€²(𝑠))𝛾]β€² 𝑑 𝑑2 𝑑𝑠 + π‘˜(1 + π‘Žπ›Ύ) ∫ 𝐺(𝑠)𝑑𝑠 𝑑 𝑑2 π‘›πœ†π‘π›Ύ ∫ 𝑄(𝑠) 𝑑 𝑑2 𝑑𝑠 ≀ πœ‰(𝑑2)(πœ”β€²(𝑑2))𝛾 + π‘Žπ›Ύπœ‰(𝜏(𝑑2))[πœ”β€²(𝜏(𝑑2))]𝛾 + π‘˜(1 + π‘Žπ›Ύ) ∫ 𝐺(𝑠) 𝑑 𝑑2 𝑑𝑠, (13) Hence by (M3) it follows from (13) ∫ 𝑄(𝑠) ∞ 𝑑2 𝑑𝑠 < ∞, contradicts (𝑀4). Theorem 2.2. Assume that π‘žπ‘–(𝑑) ≀ 0, 𝑖 = 1,2, … , 𝑛, βˆ‘ π‘Ÿπ‘–(𝑑)π‘˜ 𝑗=1 β‰₯ 0, πœ‰β€²(𝑑) > 0 on [𝑑0, ∞) . Let (M1) βˆ’ (M4) hold, and for any continuous functions𝑒(𝑑), 𝑣(𝑑), 𝑒𝑣 > 0, there exists πœ† > 0, such that (8) holds, in addition to the condition lim sup π‘‘β†’βˆž ∫ [ 1 πœ‰(𝑠) ∫ βˆ‘ |π‘žπ‘–(𝑣)|[1 βˆ’ 𝑝(𝛿𝑖(𝑣))]𝛾 𝑛 𝑖=1 𝛼(𝑠) 𝑠 𝑑𝑣] 1 𝛾 𝑑𝑠 𝛼(𝑑) 𝑑 > 1. (14) Then every solution of Equation (1) oscillates. IHJPAS. 2024, 38 (1) 412 Proof. Assume that Equation (1) has eventually positive solution π‘₯(t), that is π‘₯(𝑑) > 0, π‘₯(𝜏(𝑑)) > 0, π‘₯(𝛿𝑖(𝑑)) > 0, 𝑖 = 1,2, … , 𝑛. From Equation (1) we get (πœ‰(𝑑)(πœ”β€²(𝑑))𝛾)β€² β‰₯ 0, based on Lemma 2.2, there are two cases that need to be investigated: (a) πœ”β€²(𝑑) > 0, and lim π‘‘β†’βˆž π‘₯(𝑑) = ∞. (b) πœ”β€²(𝑑) < 0, and lim π‘‘β†’βˆž πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 = 0. Case (a) proceeding as in the proof of theorem 2.1, we conclude that (11) holds. Letting 𝑧(𝑑) = πœ‰(𝑑)(πœ”β€²(𝑑))𝛾. (15) Then 𝑧(𝑑) is positive and non-decreasing, hence (11) becomes: 𝑧′(𝑑)+π‘Žπ›Ύ[𝑧(𝜏(𝑑))] β€² + π‘›πœ†π‘„(𝑑)πœ”π›Ύ(𝛿(𝑑)) < π‘˜πΊ(𝑑)(1 + π‘Žπ›Ύ), for 𝑑 β‰₯ 𝑑2 (16) πœ”(𝑑) is positive and increasing so there exist a constant 𝑏 > 0, and 𝑑3 β‰₯ 𝑑2 such that πœ”(𝑑) β‰₯ 𝑏, 𝑑 β‰₯ 𝑑3. Therefor (16) reduce to 𝑧′(𝑑)+π‘Žπ›Ύ[𝑧(𝜏(𝑑))] β€² + π‘›πœ†π‘„(𝑑)𝑏𝛾 < π‘˜πΊ(𝑑)(1 + π‘Žπ›Ύ), for 𝑑 β‰₯ 𝑑3 (17) Integration (17) from 𝑑3 to 𝑑, where 𝑑 is sufficiently large 𝑑3, leads to ∫ 𝑧′(𝑠)𝑑𝑠 𝑑 𝑑3 +π‘Žπ›Ύ ∫ (𝑧(𝜏(𝑠)))′𝑑𝑠 𝑑 𝑑3 + π‘›πœ†π‘π›Ύ ∫ 𝑄(𝑠)𝑑𝑠 𝑑 𝑑3 < π‘˜(1 + π‘Žπ›Ύ) ∫ 𝐺(𝑠)𝑑𝑠 𝑑 𝑑3 Since 𝑧(𝑑) is non-decreasing, then the last inequality becomes: 𝑧(𝑑) βˆ’ 𝑧(𝑑3) + π‘Žπ›Ύπ‘§(𝜏(𝑑)) βˆ’ π‘Žπ›Ύπ‘§(𝜏(𝑑3)) + π‘›πœ†π‘π›Ύ ∫ 𝑄(𝑠) 𝑑 𝑑3 𝑑𝑠 < π‘˜ (1 + π‘Žπ›Ύ) ∫ 𝐺(𝑠)𝑑𝑠 𝑑 𝑑3 , There fore βˆ’π‘§(𝑑3)𝛾 βˆ’ π‘Žπ›Ύπ‘§(𝑑3) + π‘›πœ†π‘π›Ύ ∫ 𝑄(𝑠) 𝑑 𝑑3 𝑑𝑠 < π‘˜ (1 + π‘Žπ›Ύ) ∫ 𝐺(𝑠)𝑑𝑠, 𝑑 𝑑3 (18) As 𝑑 β†’ ∞ a contradiction will be got in (18). Case (b) in this case πœ”(𝑑) > 0, πœ”β€²(𝑑) < 0, lim π‘‘β†’βˆž πœ‰(𝑑)(πœ”β€²(𝑑))𝛾 = 0, (πœ‰(𝑑)(πœ”β€²(𝑑)) 𝛾 )β€² β‰₯ 0 Since πœ‰β€²(𝑑) > 0 and πœ”(𝑑) is positive decreasing, so it can be conclude that πœ”β€²β€²(𝑑) β‰₯ 0, for 𝑑 β‰₯ 𝑑2, πœ”(𝑑) > π‘₯(𝑑), π‘₯(𝑑) = πœ”(𝑑) βˆ’ 𝑝(𝑑)π‘₯(𝜏(𝑑)), π‘₯(𝛿𝑖(𝑑)) = πœ”(𝛿𝑖(𝑑)) βˆ’ 𝑝(𝛿𝑖(𝑑))π‘₯(𝜏(𝛿𝑖(𝑑))) So Equation (1) become (πœ‰(𝑑)(πœ”β€²(𝑑)) 𝛾 ) β€² + βˆ‘ π‘žπ‘–(𝑑)[πœ”(𝛿𝑖(𝑑)) βˆ’ 𝑝(𝛿𝑖(𝑑))π‘₯(𝜏(𝛿𝑖(𝑑)))]𝛾 𝑛 𝑖=1 = βˆ‘ π‘Ÿπ‘—(𝑑) π‘˜ 𝑗=1 . (19) By integrating (19) from 𝑑 to 𝛼(𝑑), where 𝛼(𝑑) > 𝑑 and 𝜏 (𝛿𝑗 (𝛼(𝛼(𝑑)))) < 𝑑, 𝛿𝑗(𝑑) = min{𝛿𝑖(𝑑) , 𝑖 = 1,2, … , 𝑛}, We get βˆ’πœ‰(𝑑)(πœ”β€²(𝑑)) 𝛾 β‰₯ βˆ’ ∫ βˆ‘ π‘žπ‘–(𝑠)[πœ”(𝛿𝑖(𝑑)) βˆ’π‘› 𝑖=1 𝛼(𝑑) 𝑑 𝑝(𝛿𝑖(𝑑))πœ”(𝜏(𝛿𝑖(𝑑)))]𝛾 𝑑𝑠, IHJPAS. 2024, 38 (1) 413 βˆ’πœ‰(𝑑)(πœ”β€²(𝑑)) 𝛾 β‰₯ βˆ’ ∫ βˆ‘ π‘žπ‘–(𝑠)πœ”π›Ύ(𝜏(𝛿𝑖(𝑑)))[1 βˆ’ 𝑝(𝛿𝑖(𝑑))]𝛾 𝑛 𝑖=1 𝛼(𝑑) 𝑑 𝑑𝑠, πœ”β€²(𝑑) ≀ πœ” (𝜏 (𝛿𝑗(𝛼(𝑑)))) [ 1 πœ‰(𝑑) ∫ βˆ‘ π‘žπ‘–(𝑠)[1 βˆ’ 𝑝(𝛿𝑖(𝑑))]𝛾 𝑛 𝑖=1 𝛼(𝑑) 𝑑 𝑑𝑠] 1 𝛾 . (20) Where 𝜏(𝛿𝑗(𝑑)) = min{𝜏(𝛿𝑖(𝑑)) , 𝑖 = 1,2, … , 𝑛], integrating (20) from 𝑑 to 𝛼(𝑑) we get πœ”(𝛼(𝑑)) βˆ’ πœ”(𝑑) ≀ πœ”(𝜏(𝛿𝑗(𝛼(𝛼(𝑑))))) ∫ [ 1 πœ‰(𝑠) ∫ βˆ‘ π‘žπ‘–(𝑣)[1 βˆ’ 𝑝(𝛿𝑖(𝑣))]𝛾 𝑛 𝑖=1 𝛼(𝑠) 𝑠 𝑑𝑣] 1 𝛾 𝑑𝑠 𝛼(𝑑) 𝑑 , 1 β‰₯ πœ”(𝑑) πœ” (𝜏 (𝛿𝑗 (𝛼(𝛼(𝑑))))) β‰₯ βˆ’ ∫ [ 1 πœ‰(𝑠) ∫ βˆ‘ π‘žπ‘–(𝑣)[1 βˆ’ 𝑝(𝛿𝑖(𝑣))]𝛾 𝑛 𝑖=1 𝛼(𝑠) 𝑠 𝑑𝑣] 1 𝛾 𝑑𝑠 𝛼(𝑑) 𝑑 . The last inequality contradicts the condition (14), thus case not valid also, hence every solution of Equation (1) oscillates. The proof is complete. 3. Examples In this section, two examples are given to illustrate the fulfillment of all necessary and sufficient conditions for the results presented in the previous section. Example 3.1. Consider the following Emden-Fowler Equation: [π‘₯(𝑑) + 1 2 π‘₯(𝑑 βˆ’ πœ‹ )] β€²β€² = βˆ’ 1 2 π‘₯(𝑑 βˆ’ 2πœ‹) βˆ’ 1 4 , 𝑑 β‰₯ 0. (21) Where πœ‰(𝑑) = 1 , 𝑝(𝑑) = 1 2 , 𝜏(𝑑) = 𝑑 βˆ’ πœ‹ , π‘ž1(𝑑) = 𝑄(𝑑) = 1 2 , π‘Ÿ1(𝑑) = βˆ’ 1 2 , and 𝛿(𝑑) = 𝑑 βˆ’ 2πœ‹, 𝑖 = 1,2, … , 𝑛, 𝛾 = 1. In reality M1 βˆ’ M3 are hold for every 𝑑 β‰₯ 𝑑0 = 0. And ∫ 𝑄(𝑠)𝑑𝑠 ∞ 𝑑0 = ∫ 𝑑𝑠 ∞ 0 = ∞. Then (M4) is holds for every 𝑑 β‰₯ 1 2 . Recall that (8) hold for Ξ» = 1. Hence all the conditions of Theorem 2.1 satisfy that according to Theorem 2.1, each solution of Equation (1) oscillates for example π‘₯(𝑑) = sin 𝑑 βˆ’ 1 2 such as this oscillation solution. See Figure (1) Figure 1. π‘₯(𝑑) = sin 𝑑 βˆ’ 1 2 IHJPAS. 2024, 38 (1) 414 Example 3. 2. Consider the following Emden-Fowler Equation: [π‘₯(𝑑) + 2π‘’βˆ’πœ‹π‘₯(𝑑 βˆ’ πœ‹ )]β€²β€² βˆ’ 4π‘’βˆ’ 3πœ‹ 2 π‘₯ (𝑑 βˆ’ 3πœ‹ 2 ) = π‘’βˆ’π‘‘, 𝑑 β‰₯ 0. (22) Where πœ‰(𝑑) = 1 , 𝑝(𝑑) = 2π‘’βˆ’πœ‹, 𝜏(𝑑) = 𝑑 βˆ’ πœ‹ , π‘Ÿ(𝑑) = π‘’βˆ’π‘‘, π‘ž1(𝑑) = 𝑄(𝑑) = βˆ’4π‘’βˆ’ 3πœ‹ 2 , and 𝛿(𝑑) = 𝑑 βˆ’ 3πœ‹ 2 , 𝛾 = 1. In reality M1 βˆ’ M3 are hold for every 𝑑 β‰₯ 𝑑0 = 0. And ∫ |𝑄(𝑠)|𝑑𝑠 ∞ 𝑑0 = ∫ 4π‘’βˆ’ 3πœ‹ 2 𝑑𝑠 ∞ 0 = ∞. Then (M4) is holds for every 𝑑 β‰₯ 0. Recall that (8) hold for Ξ» = 1. Hence all the conditions of Theorem 2.2 satisfy that according to Theorem 2.2, each solution of Equation (1) oscillates, for example π‘₯(𝑑) = π‘’βˆ’π‘‘(sin 𝑑 βˆ’ 1) is such an oscillatory solution. See Figure (2) Figure 𝟐. πœ‘(𝑑) = π‘’βˆ’π‘‘(sin 𝑑 βˆ’ 1) 4. Conclusion In this paper, we have studied the oscillation property of the solutions of neutral second-order differential Equations of the Emden-Fowler type. Some of the extracted conditions are the development of conditions known in the references, which ensure that either each solution of this Equation oscillates, or each nonoscillatory solution convergence to zero or tends to infinity as 𝑑 β†’ ∞. Some examples are presented to clarify the results obtained. Acknowledgments I would like to extend my thinks to my teacher and supervisor, Dr. Hussain Ali Mohamad, for providing guidance and feedback throughout this project. Conflict of Interest The authors declare that there are no competing interests regarding the publication of this paper. Funding This work is not supported by any the Foundation. Ethical Clearance Ethics of scientific research were carried out in accordance with international conditions. IHJPAS. 2024, 38 (1) 415 References 1. Ahmed F A , Mohamad H A. Oscillation and Asymptotic Behavior of Second Order Half Linear Neutral Dynamic Equations . Iraqi Journal of Science. 2022; 63(12):5413-542. https://doi.org/10.24996/ijs.2022.63.12.27 2. Wu Y, Yu Y, Zhang J, Xiao J. Oscillation criteria for second order Emden-Fowler functional differential Equations of neutral type. Journal of I and Applications. 2016: 1-11. https://doi.org/10.1186/s13660-016- 1268-9 3. Baculikova B, DΕΎurina J. Oscillation theorems for second order neutral differential Equations. Computers & Mathematics with Applications.2011;61(1): 94-99. https://doi.org/10.1016/j.camwa.2010.10.035 4. Mehta B N, Aris R. A note on a form of the Emden-Fowler Equation. Journal of Mathematical Analysis and Applications. 1971; 6(3): 611-621. https://doi.org/10.1016/0022-247X(71)90043-6 5. Mohamad H A, Ketab S N. Oscillation Solution for Nonlinear Third Order Neutral Differential Equations. Iraqi J. of Science. Special Issue, Part B, 2016: 412-417. https://search.emarefa.net/detail/BIM- 762600 6. Mohamad H A, Shehab L M. Oscillations of Third Order Half Linear Neutral Differential Equations. Baghdad Science Journal. 2015; 12(3): 625-631. https://doi.org/10.21123/bsj.2015.12.3.625-631 7. Moaaz O, Ramos H, Awrejcewicz J. Second-order Emden-Fowler neutral differential Equations A new precise criterion for oscillation .Applied Mathematics Letters. 2021; 118: 107-172. https://doi.org/10.1016/j.aml.2021.107172 8. Xu R, Xia Y. A note on the oscillation of second-order nonlinear neutral functional differential Equations. International Journal of Contemporary Mathematical Sciences.2008 ; 3(29-32): 1441-1450. https://doi.org/10.1016/j.camwa.2008.09.004 9. Thandapani E , Tongxing Li. On the Oscillation of Third-Order Quasi-Linear Neutral Functional Differential Equations. Archivum Mathematicum (BRNO) Tomus .2011; 47: 181–199. http://eudml.org/doc/246122 10. Hassan T S, Kong Q, El-Matary B M. Oscillation criteria for advanced half-linear Differential.Equations..of.second order .Mathematics. 2023; 11(6):1385. 11. https://doi.org/10.3390/math11061385 12. Tripathy A K, Santra S S. Necessary and sufficient conditions for oscillations to a second-order neutral differential Equations with impulses. Kragujev. J. Math. 2023; 47(12): 81–93. https://doi.org/10.46793/KgJMat2301.0 13. Mehta B N, Aris R. A Note on a Form of the Emden-Fowler Equation . Journal of Mathematical Analysis and Applications. 1971; 63: 611-621. https://doi.org/10.1016/0022-247X(71)90043-6 14. Vidhyaa K S, Graef J R, Thandapani E. New oscillation results for third-order half-linear neutral differential Equations. Mathematics. 2020; 8(3): 325. https://doi.org/10.3390/math8030325 15. Dassios I, Liu M, Milano F. On the stability analysis of systems of neutral delay differential Equations.Circuits, Systems, and Signal Processing. 2019; 38: 1639-1653. https://doi.org/10.1007/s00034- 018-0943-0 16. Li T, Rogovchenko YV. Oscillation criteria for second-order superlinear Emden–Fowler neutral differential Equations. Monatshefte FΓΌr Math.2017; 184: 489–500.https://doi.org/10.1007/s00605-017- 1039-9 17. Marappan S K, Almutairi A, Iambor LF, Bazighifan O. Oscillation of Emden–Fowler-Type Differential Equations with Non-Canonical Operators and Mixed Neutral Terms. Symmetry. 2023; 15(2): 553. https://doi.org/10.3390/sym15020553 https://doi.org/10.1016/j.camwa.2010.10.035 https://doi.org/10.1016/0022-247X(71)90043-6 https://doi.org/10.21123/bsj.2015.12.3.625-631 https://doi.org/10.1016/j.aml.2021.107172 https://doi.org/10.1016/j.camwa.2008.09.004 https://doi.org/10.3390/math11061385 file:///C:/doi.org/10.46793/KgJMat2301.0 https://doi.org/10.1016/0022-247X(71)90043-6 https://doi.org/10.3390/math8030325 https://doi.org/10.1007/s00034-018-0943-0 https://doi.org/10.1007/s00034-018-0943-0 https://doi.org/10.3390/sym15020553 IHJPAS. 2024, 38 (1) 416 18. Baty H. Solving higher-order Lane- Emden Fowler type Equations using physics-informed neural networks: benchmark tests comparing soft and hard constraints. 2023; arXiv preprint arXiv:2307. 07302. https://doi.org/10.48550/arXiv.2307.07302 19. Naeif A A, Mohamad H A. Almost Oscillatory Solutions of a Three-Dimensional Half-Linear Neutral differential System of the Second Order. Journal of University of Babylon for Pure and Applied Sciences. 2023; 31)1); 48–71. https://doi.org/10.29196/jubpas.v31i1.4528 20. Suha N, Singh R. An efficient new numerical algorithm for solving Emden–FowleSr pantograph differential Equation using Laguerre polynomials. Journal of Computational Science. 2023; 72:102108 .https://doi.org/10.1016/j.jocs.2023.102108 21. Rufai M A, Ramos H. Numerical integration of third-order singular boundary-value problems of Emden– Fowler type using hybrid block techniques. Communications in Nonlinear Science and Numerical Simulation. 2022; 105:106069 .https://doi.org/10.1016/j.cnsns.2021.106069 22. Mahdy A M S. A numerical method for solving the nonlinear Equations of Emden-Fowler models. Journal of Ocean Engineering and Science. 2022. https://doi.org/10.1016/j.joes.2022.04.019 23. Singh R, Singh M. An optimal decomposition method for analytical and numerical solution of third-order Emden–Fowler type Equations .Journal of Computational Science. 2022; 63: 101790. https://doi.org/10.1016/j.jocs.2022.101790 https://doi./ https://doi.org/10.29196/jubpas.v31i1.4528 https://doi.org/10.1016/j.jocs.2023.102108 https://doi.org/10.1016/j.cnsns.2021.106069 https://doi.org/10.1016/j.joes.2022.04.019 https://doi.org/10.1016/j.jocs.2022.101790