Abstract
Obtaining a solution of a given SDE is essential in neuroscience, especially, in modeling transmission of nerve impulses between neurons through myelin substance. This paper analyzes a particular scalar differential equation (SDE). The current scalar model involves two categories of differential equations–advanced and delayed–based on the domain of the independent variable. The results are consistent with existing literature as the advance/delay parameter approaches unity. Theoretical and graphical analyses of the solution’s properties are presented. To the best of our knowledge, this is the first study to analyze this form of SDE.
1 Introduction
In ordinary differential equations (ODEs), the equation is typically classified as a delay differential equation (DDE) in the domain , since for any , with serving as the delay parameter. Conversely, the equation is considered an advanced differential equation (ADE) in the domain , as for all , with interpreted as the advance parameter. However, if an ODE involves both delay and advance terms in distinct but connected domains, it is more appropriately classified as a scalar differential equation (SDE). In the examples above, the terms and involve positive coefficients of the independent variable , allowing straightforward classification of the respective equations as DDE and ADE.
A question arises here: what is the type of the second ODE if is changed to ? Answering this question requires two steps to determine the domains of for which (delay) and (advance). The first step, implies , while the second step, leads to . Based on this, the ODE can be classified as an ADE in the domain , and as a DDE in the domain . Hence, we may refer to the ODE as an SDE because it involves both types of advance and delay equations, as pointed out in Refs. [, ]. Another important observation concerns the central point connecting the two domains, which is . This central point plays a fundamental role in deriving the analytical solutions of a given SDE, as will be demonstrated later. It is also useful to distinguish between proportional delay parameters and pure delay parameters. In the DDE , the parameter is referred to as a pure delay parameter. However, other types of DDEs involve proportional delay parameters, such as in the pantograph equation (PE) , [, ]. The PE has applications in modeling the behavior of overhead catenary systems for railway electrification [–], the dynamic response of trolley wire overhead contact systems for electric railways [], and current collection systems in electric locomotives []. Several authors have analyzed the PE in detail [–]. Another notable example is the Ambartsumian equation (AE), given by , where . This equation has practical significance in astronomy, particularly in studies of surface brightness in the Milky Way [–]. In these models, and are considered proportional delay parameters. For the pantograph model, implies indicating a delay for all . Similarly, the Ambartsumian model also represents a delay.
In this paper, we consider the following general form of the SDE:where , , and are real constants. It can be readily shown that Equation 1 represents an advanced equation in the domain , while it becomes a delayed equation for . Finding a solution to Equation 1 poses a significant challenge and, to the best of our knowledge, may be considered for the first time. Moreover, standard methods such as the Adomian decomposition method (ADM) [–], the homotopy perturbation method [, ], and the Laplace transform (LT) [–] may encounter difficulties when applied to such problems.
To address this, a direct series approach is developed to solve the advanced equation. A closed-form expression of the series is obtained, and its convergence is established theoretically. These results are then used to construct the solution for the delayed equation. Several existing results in the literature can be recovered as special cases of the present findings. In addition, the properties of the obtained solutions are analyzed both theoretically and graphically. Finding a solution for a SDE is helpful for understanding the transmission of nerve impulses between neurons through myelin substance which covers all the nerves in the brain and nervous system in humans []. Other areas of applications can be further extended to involve some recent dynamical systems [–] and relatively new physical phenomena [, ].
2 Advanced equation
In the domain , SDE (1) becomes an advanced equation since , . Moreover, in the advanced equation domain, , see Refs. [, ] for details. Accordingly, the value of the function is unknown, which prevents the application of the step method to solve the advanced equation: Before discussing the main objective of this section, it is important to note that the condition:must be satisfied by any solution to Equation 2 in addition to the initial condition (IC) .
2.1 Closed-form series solution
An effective solution for a given model can be derived as a closed-form series solution. The solution in such a form facilitates numerical calculations and also leads to easier analysis to study the properties/behavior of the physical system. Let us attempt a series of solutions in the form of: This assumption yields:andSubstituting Equations 4–6 into Equation 2 leads to:Hence,From (4), we can write:Employing (8), we obtain:orwhereBy applying IC to Equation 11, we get:Substituting (13) into (11) yields:Equation 14 declares that IC is satisfied automatically. Let us now check the satisfaction of condition (3). For this purpose, from the solution (14), we obtain:andwhere Equation 12 is implemented to calculate . The last two equations show that condition (3) is also satisfied. The next step is to examine the convergence of the obtained series solution, which is discussed in the following subsection.
2.2 Convergence analysis
To provide a theoretical proof of the convergence of the series solution (14), it is sufficient to prove the convergence of the series in the domain .
For, the series:converge in the domain.
Proof.Let us define,Applying the ratio test:Hence,The limit on the right-hand side of the last equation tends toward zero as for every and . At , the value according to . In this case, Equation 20 becomes:which also tends toward zero, thereby completing the proof.
Remark 1Through a similar analysis, we can easily prove that the seriesis convergent for all.
2.3 Special case and exact solution
In this section, we show that the obtained series solution in Section 2.1 converges to the exact hyperbolic and trigonometric forms when under the conditions and , respectively. We consider in Equation 2 and then extract the solution of the corresponding advanced equation: In this case, the solution given by Equation 14 reads:where in Equation 12 becomes:This equation can be used to generate the following equations for the even-order coefficients and odd-order coefficients as follows:The numerator of solution (23) can be written as follows:Similarly, the denominator of solution (23) can be written as follows:Substituting (26) and (27) into (23), we obtain the exact hyperbolic solution:Moreover, if we rewrite the coefficients and as follows:then, we can arrive at the exact periodic solution:Solution (30) agrees with the corresponding values obtained in Ref. [] for the advanced Equation 22.
3 Delay equation
It may be useful to divide the domain into two intervals, and . This is simply because the value of in each of the above two intervals can be assigned a certain value, as described in the next subsections. To achieve our target, we first denote as the solution in the interval ; hence,where is given by Equation 11.
3.1 Solution in interval
In the interval we find that and accordingly Equation 31 gives:The advanced equation in this interval takes the form:subject toSubstituting (32) into (33) results in the following ODE:Solving this ODE under Condition (34) yieldswhereThis integral appears complex; however, it can be evaluated analytically in terms of the generalized incomplete gamma function defined by:The integral (37) can be determined by:Therefore, the solution (36) takes the following form:The series on the right-hand side of this equation must also be checked for convergence, which is discussed in the next theorem.
For, the seriesconverges in the domain.
Proof.Assume that,Applying the ratio test:i.e.,We obtain:Therefore,The limit on the right-hand side tends to zero as for every , thus completing the proof.
3.2 Solution in interval
Let us define as the solution in the previous interval . At , we get , whereIn the interval , we have which yields . Therefore, the delay equation is reduced to:The solution to this ODE is as follows:
4 Results
The objective of this section is to extract the numerical results for the convergence of the obtained series solutions for the advanced equation in the interval and for the delay equation in the intervals and . Since the obtained solutions are expressed in terms of an infinite series, which was proven theoretically for convergence, one may replace infinity with a finite number. Let us denote , , and as the -term approximate solutions for the obtained solutions in the intervals , , and , respectively. Accordingly, we obtain:andwhile can be written as follows:Figures 1–3 show the curves of the approximations , , and at when , , , , and . It can be seen in these figures that the convergence of the solutions in the above three intervals is achieved using few terms. The same conclusion applies to the curves shown in Figures 4–6 when , , , , and .
FIGURE 1
FIGURE 2
FIGURE 3
FIGURE 4
FIGURE 5
FIGURE 6
The behavior of the solution in the full domain is depicted in Figures 7, 8 for the same set of values of the constants used to generate Figures 1, 4, respectively. It should be noted that the solutions plotted in Figures 7, 8 are produced using the terms in series (50)–(52).
FIGURE 7
FIGURE 8
The two black dots shown in Figures 7, 8 represent the three intervals , and . In addition, these dots represent the approximate values of and . However, Figures 7, 8 indicate that the solution is continuous at the joint points, where and .
Regarding the continuity of the derivative , we can prove that is continuous at but discontinuous at .
This conclusion can be explained theoretically as follows. At , we obtain Equation 2. The left derivative is , and the right derivative is derived from Equation 33 as . Since , then ; hence, is always continuous at .
At , Equation 33 gives the left derivative as , i.e., . Equation 48 yields the right derivative at as .
Since , then , which leads to , where and are assumed.
5 Conclusion
A new type of differential equation was addressed and solved in this study. The model took the form of SDE, , where and . The SDE splits into an advanced equation and delay equation in the domains and , respectively. The solution of the advanced equation was obtained in a closed series form, for which convergence was theoretically proven. As tended toward unity, the series solution for the advanced equation transformed into exact hyperbolic and trigonometric forms for and , respectively. The solution of the delay equation was explicitly determined in terms of the incomplete gamma function using a stepwise method. The results agreed with those in the literature when tended toward unity. The properties of the solutions were analyzed both theoretically and graphically. The results showed that the solution was continuous over the full domain of the problem. Additionally, the derivative remained continuous at the point . It was also indicated that is discontinuous at provided that or did not vanish. The proposed approach is promising and can be further extended to include additional SDEs of more complex types. Thus, it maybe interested to extend this work to the domain of distributed parameter systems as in Refs. [–].
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
LS: Conceptualization, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Validation, Writing – review and editing. EE-Z: Conceptualization, Formal Analysis, Investigation, Methodology, Validation, Writing – original draft. AE: Formal Analysis, Investigation, Methodology, Validation, Writing – review and editing. AA: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Validation, Writing – original draft.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. The authors extend their appreciation to Prince Sattam bin Abdulaziz University for funding this research through project number (PSAU/2024/01/31342).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
AlshomraniNAMEbaidAAldosariFAljoufiMD. On the exact solution of a scalar differential equation via a simple analytical approach. Axioms (2024) 13:129. 10.3390/axioms13020129
2.
AlbidahABAlsulamiIMEl-ZaharEREbaidA, Advances in mathematical analysis for solving inhomogeneous scalar differential equation[J], AIMS Mathematics, (2024), 9(9): 23331–43. 10.3934/math.20241134
3.
El-ZaharEREbaidA. Analytical and numerical simulations of a delay model: the pantograph delay equation. Axioms (2022) 11:741. 10.3390/axioms11120741
4.
AlenazyAHSEbaidAAlgehyneEAAl-JeaidHK. Advanced study on the delay differential equation y′(t) = ay(t) + by(ct). Mathematics (2022) 10:4302. 10.3390/math10224302
5.
AndrewsHI. Third paper: calculating the behaviour of an overhead catenary system for railway electrification. Proc Inst Mech Eng (1964) 179:809–46. 10.1243/pime_proc_1964_179_050_02
6.
GilbertGDavtcsHEH. Pantograph motion on a nearly uniform railway overhead line. Proc Inst Electr Eng (1966) 113:485–92.
7.
CainePMScottPR. Single-wire railway overhead system. Proc Inst Electr Eng (1969) 116:1217–21. 10.1049/piee.1969.0226
8.
AbbottMR. Numerical method for calculating the dynamic behaviour of a trolley wire overhead contact system for electric railways. Comput J (1970) 13:363–8. 10.1093/comjnl/13.4.363
9.
OckendonJTaylerAB. The dynamics of a current collection system for an electric locomotive, proc. R. Soc. Lond. A math. Phys Eng Sci (1971) 322:447–68.
10.
FoxLMayersDOckendonJRTaylerAB. On a functional differential equation. IMA J Appl Math (1971) 8:271–307. 10.1093/imamat/8.3.271
11.
KatoTMcLeodJB. The functional-differential equation y′(x) = ay(λx) + by(x). Bull Am Math Soc (1971) 77:891–935.
12.
IserlesA. On the generalized pantograph functional-differential equation. Eur J Appl Math (1993) 4:1–38. 10.1017/s0956792500000966
13.
AmbartsumianVA. On the fluctuation of the brightness of the milky way. Doklady Akad Nauk USSR (1994) 44:223–6.
14.
PatadeJBhalekarS. On analytical solution of Ambartsumian equation. Natl Acad Sci Lett (2017) 40:291–3. 10.1007/s40009-017-0565-2
15.
AlyoubiRMSEbaidAEl-ZaharERAljoufiMD. A novel analytical treatment for the Ambartsumian delay differential equation with a variable coefficient[J]. AIMS Mathematics, (2024), 9(12): 35743–58. 10.3934/math.20241696
16.
EbaidACattaniCAl JuhaniASEl-ZaharER. A novel exact solution for the fractional Ambartsumian equation. Adv Differ Equ (2021) 88:88. 10.1186/s13662-021-03235-w
17.
AdomianG. Solving frontier problems of physics: the decomposition method. Boston, MA, USA: Kluwer Acad (1994).
18.
WazwazAM. Adomian decomposition method for a reliable treatment of the Bratu-type equations. Appl Math Comput (2005) 166:652–63. 10.1016/j.amc.2004.06.059
19.
DuanJSRachR. A new modification of the Adomian decomposition method for solving boundary value problems for higher order nonlinear differential equations. Appl Math Comput (2011) 218:4090–118. 10.1016/j.amc.2011.09.037
20.
LiWPangY. Application of Adomian decomposition method to nonlinear systems. Adv Differ Equ (2020) 2020:67. 10.1186/s13662-020-2529-y
21.
AyatiZBiazarJ. On the convergence of homotopy perturbation method. J Egypt Math Soc (2015) 23:424–8. 10.1016/j.joems.2014.06.015
22.
AlbidahABKanaanNEEbaidAAl-JeaidHK. Exact and numerical analysis of the pantograph delay differential equation via the homotopy perturbation method. Mathematics (2023) 11:944. 10.3390/math11040944
23.
KhaledSM. The exact effects of radiation and joule heating on magnetohydrodynamic marangoni convection over a flat surface. Therm Sci (2018) 22:63–72. 10.2298/tsci151005050k
24.
KhaledSM. Exact solution of the one-dimensional neutron diffusion kinetic equation with one delayed precursor concentration in Cartesian geometry. AIMS Mathematics, (2022), 7(7): 12364–73. 10.3934/math.2022686
25.
AlrebdiRAl-JeaidHK. Accurate solution for the pantograph delay differential equation via laplace transform. Mathematics (2023) 11:2031. 10.3390/math11092031
26.
AlrebdiRAl-JeaidHK. Two different analytical approaches for solving the pantograph delay equation with variable coefficient of exponential order. Axioms (2024) 13:229. 10.3390/axioms13040229
27.
ChiHBellJHassardB. Numerical solution of a nonlinear advance-delay differential equation from nerve conduction theory. J Math Biol (1986) 5:583–601. 10.1007/BF00275686
28.
SunWJinYLuG. Genuine multipartite entanglement from a thermodynamic perspective. Phys Rev A (2024) 109(4):042422. 10.1103/PhysRevA.109.042422
29.
YangYLiH. Neural ordinary differential equations for robust parameter estimation in dynamic systems with physical priors. Appl Soft Comput (2024) 169:112649. 10.1016/j.asoc.2024.112649
30.
MaPYaoTLiuWPanZChenYYangH. A 7-kW narrow-linewidth fiber amplifier assisted by optimizing the refractive index of the large-mode-area active fiber. High Pow Laser Sci Eng (2024). 12, e67. 10.1017/hpl.2024.41
31.
ShiQWangSYangJ. Small parameters analysis applied to two potential issues for maxwell-schrödinger system. J Differential Equations (2025) 424:637–59. 10.1016/j.jde.2025.01.049
32.
JiangNFengQYangXHeJ-RLiB-Z. The octonion linear canonical transform: properties and applications. Chaos, Solitons and Fractals (2025) 192:116039. 10.1016/j.chaos.2025.116039
33.
RekkabSBenhadidSAl-SaphoryR. An asymptotic analysis of the gradient remediability problem for disturbed distributed linear systems. Baghdad Sci J (2022) 19(6):1623–35. 10.21123/bsj.2022.6611
34.
Al-SaphoryRKhalidZEl JaiA. Regional boundary gradient closed loop control system and Γ*AGFO-Observer. J Phys Conf Ser (2020) 1664(1):012061. 10.1088/1742-6596/1664/1/012061
35.
Al-SaphoryRAl-ShayaARekkabS. Regional boundary asymptotic gradient reduced order observer. J Phys Conf Ser (2020) 1664:012101–19. 10.1088/1742-6596/1664/1/012101
Summary
Keywords
scalar differential equation, ordinary differential equation, delayed differential equation, advanced differential equation, series. MSC, 34K06, 34K07, 65L03
Citation
Seddek LF, El-Zahar ER, Ebaid A and Al Qarni AA (2025) Treatment of a generalized scalar differential equation: analysis and explicit solution. Front. Phys. 13:1611846. doi: 10.3389/fphy.2025.1611846
Received
14 April 2025
Accepted
04 August 2025
Published
01 September 2025
Volume
13 - 2025
Edited by
Khursheed Alam, Sharda University, India
Reviewed by
Njitacke Tabekoueng Zeric, University of Buea, Cameroon
Raheam Al-Saphory, Mustansiriyah University, Iraq
Kameshwar Sahani, Kathmandu University School of Engineering Dhulikhel Nepal, Nepal
Updates
Copyright
© 2025 Seddek, El-Zahar, Ebaid and Al Qarni.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Abdelhalim Ebaid, aebaid@ut.edu.sa
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.