Modeling the effect of distancing and wearing of face masks on transmission of Covid19 infection dynamics
Nurudeen O. Lasisi^{1} , Kolawole A. Adeyemo^{2}
^{1}Department of Statistics, Federal Polytechnic, Kaura Namoda, Nigeria
^{2}Department of Computer Science and Mathematics, Nigeria Police Academy, Wudil, Nigeria
^{1}Corresponding author
Journal of Complexity in Health Sciences, Vol. 4, Issue 1, 2021, p. 1020.
https://doi.org/10.21595/chs.2021.21976
Received 4 April 2021; received in revised form 25 April 2021; accepted 4 May 2021; published 11 June 2021
JVE Conferences
The COVID19 is an infection caused by corona virus which has been designated as pandemic in the word. In this paper, we proposed a model to study the effect of distancing and wearing of mask on transmission of Covid19 infection dynamics. The invariant region of the model is established. The Covid19 free equilibrium and the reproduction number of the model were obtained. The local and global stability of the model is determined using linearization technique method and Lyapunov method. It was found that, Covid19 free equilibrium state is locally asymptotically stable in feasible region $\mathrm{\Omega}$ if ${R}_{0}<1$ and globally asymptomatically stable if${R}_{0}1$, otherwise unstable if ${R}_{0}>1$. More so, numerical analysis and simulations of the dynamics of the Covid19 infection are presented. It was found that as effectiveness of face mask usage increased it resulted in a decrease in secondary cases and as physical distancing is maintained shown reduction in reproduction number. It is also found that, as infection rate increases, it translated to increase in secondary cases and as number of isolation individual increases, it decreased the reproduction number. More so, varying the treatment rate, means as treatment rate increases, it shown that an individuals in the isolation decreases.
Keywords: distancing, reproduction number, wearing of mask, local and global stability, modelling, transmission.
1. Introduction
The Covid19 disease is an infectious disease caused by a newly revealed Corona virus [1]. Covid19 characterised as asymptomatic or mild symptom to severe illness and mortality followed [25]. Corona viruses are large group of virus that cause mild to moderate highrespiratory tract illness, the main symptom are, fever, coughing, shortness of breath, trouble, breathing, fatigue, chills and shaking, body aches, headache, sore throat, loss of smell or taste, nausea, diarrhea [6]. The individuals infected with Covid19 virus experienced mild to nonextreme respiratory illness and recover without requiring particular treatment [710]. People with ages and those with underlying health problems such as diabetes illness, cancer, and cardiovascular disease are more likely to develop serious illness [11, 12]. The most excellent way to prevent and slow down transmission is be well informed and educating about the Covid19 virus, the disease it causes and how it spreads, it found that for people with mild disease, recover time is about two weeks, while people with severe or critical disease recover within three to six weeks [1317]. As at now, there are no specific vaccinations or treatments for Covid19. Meanwhile, there are numerous ongoing clinical trials evaluating potential treatments.
Mathematical models have played an indispensable role for understanding the dynamics of disease. Among the most popular are Ebola Virus [18] and The Human Immuno deficiency Virus (HIV) [19]. These models have been helpful to study the control and the transmission of the virus kinetics in order to provide a quantitative understanding and create public awareness of the virus, while [20] have designed mathematical models to explore the transmission dynamics and control of the MonkeyPox infection. Hence, with the scarcity of data, the calculation, assumption and estimation of the values of the parameters of the model would be considered.
XGIJR model is an extension of standard SEIR model which we considered as more suitable for Covid19 infection. In SEIR model, E is latent compartment denoted the individuals who are infected but not infectious. XGIJR model is the modification of SEIR model. In [21], the authors used SEIR model equations by incorporating time depend control and optimal control analysis. They used SEIR compartments to considered limited parameters, because they only focused at the trends, feeling, economic and political impacts caused by Covid19 pandemic. The reality of Covid19 was not considered as a compartment in SEIR model. In this case, our model considered the reality and mechanism of Covid19 infection, which are Susceptible to the infection (S), Contact tracing (G) of individual, Infected individual (I), diagnosed and isolation of individual (J), and Recovered individual (R) and parameters such as, effectiveness of Mask usage and physical distancing. Even though there are still limitations on our work due to availability of data. Thus, give room for other to research and extend on.
2. Model formulation
A nonlinear mathematical model is proposed and analysed to study the effect of distancing and wearing of mask on transmission of Covid19 infection dynamics. The proposed model subdivides the population of interest into five epidemiological compartments. In the modelling the Covid19 dynamics, the population is divided into five subgroups: susceptible $X\left(t\right)$, infective $I\left(t\right)$, contact tracing $G\left(t\right)$, isolation $J\left(t\right)$, and recovery $R\left(t\right)$. The susceptible are individuals that have not contracted the Covid19 disease; Infective are individuals that have contracted the Covid19 infection; Contact Tracing are individuals that have contact with infected person but have not being confirm of the infection; Isolation are individuals that have been confirmed of the infection and move to Isolation in order not to spread the disease, as shown in Fig. 1.
Fig. 1. Schematic representation for the model equations
The transfer rates between the subgroups indicating with tick arrows are collection of several epidemiological parameters. The susceptible population increase by recruitment rate $\mathrm{\Lambda}$ and having the natural death rate $\mu $, Infection or Attack rate $\alpha $ is the rate at which $X\left(t\right)$ is infected, $\sigma $ is the compliance of face mask usage and $\epsilon $ is the practicing of the social and physical distancing, $\omega $ is the unconfirmed contract tracing rate, $\varphi $ is the uninfected contact tracing return back to susceptible and ($\varphi 1$) is the proportion of contact tracing enter Isolation as a result of infection, $I\left(t\right)$ enter $J\left(t\right)$ with the rate of $\delta $ and $J\left(t\right)$ enter $R\left(t\right)$ with the recovery parameters $\gamma $ and $\tau $. The rate of death due to Infection without isolation is ${\mu}_{I}$ and Death rate due to Isolation is ${\mu}_{J}$.
Based on the above schematic representation and assumptions of the models, the equations governing the dynamics of the Covid19 are given as:
where, $e=(1\epsilon )(1\sigma )$.
3. Basic properties of the model
3.1. Invariant region
To obtain the invariant region, we considered the total human population ($N$), where $N=X+I+G+Q+R$. Then, the differentiation of N with respect to time leading to:
$\left(1\varphi \right)\frac{\alpha IG}{N}\mu G+{\delta}_{1}I+\left(1\varphi \right)\frac{\alpha IG}{N}\left(\gamma +\tau \right)J\left(\mu +{\mu}_{J}\right)J+\gamma J+\tau J\mu R.$
Eq. (6) gives as:
In the absence of mortality due to Covid19 infection (${\mu}_{I}={\mu}_{J}=0$), Eq. (7) yields:
It follows from Eq. (8), the Gronwall inequality, that:
$N\left(t\right)\le \mathrm{\Lambda}/\mu $ if $N\left(0\right)\le \mathrm{\Lambda}/\mu $. Thus, $\mathrm{\Omega}$ is positively invariant. Hence, it is sufficient to consider the model dynamics Eqs. (15) in $\mathrm{\Omega}$.
Theorem 1: (Invariant Region) the following biological feasible region of the model Eqs. (15), $\mathrm{\Omega}=\left\{\left(X,I,G,J,R\right)\in {\mathfrak{R}}^{5}:X\left(t\right)+I\left(t\right)+G\left(t\right)+J\left(t\right)+R\left(t\right)\le \frac{\mathrm{\Lambda}}{\mu}\right\}$ is positive invariant and attracting.
3.2. Covid19 free equilibrium (CFE)
The Covid19 free equilibrium denoted by ${E}_{0}$ of the model system Eqs. (15) was obtained as:
3.3. The basic reproduction number, ${\mathit{R}}_{0}$
Reproduction number ${R}_{0}$ is a process of measuring the transmission possible of disease. It is an average number of secondary infections produced by individual case of an infection in a population where everyone is susceptible. The basic reproduction number, ${R}_{0}$ of the model system Eqs. (15) is obtained by using the next generation operator method, approach of [21] and is given below:
3.4. Local stability of Covid19 free equilibrium
To determine the local stability of Covid19 free equilibrium of the model system Eqs. (15) corresponding to Covid19 free ${E}_{0}$ is obtained.
Theorem 2: If ${R}_{0}<1$, then the Covid19 free equilibrium ${E}_{0}$ of the model system Eqs. (15) is locally asymptotically stable.
Proof: We used the linearization technique to determine the local stability of the Covid19 free equilibrium ${E}_{0}$, we have:
At Covid19 free equilibrium, we use Gaussian elimination row operation on Eq. (13) and the characteristic equation $\leftJ\left({E}_{0}\right)\lambda I\right$ is given as:
where, ${P}_{1}=\omega +\mu $, ${P}_{2}={\delta}_{1}+\mu +{\mu}_{I}$, ${P}_{3}=\varphi +\mu $, ${P}_{4}=\gamma +\tau +\mu +{\mu}_{J}$.
The characteristics equation corresponding to ${H}_{0}$ is:
${\lambda}_{1}={P}_{1}=(\omega +\mu )<0,$
${\lambda}_{2}=\frac{e\alpha {X}^{0}{N}^{0}{P}_{2}}{{N}^{0}}=\frac{e\alpha {X}^{0}{N}^{0}({\delta}_{1}+\mu +{\mu}_{I})}{{N}^{0}},$
To confirm that ${\lambda}_{2}<0$, we have:
${\lambda}_{2}=\frac{\left[\mathrm{\Lambda}\right(\varphi +\mu )+\omega \mathrm{\Lambda}]({\delta}_{1}+\mu +{\mu}_{I})[{R}_{0}1]}{\mathrm{\Lambda}(\varphi +\mu )+\omega \mathrm{\Lambda}},$
${\lambda}_{2}=\left({\delta}_{1}+\mu +{\mu}_{I}\right)\left({R}_{0}1\right)<0,\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{R}_{0}1.$
Therefore, Covid19 Free Equilibrium is locally asymptomatically stable if ${R}_{0}<1$ and $(\omega +\mu )(\varphi +\mu )>\varphi \omega $. So that infections do not persist in the population and under this condition the endemic equilibrium does not exist. It is unstable for ${R}_{0}>1$ and endemic equilibrium exists and the Covid19 infection is maintained in the population.
3.5. Global stability of Covid19 free equilibrium
To prove global stability of the Covid19 free equilibrium ${E}_{0}$ of the model Eqs. (15), we apply the approach by CastilloChavez (2002), by rewriting the model system Eqs. (15) as:
$\frac{d{K}_{2}}{dt}=H\left({K}_{1},{K}_{2}\right),\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}H({K}_{1},0)=0,$
where ${K}_{1}=({X}^{0},{G}^{0},{R}^{0})$ and${K}_{2}=({I}^{0},{J}^{0})$, with the components of ${K}_{1}\in {\mathfrak{R}}^{3}$ denoting the uninfected groups and the components of ${K}_{2}\in {\mathfrak{R}}^{2}$ denoting the infected groups
The Covid19 free equilibrium is now denoted as ${E}_{0}=({K}_{1}^{\mathrm{*}},0)$, where ${K}_{1}^{\mathrm{*}}=({N}^{0},0)$.
The two Conditions are:
– L1: $\frac{d{K}_{1}}{dt}=F({K}_{1},0)$.
– L2: $\widehat{H}({K}_{1},{K}_{2})=A{K}_{2}H({K}_{1},{K}_{2})$.
Theorem 3: If ${R}_{0}<1$, the Covid19 free equilibrium ${E}_{0}$ of the model system Eqs. (15) is globally asymptotically stable (GAS) in $\mathrm{\Omega}$.
Proof: To establish the global stability of Covid19 free equilibrium ${E}_{0}$, the two conditions (L1) and (L2) as in CastilloChavez (2002) must hold for ${R}_{0}<1$.
Now, for the first condition, that is global asymptotic stability of ${K}_{1}^{0}$, gives:
Solving linear differential equations of Eq. (18), we have:
${G}^{0}\left(t\right)=\frac{\omega {X}^{0}}{\left(\varphi +\mu \right)}\frac{\omega {X}^{0}}{\left(\varphi +\mu \right)}{e}^{\left(\varphi +\mu \right)t}+{G}^{0}\left(0\right){e}^{\left(\varphi +\mu \right)t},$
${R}^{0}\left(t\right)={R}^{0}\left(0\right){e}^{\mu t}.$
Now, clearly Eq. (10) gives $\left[{X}^{0}\right(t)+,{G}^{0}(t)+,{R}^{0}(t\left)\right]\to {N}^{0}\left(t\right)$ as $t\to \infty $ regardless of the value of ${X}^{0}\left(t\right)$, ${G}^{0}\left(t\right)$ and ${R}^{0}\left(t\right)$. Thus, ${K}_{1}^{\mathrm{*}}=({N}^{0},0)$ is globally asymptotically stable.
Next, for the second condition, that is $\widehat{H}({K}_{1},{K}_{2})=A{K}_{2}H({K}_{1},{K}_{2})$, gives:
$H\left({K}_{1},{K}_{2}\right)=\left(\begin{array}{c}\frac{e\alpha {I}^{0}{X}^{0}}{{N}^{0}}\left({\delta}_{1}+\mu +{\mu}_{I}\right){I}^{0}\\ {\delta}_{1}{I}^{0}+\left(1\varphi \right)\frac{\alpha {I}^{0}{G}^{0}}{{N}^{0}}\left(\gamma +\tau +\mu +{\mu}_{J}\right){J}^{0}\end{array}\right).$
Therefore:
$=\left(\begin{array}{ll}\frac{e\alpha {X}^{0}}{{N}^{0}}({\delta}_{1}+\mu +{\mu}_{I})& 0\\ {\delta}_{1}+(1\varphi )\frac{\alpha {G}^{0}}{{N}^{0}}& (\gamma +\tau +\mu +{\mu}_{J})\end{array}\right)\left(\begin{array}{l}{I}^{0}\\ {J}^{0}\end{array}\right)$
$\left(\begin{array}{c}\frac{e\alpha {I}^{0}{X}^{0}}{{N}^{0}}({\delta}_{1}+\mu +{\mu}_{I}){I}^{0}\\ {\delta}_{1}{I}^{0}+(1\varphi )\frac{\alpha {I}^{0}{G}^{0}}{{N}^{0}}(\gamma +\tau +\mu +{\mu}_{J}){J}^{0}\end{array}\right)$
$=\left(\begin{array}{cc}\frac{e\alpha {X}^{0}{I}^{0}}{{N}^{0}}({\delta}_{1}+\mu +{\mu}_{I}){I}^{0}& 0\\ {\delta}_{1}+(1\varphi )\frac{\alpha {I}^{0}{G}^{0}}{{N}^{0}}& (\gamma +\tau +\mu +{\mu}_{J}){J}^{0}\end{array}\right)$
$\left(\begin{array}{c}\frac{e\alpha {I}^{0}{X}^{0}}{{N}^{0}}({\delta}_{1}+\mu +{\mu}_{I}){I}^{0}\\ {\delta}_{1}{I}^{0}+(1\varphi )\frac{\alpha {I}^{0}{G}^{0}}{{N}^{0}}(\gamma +\tau +\mu +{\mu}_{J}){J}^{0}\end{array}\right)$
$=\left(\begin{array}{l}0\\ 0\end{array}\right)\widehat{H}\left({K}_{1},{K}_{2}\right)=A{K}_{2}H\left({K}_{1},{K}_{2}\right)=\left(\begin{array}{l}0\\ 0\end{array}\right),$
i.e $\widehat{H}({K}_{1},{K}_{2})={\left(\mathrm{0,0}\right)}^{T}$.
It is thus obvious that $\widehat{H}({K}_{1},{K}_{2})=0$. Hence, the proof is complete. Therefore, Covid19 free equilibrium is globally asymptomatically Stable in the feasible region $\mathrm{\Omega}$ if ${R}_{0}<1$.
4. Numerical simulations
According to Nigeria Centre for Disease Control (NCDC) reported on 17th June, 2020 on Covid19 updated, total cumulative tested person is 103,799; total cumulative infected cases is 17,735; total cumulative recovered cases is 5,967; total cumulative death is 469 and we assumed the number of the contact tracing is 200,000. We calculate and estimate the parameters values based on the availability of information from the aforementioned organization. The reasons for the values are obviously explained in details in the following statements. According to the available data, 80 % of total death resulted to death due to Disease without quarantine in Nigeria, however, due to the overwhelmed and poor medical infrastructure and 20 % of total death resulted to death due to quarantine.
Total population = $X\left(t\right)+\left[I\left(t\right)+J\left(t\right)\right]+G\left(t\right)+R\left(t\right)$, where, $X\left(t\right)=$200,000,000 – [7,735 + 10,000] – 200,000 – (5,967 + 469) = 199,775,829.
Recovery rate ($\gamma $) = number of recovery / number of cases.
Infection or attack rate ($\alpha $) = number of new cases in population at risk / number of persons at risk in the population.
Disease death rate = number of death from Covid19 infection/ total number of cases death rate due to infection (${\mu}_{I}$).
Contact tracing rate ($\omega $) = number of tracing or sample tested / total population of the area.
Uninfected contact tracing rate ($\varphi \omega $) = number of contact tracing tested but uninfected / total number tested for Covid19.
Table 1. Values for variables and parameters of the model
Parameters

