Mathematical modelling and optimal cost-effective control of COVID-19 transmission dynamics

The novel coronavirus disease (COVID-19) caused by a new strain of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remains the current global health challenge. In this paper, an epidemic model based on system of ordinary differential equations is formulated by taking into account the transmission routes from symptomatic, asymptomatic and hospitalized individuals. The model is fitted to the corresponding cumulative number of hospitalized individuals (active cases) reported by the Nigeria Centre for Disease Control (NCDC), and parameterized using the least squares method. The basic reproduction number which measures the potential spread of COVID-19 in the population is computed using the next generation operator method. Further, Lyapunov function is constructed to investigate the stability of the model around a disease-free equilibrium point. It is shown that the model has a globally asymptotically stable disease-free equilibrium if the basic reproduction number of the novel coronavirus transmission is less than one. Sensitivities of the model to changes in parameters are explored, and safe regions at certain threshold values of the parameters are derived. It is revealed further that the basic reproduction number can be brought to a value less than one in Nigeria, if the current effective transmission rate of the disease can be reduced by 50%. Otherwise, the number of active cases may get up to 2.5% of the total estimated population. In addition, two time-dependent control variables, namely preventive and management measures, are considered to mitigate the damaging effects of the disease using Pontryagin's maximum principle. The most cost-effective control measure is determined through cost-effectiveness analysis. Numerical simulations of the overall system are implemented in ${\text{MatLab}}^{\circledR }$ for demonstration of the theoretical results.