Numerical solution by LMMs of stiff delay differential systems modelling an immune response

Summary. We consider the application of linear multistep methods (LMMs) for the numerical solution of initial value problem for stiff delay differential equations (DDEs) with several constant delays, which are used in mathematical modelling of immune response. For the approximation of delayed variables the Nordsieck's interpolation technique, providing an interpolation procedure consistent with the underlying linear multistep formula, is used. An analysis of the convergence for a variable-stepsize and structure of the asymptotic expansion of global error for a fixed-stepsize is presented. Some absolute stability characteristics of the method are examined. Implementation details of the code DIFSUB-DDE, being a modification of the Gear's DIFSUB, are given. Finally, an efficiency of the code developed for solution of stiff DDEs over a wide range of tolerances is illustrated on biomedical application model.

[1]  Christopher T. H. Baker,et al.  DELSOL: a numerical code for the solution of systems of delay-differential equations , 1992 .

[2]  Guriĭ Ivanovich Marchuk.  Mathematical models in immunology , 1983 .

[3]  Christopher T. H. Baker,et al.  Issues in the numerical solution of evolutionary delay differential equations , 1995, Adv. Comput. Math..

[4]  R. Mohler,et al.  Distribution models of recirculating lymphocytes , 1989, IEEE Transactions on Biomedical Engineering.

[5]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[6]  G. Marchuk,et al.  Mathematical model of antiviral immune response. II. Parameters identification for acute viral hepatitis B. , 1991, Journal of theoretical biology.

[7]  Russian Federation..  Numerical treatment of the parameter identification problem for delay-differential systems arising in immune response modelling , 1994 .

[8]  G. Adomian,et al.  Solution of the Marchuk model of infectious disease and immune response , 1986 .

[9]  The Stability of Variable-Stepsize Nordsieck Methods , 1983 .

[10]  Rolf Dieter Grigorieff,et al.  Numerical Treatment of Differential Equations , 1978 .

[11]  Jesper Oppelstrup,et al.  The RKFHB4 method for delay — Differential equations , 1978 .

[12]  Herbert Arndt.  Numerical solution of retarded initial value problems: Local and global error and stepsize control , 1984 .

[13]  H. Stetter.  Analysis of Discretization Methods for Ordinary Differential Equations , 1973 .

[14]  K. J. in 't Hout,et al.  A new interpolation procedure for adapting Runge-Kutta methods to delay differential equations , 1992 .

[15]  H. J. Pesch,et al.  Numerical treatment of delay differential equations by Hermite Interpolation , 1981 .

[16]  S. Thompson,et al.  Software for the numerical solution of systems of functional differential equations with state-dependent delays , 1992 .

[17]  Irving R. Epstein,et al.  Delay effects and differential delay equations in chemical kinetics , 1992 .

[18]  M. N. Spijker,et al.  The stability of the θ-methods in the numerical solution of delay differential equations , 1990 .

[19]  C. W. Gear,et al.  Algorithm 407 DIFSUB for Solution of Ordinary Differential Equations , 1971 .

[20]  E. Hairer,et al.  Solving Ordinary Differential Equations I , 1987 .

[21]  P Hogeweg,et al.  Tumor escape from immune elimination: simplified precursor bound cytotoxicity models. , 1985, Journal of theoretical biology.

[22]  Improved Absolute Stability of Predictor-Corrector Methods for Retarded Differential Equations , 1983 .

[23]  Robert D. Skeel.  Analysis of Fixed-Stepsize Methods , 1976 .

[24]  C. W. Geart,et al.  THE EFFECT OF VARIABLE MESH SIZE ON THE STABILITY OF MULTISTEP METHODS , 1974 .

[25]  J. Lambert.  Computational Methods in Ordinary Differential Equations , 1973 .

[26]  A S Perelson,et al.  A mechanism of immune escape by slow-replicating HIV strains. , 1992, Journal of acquired immune deficiency syndromes.

[27]  G. Bocharov,et al.  Mathematical model of antiviral immune response. III. Influenza A virus infection. , 1994, Journal of theoretical biology.

[28]  Richard Bellman,et al.  Differential-Difference Equations , 1967 .

[29]  B. Sulzer,et al.  Memory B Cells Stabilize Cycles in a Repressive Network , 1992 .

[30]  Theodore A. Bickart.  P-stable andP[α, β]-stable integration/interpolation methods in the solution of retarded differential-difference equations , 1982 .

[31]  Theodore A. Bickart,et al.  F-stable and F[α,β]-stable integration/interpolation methods in the solution of retarded differential-difference equations , 1982 .

[32]  Ernst Hairer,et al.  Asymptotic expansions of the global error of fixed-stepsize methods , 1984 .

[33]  M. R. Osborne,et al.  The numerical solution of differential equations , 1961 .

[34]  D. S. Watanabe,et al.  THE STABILITY OF DIFFERENCE FORMULAS FOR DELAY DIFFERENTIAL EQUATIONS , 1985 .