Numerical bifurcation analysis of renewal equations via pseudospectral approximation

We propose an approximation of nonlinear renewal equations by means of ordinary differential equations. We consider the integrated state, which is absolutely continuous and satisfies a delay differential equation. By applying the pseudospectral approach to the abstract formulation of the differential equation, we obtain an approximating system of ordinary differential equations. We present convergence proofs for equilibria and the associated characteristic roots, and we use some models from ecology and epidemiology to illustrate the benefits of the approach to perform numerical bifurcation analyses of equilibria and periodic solutions. The numerical simulations show that the implementation of the new approximating system is ten times more efficient than the one originally proposed in [Breda et al, SIAM Journal on Applied Dynamical Systems, 2016], as it avoids the numerical inversion of an algebraic equation.

[1]  L. Trefethen Spectral Methods in MATLAB , 2000 .

[2]  Anke Schmid,et al.  An Introduction To Complex Analysis , 2016 .

[3]  Odo Diekmann,et al.  Twin semigroups and delay equations , 2019, Journal of Differential Equations.

[4]  Odo Diekmann,et al.  Pseudospectral Approximation of Hopf Bifurcation for Delay Differential Equations , 2020, SIAM J. Appl. Dyn. Syst..

[5]  M. Gyllenberg,et al.  Stability analysis of a state-dependent delay differential equation for cell maturation: analytical and numerical methods , 2019, Journal of Mathematical Biology.

[6]  G. Samaey,et al.  DDE-BIFTOOL v. 2.00: a Matlab package for bifurcation analysis of delay differential equations , 2001 .

[7]  Dimitri Breda,et al.  Stability of Linear Delay Differential Equations: A Numerical Approach with MATLAB , 2014 .

[8]  Odo Diekmann,et al.  Stability and Bifurcation Analysis of Volterra Functional Equations in the Light of Suns and Stars , 2007, SIAM J. Math. Anal..

[9]  Dimitri Breda,et al.  Pseudospectral Differencing Methods for Characteristic Roots of Delay Differential Equations , 2005, SIAM J. Sci. Comput..

[10]  O. Diekmann,et al.  Pseudospectral discretization of delay differential equations in sun-star formulation: Results and conjectures , 2020, Discrete & Continuous Dynamical Systems - S.

[11]  Kuan Xu,et al.  The Chebyshev points of the first kind , 2016 .

[12]  Dimitri Breda,et al.  Pseudospectral Discretization of Nonlinear Delay Equations: New Prospects for Numerical Bifurcation Analysis , 2016, SIAM J. Appl. Dyn. Syst..

[13]  D. Breda,et al.  Collocation Techniques for Structured Populations Modeled by Delay Equations , 2020 .

[14]  Lloyd N. Trefethen,et al.  Barycentric Lagrange Interpolation , 2004, SIAM Rev..

[15]  P. Revesz Interpolation and Approximation , 2010 .

[16]  S. P. Blythe,et al.  Nicholson's blowflies revisited , 1980, Nature.

[17]  Mats Gyllenberg,et al.  Equations with infinite delay: Numerical bifurcation analysis via pseudospectral discretization , 2018, Appl. Math. Comput..

[18]  O. Diekmann,et al.  The Second Half - With a Quarter of a Century Delay , 2008 .

[19]  D. Breda,et al.  Floquet Theory and Stability of Periodic Solutions of Renewal Equations , 2020 .

[20]  Dimitri Breda,et al.  Numerical continuation and delay equations: A novel approach for complex models of structured populations , 2020 .

[21]  A. D. Roos,et al.  PSPManalysis: Steady‐state and bifurcation analysis of physiologically structured population models , 2020, Methods in Ecology and Evolution.

[22]  Odo Diekmann,et al.  Equations with infinite delay: blending the abstract and the concrete , 2012 .

[23]  Y. Kuznetsov,et al.  New features of the software MatCont for bifurcation analysis of dynamical systems , 2008 .

[24]  Giuseppe Mastroianni,et al.  Some new results on Lagrange interpolation for bounded variation functions , 2010, J. Approx. Theory.

[25]  R. Cooke Real and Complex Analysis , 2011 .

[26]  Odo Diekmann,et al.  UvA-DARE ( Digital Academic Repository ) Daphnia revisited : local stability and bifurcation theory for physiologically structured population models explained by way of an example , 2010 .

[27]  O. Diekmann,et al.  Numerical Bifurcation Analysis of Physiologically Structured Population Models via Pseudospectral Approximation , 2020, Vietnam Journal of Mathematics.

[28]  Dimitri Breda,et al.  A numerical approach for investigating the stability of equilibria for structured population models , 2013, Journal of biological dynamics.

[29]  Satish C. Reddy,et al.  A MATLAB differentiation matrix suite , 2000, TOMS.

[30]  O. Diekmann,et al.  Prelude to hopf bifurcation in an epidemic model: Analysis of a characteristic equation associated with a nonlinear Volterra integral equation , 1982, Journal of mathematical biology.

[31]  Dimitri Breda,et al.  Computing the Eigenvalues of Realistic Daphnia Models by Pseudospectral Methods , 2015, SIAM J. Sci. Comput..