Abstract
A computational method is developed for a family of functional differential equations in one independent variable with a single deviating argument, assumed to be a delay. These equations are mapped into a finite system of ODEs, a subsystem of which involves only the function controlling the delay. Key features include efficiency during method of steps, freedom from Jacobians and root finding techniques, and computing a continuous approximation. The nonlinear differential equations may be retarded, neutral or advanced. The method is established for state dependent delays, stiff equations, discontinuous initial history, a specific loss of monotonicity in the delay and extended naturally to distributed delays. When restricted to ODEs, the method extends naturally to PDEs; the hope is the delayed version here can eventually be extended to delayed PDEs. Conditions for convergence of the approximation are established, and results of numerical experiments are reported to indicate robustness of the implementation.
Similar content being viewed by others
References
Hale, J.K., Lunel, S.M.V.: Introduction to Functional Differential Equations. Springer-Verlag, New York (1993)
Carothers, D.C., Parker, G.E., Sochacki, J.S., Warne, P.G.: Some properties of solutions to polynomial systems of equations. Electron. J. Differ. Equ. 2005(40), 1–17 (2005)
Parker, G.E., Sochacki, J.S.: Implementing the Picard iteration. Neural Parallel Sci. Comput. 4(1), 97–112 (1996)
Sochacki, J.S.: Polynomial ODEs—examples, solutions and properties. Neural Parallel Sci. Comput. 18, 441–449 (2010)
Hu, Q.: A model of regulatory dynamics with threshold-type state-dependent delay. Math. Biosci. Eng. 15(4), 863–882 (2018)
Isaia, V.M.: Nonlinear differential equations with deviating arguments and their approximation via a Parker–Sochacki approach. Electron. J. Differ. Equ. 2017(68), 1–23 (2017)
Baker, C.T.H., Paul, C.A.H., Willé, D.R.: Issues in the numerical solution of evolutionary delay differential equations. Adv. Comput. Math. 3, 171–196 (1995)
Fehlberg, E.: Numerical Integration of Differential Equations by Power Series Expansions, Illustrated by Physical Examples. Technical Report NASA-TN-D-2356, NASA (1964)
Parker, G.E., Sochacki, J.S.: A Picard–Maclaurin theorem for initial value PDEs. Abstr. Appl. Anal. 5(1), 47–63 (2000)
Norkin, S.B.: Differential Equations of the Second Order with Retarded Argument. Translations of Mathematical Monographs (L. J. Grimm) vol. 31, AMS, Providence (1972)
Bellen, A., Guglielmi, N.: Solving neutral delay differential equations with state-dependent delays. J. Comput. Appl. Math. 229, 350–362 (2009)
Bellen, A., Zennaro, M.: Numerical Methods for Delay Differential Equations. Oxford University Press, Oxford (2003)
Castleton, R.N., Grimm, L.J.: A first order method for differential equations of neutral type. Math. Comput. 27(123), 571–577 (1973)
Bellman, R.E., Buell, J.D., Kalaba, R.E.: Numerical integration of a differential-difference equation with a decreasing time-lag. Commun. ACM 8(4), 227–228 (1965)
Warne, P.G., Warne, D.A.P., Sochacki, J.S., Parker, G.E., Carothers, D.C.: Explicit a-priori error bounds and adaptive error control for approximation of nonlinear initial value differential systems. Comput. Math. Appl. 52, 1695–1710 (2006)
Yeung, M.K.S., Strogatz, S.H.: Time delay in the Kuramoto model of coupled oscillators. Phys. Rev. Lett. 82(3), 648–651 (1996)
Bellen, A., Guglielmi, N., Ruehli, A.E.: Methods for linear systems of circuit delay differential equations of neutral type. IEEE Trans. Circuits Syst. I Fundam. Theory Appl. 46(1), 212–215 (1999)
de Pillis, L., Gallegos, A., Radunskaya, A.: A model of dendritic cell therapy for melanoma. Front. Oncol. 3(56), 1–14 (2013)
Guglielmi, N., Hairer, E.: Computing breaking points in implicit delay differential equations. Adv. Comput. Math. 29(3), 229–247 (2008)
Bellman, R.E., Buell, J.D., Kalaba, R.E.: Mathematical experimentation in time-lag modulation. Commun. ACM 9(10), 752–754 (1966)
Acknowledgements
The author would like to thank Dr. James Sochacki and Dr. Anthony Tongen for their contributions which enabled this manuscript to be produced. The author would also like to thank Dr. Lisette de Pillis, who made the author aware of the blood, spleen and tumor problem. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. The author would also like to thank the (hard working) referees for their helpful suggestions, which ultimately made this a better paper.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The author declares that they have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix
Appendix
For a full discussion of the BST model, see [19] and references therein. This is the list of 9 components from the full system of DEs in the BST problem, see page 35, which are all functions of t
This is the list of parameters (all constant) from the same system:
The functions \(v_{blood}(t)\) and \(v_{tumor}(t)\) are nonhomogeneous (forcing) terms. This is a list of known functions in the same system:
The full system of equations, using the above and prime for derivative, is
Auxiliary variables were used for \(\mu _{SB}\), \(\mu _{BTE}\), \(E^a_T T^{-1}\), \({\mathscr {D}}\), \(1 - e^{\frac{-\mu _{BS} D_B}{\max D}}\), \((\theta _D + D_S(t-\tau _D))^{-1}\), \((q + T)^{-1}\), and \((\cdot )^{l}\) since \(l = 2/3\), to polynomial ones.
Rights and permissions
About this article
Cite this article
Isaia, V.M. A Fixed Point Approach to Simulation of Functional Differential Equations with a Delayed Argument. J Dyn Diff Equat 35, 1047��1082 (2023). https://doi.org/10.1007/s10884-021-10081-7
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10884-021-10081-7