Computation of periodic solutions to models of infectious disease dynamics and immune response. (English) Zbl 1466.92192

Summary: The paper is focused on computation of stable periodic solutions to systems of delay differential equations modelling the dynamics of infectious diseases and immune response. The method proposed here is described by an example of the well-known model of dynamics of experimental infection caused by lymphocytic choriomeningitis viruses. It includes the relaxation method for forming an approximate periodic solution, a method for estimating the approximate period of this solution based on the Fourier series expansion, and a Newton-type method for refining the approximate period and periodic solution. The results of numerical experiments are presented and discussed. The proposed method is compared to known ones.


92D30 Epidemiology
34K13 Periodic solutions to functional-differential equations
34L16 Numerical approximation of eigenvalues and of other parts of the spectrum of ordinary differential operators
Full Text: DOI


[1] G. A. Bocharov, Modelling the dynamics of LCMV infection in mice: conventional and exhaustive CTL responses. J. Theor. Biol. 192 (1998), No. 3, 283-308.
[2] K. Engelborghs, T. Luzyanina, K. J. in ’t Hout, D. Roose, Collocation methods for the computation of periodic solutions of delay differential equations. SIAM J. Sci. Comp. 22 (2001), No. 5, 1593-1609. · Zbl 0981.65082
[3] K. Engelborghs, T. Luzyanina, and D. Roose, Numerical bifurcation analysis of delay differential equations. J. Comp. Appl. Math. 125 (2000), No. 1-2, 265-275. · Zbl 0965.65099
[4] K. Engelborghs, T. Luzyanina, and D. Roose, Numerical bifurcation analysis of delay differential equations using DDE-BIFTOOL. ACM Trans. Math. Softw. 28 (2002), No. 1, 1-21. · Zbl 1070.65556
[5] G. E. Forsythe, M. A. Malcolm, and C. B. Moler, Computer Methods for Mathematical Computations. Prentice-Hall, Englewood Cliffs, 1977. · Zbl 0361.65002
[6] E. Hairer and G. Wanner, Solving Ordinary Differential Equations. Springer-Verlag, Berlin, 1996. · Zbl 0859.65067
[7] J. K. Hale, Theory of Functional Differential Equations. Springer-Verlag, New York, 1977.
[8] D. Kahaner, C. Moler, and S. Nash, Numerical Methods and Software. Prentice-Hall, Englewood Cliffs, 1977.
[9] G. I. Marchuk, Mathematical Models in Immunology. Optimization Software Inc. Publications Division, New York, 1983. · Zbl 0685.92002
[10] G. I. Marchuk, R. V. Petrov, A. A. Romanyukha, and G. A. Bocharov, Mathematical model of antiviral immune response. I. Data analysis, generalized picture construction and parameters evaluation for hepatitis B. J. Theor. Biol. 151 (1991), No. 1, 1-40.
[11] Yu. M. Nechepurenko, M. Yu. Khristichenko, D. S. Grebennikov, and G. A. Bocharov, Bistability analysis of virus infection models with time delays. Discrete and Continuous Dynamical Systems Series S 13 (2020), No. 9, 2385-2401. · Zbl 1451.92303
[12] D. Roose and R. Szalai, Continuation and Bifurcation Analysis of Delay Differential Equations/ Numerical Continuation Methods for Dynamical Systems: Path Following and Boundary Value Problems. Springer Netherlands, Dordrecht, 2007.
[13] R. Seydel, Practical Bifurcation and Stability Analysis. Springer-Verlag, New York, 2009.
[14] J. Sieber, K. Engelborghs, T. Luzyanina, G. Samaey, and D. Roose, DDE-BIFTOOL Manual — Bifurcation Analysis of Delay Differential Equations. ArXiv:1406.7144 [Math], June 27, 2014.
[15] E. V. Sklyarova, Yu. M. Nechepurenko, and G. A. Bocharov, Numerical steady state analysis of the Marchuk-Petrov model of antiviral immune response. Russ. J. Numer. Anal. Math. Modelling 35 (2020), No. 2, 95-110. · Zbl 07207776
[16] G. W. Stewart, Simultaneous iteration for computing invariant subspaces of non-Hermitian matrices. Numer. Math. 25 (1976), No. 2, 123-126. · Zbl 0328.65025
[17] G. Szego, Orthogonal Polynomials, American Mathematical Society Colloquium Publications. AMS, New York, 1939. · JFM 65.0286.02
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.