Definitions

Values

$\Lambda $

Recruitment rate

271.23/day

$\sigma $

Usage of face mask

{0, 1}

$\epsilon $

Practicing social and physical distancing

{0, 1}

$\varphi $

Uninfected contact tracing rate

0.83/day

$\omega $

Contact tracing rate

0.00052/day

$\alpha $

Infection or attack rate

0.001/day

$\gamma $

Recovery rate

0.336/day

$\tau $

Booster

{0,0.5}

$\delta $

Movement rate from $I$ to $J$

Assumed (0.6)

$\mu $

Natural death rate

0.0017/day

${\mu}_{I}$

Death rate due to infection

0.0212/day

${\mu}_{J}$

Death rate at isolation centre

0.0053/day

$N$

Total population number

200,000,000

The data for numerical simulations of ${R}_{0}$ with respect to each of the epidemic parameters are given in Table 1 and Using Maple 17 Software for the graphical representation of reproduction numbers and model simulations with parameter values of the model equations are sin Figs. 23.
Fig. 2 shows the effect of face masks usage in reducing the spread of Covid19 infection. It is shown that as the use of face mask increases it is resulted to decrease in reproduction number (secondary cases). Thus, varying the social and physical distancing, we observed that reproduction number reduced significantly.
Fig. 3 showed that as uninfected contract tracing increases, it is also increased the reproduction number, this is because, as the uninfected contact tracing reappear into population, it increased the secondary cases and make individual more susceptible to the infection. Meanwhile, varying the contact rate we found that as contact rate increases, it increased the secondary cases.
It is observed from Fig. 4 that as the infection rate (attack rate) increases, it is translated to increase on reproduction number (secondary cases), also as contact tracing rate increase, the reproduction number decreased.
Fig. 5 shows that as rate of isolation increases, it decreased the reproduction number, varying the death rate due to Covid19 infection without isolation, we found that as death rate due to Covid19 infection without isolation increases, it decreased the reproduction number.
Fig. 2. Effect of face masks usage on reproduction number on reproduction number $\left({R}_{0}\right)$
Fig. 3. Effect of uninfected contact tracing on reproduction number $\left({R}_{0}\right)$
Fig. 4. Effect of infection rate on reproduction number $\left({R}_{0}\right)$
Fig. 5. Effect of rate of $I$ to $J$ on reproduction number $\left({R}_{0}\right)$
Fig. 6. Variation of proportion of contact tracing for different values of $\omega $, $\epsilon $
a)
b)
Fig. 6(a) shows that as the proportion of contract tracing increases, it is translated to decrease in force of infection. It is also observed from Fig. 6(b) that as number of complier in physical distancing increases, it is resulted in decrease in force of infection.
Fig. 7, shown the dynamics of population of isolation, it is observed that the number of isolation individuals increases with time, varying treatment rate. It means as treatment rate increased, it is seen the individuals in the isolation decreases.
Fig. 8 shown the dynamics of Recovered person with time, with different rate, it is observed from the simulation, that as recovery rate increases, it is translate to increase in recovery individual.
Fig. 9 shows the proportion of three populations, it is observed that the proportion of infective population increase to a pick before started decreasing as a result of isolation of infected person. The process of contact tracing is not encouraging as it shown in simulation.
It is observed from Fig. 10, that as infected person reached it pick, the isolation and recovery individual started increasing which resulted to decreased in infected person.
Fig. 7. Population of isolation for different values of treatment,$\tau $
Fig. 8. Recovery individual population for different values of $\gamma $
Fig. 9. Variation of proportion of total population in different classes (infected, contact tracing and isolation)
Fig. 10. Variation of total population in different classes (infected, recovered and isolation)
5. Conclusions
A nonlinear mathematical model has been proposed and analysed to study the effect of distancing and wearing of masks on transmission of Covid19 infection dynamics in a population. The invariant region of the model was established. The Covid19 free equilibrium state was obtained and the stability of Covid19 free equilibrium was investigated. The model showed that the Covid19 free equilibrium is locally stable using linearization technique method at threshold parameter less than unity and unstable at threshold parameter greater than unity. The global stability of Covid19 free equilibrium state was established, by the construction of Lyapunov method, it was showed that Covid19 free equilibrium is globally asymptomatically stable in the feasible region $\mathrm{\Omega}$ if ${R}_{0}<1$. Based on the results of the study and simulations, we conclude that the most effective way to reduce the transmission of Covid19 Infection is to educate individuals on the use of face masks and sustain physical distancing, because use of face mask reduced the spread of Covid19 infection and physical distancing also reduced the number of secondary cases. More so, as isolation rate increased, it depleted the reproduction number and as treatment rate increased, it is seen that the individuals in the isolation decreased. Finally, there is need for further research work on the effect of asymptomatic and symptomatic on the spread of the disease, vaccination and other control strategies should be consider in controlling the spread of Covid19 infection.
References
 Y. Li, et al., “Mathematical modeling and epidemic prediction of COVID19 and its significance to epidemic prevention and control measures,” Annals of Infectious Disease and Epidemiology, Vol. 5, No. 1, 2020, p. 1052. [Search CrossRef]
 Z. J. Cheng and J. Shan, “2019 Novel coronavirus: where we are and what we know”, Infection Journal, Vol. 48, No. 2, pp. 155–163, 2020, https://doi.org/10.1007/s1501002001401y [Publisher]
 Q. Li, et al., “Early transmission dynamics in Wuhan, China, of novel coronavirusinfected pneumonia,” New England Journal of Medicine, Vol. 382, No. 13, pp. 11991207, 2020, https://doi.org/10.1056/NEJMoa2001316 [Publisher]
 C. Rothe, et al. “Transmission of 2019nCoV Infection from an asymptomatic contact in Germany”, New England Journal of Medicine, Vol. 382, No. 10, pp. 970–971, 2020, https://doi.org/10.1056/NEJMc2001468 [Publisher]
 A. Sahin, et al. “2019 novel coronavirus (COVID19) outbreak: A review of the current literature”, Eurasian Journal of Medicine and Oncology, Vol. 4, No. 1, pp. 1–7, 2020, https://doi.org/10.14744/ejmo.2020.12220 [Publisher]
 “WHO statement regarding cluster of pneumonia cases in Wuhan, China,” 2020, https://www.who.int/china/news/detail/09012020\whostatementregardingclusterofpneumoniacasesinwuhanchina [Search CrossRef]
 J. M. Read, J. R. E. Bridgen, D. A. T. Cummings, A. Ho, and C. P. Jewell, “Novel coronavirus 2019nCoV: early estimation of epidemiological parameters and epidemic predictions”, medRxiv, Vol. 1, No. 1, 2020, p. 234, https://doi.org/10.1101/2020.01.23.20018549 [Publisher]
 B. Tang, et al., “Estimation of the transmission risk of 2019nCoV and its implication for public health interventions”, Journal of Clinical Medicine, Vol. 9, No. 2, 2020, p. 462, https://doi.org/10.3390/jcm9020462 [Publisher]
 N. Imai, et al., “Report 3: Transmissibility of 2019nCoV”, Reference Source, 2020, https://www.imperial.ac.uk/mrcglobalinfectiousdiseaseanalysis/newswuhancoronavirus [Search CrossRef]
 H. Zhu, et al., “Host and infectivity prediction of Wuhan 2019 novel coronavirus using deep learning algorithm”, Mathematical Biosciences and Engineering Vol. 17, No. 3, pp. 2708–2724, 2020. [Search CrossRef]
 “WHO novel coronavirus (2019nCoV) situation reports”, https://www.who.int/emergencies/diseases/novelcoronavirus2019/situationreports [Search CrossRef]
 “Centers for Disease Control and Prevention: 2019 novel coronavirus”, https://www.cdc.gov/coronavirus/2019ncov [Search CrossRef]
 P. Zhou, et al., “Discovery of a novel coronavirus associated with the recent pneumonia outbreak in humans and its potential bat origin”, bioRxiv, 2020, https://doi.org/10.1101/2020.01.22.914952 [Publisher]
 L. E. Gralinski, and V. D. Menachery, “Return of the coronavirus: 2019nCoV,” Viruses, Vol. 12, No. 12, 2020, p. 135, https://doi.org/10.3390/v12020135 [Publisher]
 V. J. Munster, M. Koopmans, N. V. Doremalen, D. V. Riel, and E. D. Wit, “A novel coronavirus emerging in China – key questions for impact assessment”, New England Journal of Medicine, Vol. 382, pp. 692–694, 2020, https://doi.org/10.1056/NEJMp2000929 [Publisher]
 C. Geller, M. Varbanov, and R. E. Duval, “Human coronaviruses: Insights into environmental resistance and its influence on the development of new antiseptic strategies”, Viruses, Vol. 4, pp. 3044–3068, 2012, https://doi.org/10.3390/v4113044 [Publisher]
 J.T. Wu, K. Leung, and G. M Leung, “Nowcasting and forecasting the potential domestic and international spread of the 2019nCoV outbreak originating in Wuhan, China: a modelling study”, Lancet, Vol. 395, pp. 689–697, 2020, https://doi.org/10.1016/S01406736(20)302609 [Publisher]
 N. O. Lasisi et. al., “Mathematical model for Ebola virus infection in Human with effectiveness of drug usage”, Journal of Applied Science Environmental Management, Vol. 22, No. 7, pp. 1089–1095, 2018, https://dx.doi.org/10.4314/jasem.v22i7.16 [Publisher]
 N. O. Lasisi, “Effect of public awareness, behaviours and treatment on infectionagestructured of mathematical model for HIV/AIDS dynamics”, Journal of Mathematical Models in Engineering, Vol. 6, No. 2, p. 103–121, 2020, https://doi.org/10.21595/mme.2020.21249 [Publisher]
 N. O. Lasisi, N. I. Akinwande, and F. A. Oguntolu, “Development and exploration of a Mathematical Model for Transmission of MonkeyPox in Humans”, Journal of Mathematical Models in Engineering, Vol. 6, No. 1, pp. 23–33, 2020, https://doi.org/10.21595/mme.2019.21234 [Publisher]
 N. O. Lasisi, N. I. Akinwande, and S. Abdulrahaman, “Optimal Control and Effect of Poor Sanitation on Modeling the Acute Diarrhea Infection”, Journal of Complexity in Health Sciences, Vol. 3, No. 1, pp. 91–103, 2020, https://doi.org/10.21595/chs.2020.21409 [Publisher]
 H. T. Alemneh and G. T. Tilahun, “Mathematical modeling and optimal control analysis of COVID19 in Ethiopia,” medRxiv, July 24, 2020, https://doi.org/10.1101/2020.07.23.20160473 [Publisher]