R0_PeriodicODEs

Project information

This repository contains the MATLAB codes for the numerical approximation of the (basic and type) reproduction numbers of linear time-periodic ODE systems. This package allows full flexibility in splitting the model coefficients into ''birth'' and ''transition'' processes. R0_PeriodicODEs computes the pseudospectral discretization of the birth and transition operators B and M, namely BN and MN, obtained from the original time-periodic ODE system by applying the theory of evolution semigroups. Then, the approximation of the relevant reproduction number is obtained through the spectral radius of the approximating ''Next-Generation matrix'' BN*MN^(-1).