Published: 11 June 2021

Modeling the effect of distancing and wearing of face masks on transmission of Covid-19 infection dynamics

Nurudeen O. Lasisi1
Kolawole A. Adeyemo2
1Department of Statistics, Federal Polytechnic, Kaura Namoda, Nigeria
2Department of Computer Science and Mathematics, Nigeria Police Academy, Wudil, Nigeria
Corresponding Author:
Nurudeen O. Lasisi
Views 80
Reads 37
Downloads 851

Abstract

The COVID-19 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 Covid-19 infection dynamics. The invariant region of the model is established. The Covid-19 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, Covid-19 free equilibrium state is locally asymptotically stable in feasible region Ω if R0<1 and globally asymptomatically stable if R0<1, otherwise unstable if R0>1. More so, numerical analysis and simulations of the dynamics of the Covid-19 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.

1. Introduction

The Covid-19 disease is an infectious disease caused by a newly revealed Corona virus [1]. Covid-19 characterised as asymptomatic or mild symptom to severe illness and mortality followed [2-5]. Corona viruses are large group of virus that cause mild to moderate high-respiratory 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 Covid-19 virus experienced mild to non-extreme respiratory illness and recover without requiring particular treatment [7-10]. 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 Covid-19 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 [13-17]. As at now, there are no specific vaccinations or treatments for Covid-19. 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 Monkey-Pox 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 Covid-19 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 Covid-19 pandemic. The reality of Covid-19 was not considered as a compartment in SEIR model. In this case, our model considered the reality and mechanism of Covid-19 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 Covid-19 infection dynamics. The proposed model subdivides the population of interest into five epidemiological compartments. In the modelling the Covid-19 dynamics, the population is divided into five subgroups: susceptible Xt, infective It, contact tracing Gt, isolation Jt, and recovery Rt. The susceptible are individuals that have not contracted the Covid-19 disease; Infective are individuals that have contracted the Covid-19 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. 1Schematic representation for the model equations

Schematic representation for the model equations

The transfer rates between the sub-groups indicating with tick arrows are collection of several epidemiological parameters. The susceptible population increase by recruitment rate Λ and having the natural death rate μ, Infection or Attack rate α is the rate at which Xt is infected, σ is the compliance of face mask usage and ε is the practicing of the social and physical distancing, ω is the unconfirmed contract tracing rate, ϕ is the uninfected contact tracing return back to susceptible and (ϕ-1) is the proportion of contact tracing enter Isolation as a result of infection, It enter Jt with the rate of δ and Jt enter Rt with the recovery parameters γ and τ. The rate of death due to Infection without isolation is μI and Death rate due to Isolation is μJ.

Based on the above schematic representation and assumptions of the models, the equations governing the dynamics of the Covid-19 are given as:

1
dXdt=Λ+ϕG-eαIXN-ωX-μX,
2
dIdt=eαIXN-δ1I-μI-μII,
3
dGdt=ωX-ϕG-1-ϕαIGN-μG,
4
dJdt=δ1I+1-ϕαIGN-γ+τJ-μ+μJJ,
5
dRdt=γJ+τJ-μR,

where, e=(1-ε)(1-σ).

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:

6
dNdt=Λ+ϕG-eαIXN-ωX-μX+eαIXN-δ1I-μI-μII+ωX-ϕG
-1-ϕαIGN-μG+δ1I+1-ϕαIGN-γ+τJ-μ+μJJ+γJ+τJ-μR.

Eq. (6) gives as:

7
dNdt=Λ-μN-μII-μJJ.

In the absence of mortality due to Covid-19 infection (μI=μJ=0), Eq. (7) yields:

8
dNdtΛ-μN.

It follows from Eq. (8), the Gronwall inequality, that:

9
NtN0e-μt+Λμ1-e-μt.

N(t)Λ/μ if N(0)Λ/μ. Thus, Ω is positively invariant. Hence, it is sufficient to consider the model dynamics Eqs. (1-5) in Ω.

Theorem 1: (Invariant Region) the following biological feasible region of the model Eqs. (1-5), Ω=X,I,G,J,RR5:Xt+It+Gt+Jt+RtΛμ is positive invariant and attracting.

3.2. Covid-19 free equilibrium (CFE)

The Covid-19 free equilibrium denoted by E0 of the model system Eqs. (1-5) was obtained as:

10
X0,I0,G0,J0,R0=Λϕ+μω+μϕ+μ-ωϕ,0,0,ωΛϕ+μω+μ-ωϕ,0,
11
N0=Λ(ϕ+μ)+ωΛ(ω+μ)(ϕ+μ)-ωϕ.

3.3. The basic reproduction number, R0

Reproduction number R0 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, R0 of the model system Eqs. (1-5) is obtained by using the next generation operator method, approach of [21] and is given below:

12
R0=(1-ε)(1-σ)α(ϕ+μ)(δ1+μ+μI)[(ϕ+μ)+ω].

3.4. Local stability of Covid-19 free equilibrium

To determine the local stability of Covid-19 free equilibrium of the model system Eqs. (1-5) corresponding to Covid-19 free E0 is obtained.

Theorem 2: If R0<1, then the Covid-19 free equilibrium E0 of the model system Eqs. (1-5) is locally asymptotically stable.

Proof: We used the linearization technique to determine the local stability of the Covid-19 free equilibrium E0, we have:

13
H0=J(X0,I0,G0,J0,R0) =-eαIN-ω-μ-eαXNϕ00eαINeαXN-δ1-μ-μI000ω-(1-ϕ)αGN-ϕ-(1-ϕ)αIN-μ000δ1+(1-ϕ)αGN(1-ϕ)αIN(-γ-τ-μ-μJ)0000γ+τ-μ.

At Covid-19 free equilibrium, we use Gaussian elimination row operation on Eq. (13) and the characteristic equation J(E0)-λI is given as:

14
J(E0)-λI=-P1-λ-eαX0N0ϕ000eαX0-N0P2N0-λ00000ϕω-P1P3P1-λ00000-P4-λ00000-μ-λ,

where, P1=ω+μ, P2=δ1+μ+μI, P3=ϕ+μ, P4=γ+τ+μ+μJ.

The characteristics equation corresponding to H0 is:

15
J(E0)-λI=-P1-λeαX0-N0P2N0-λϕω-P1P3P1-λ(-P4-λ)(-μ-λ)=0,
λ1=-P1=-(ω+μ)<0,
λ2=eαX0-N0P2N0=eαX0-N0(δ1+μ+μI)N0,
λ3=ϕω-(ω+μ)(ϕ+μ)(ω+μ)<0, ω+μϕ+μ>ϕω,
λ4=-P4=-(γ+τ+μ+μJ)<0,
λ5=-μ<0.

To confirm that λ2<0, we have:

16
λ2=eαX0-N0(δ1+μ+μI)N0 =eαΛ(ϕ+μ)(ω+μ)(ϕ+μ)-ωϕ-Λ(ϕ+μ)+ωΛ(ω+μ)(ϕ+μ)-ωϕ(δ1+μ+μI)Λ(ϕ+μ)+ωΛ(ω+μ)(ϕ+μ)-ωϕ,
λ2=[Λ(ϕ+μ)+ωΛ](δ1+μ+μI)[R0-1]Λ(ϕ+μ)+ωΛ,
λ2=δ1+μ+μIR0-1<0, R0<1.

Therefore, Covid-19 Free Equilibrium is locally asymptomatically stable if R0<1 and (ω+μ)(ϕ+μ)>ϕω. So that infections do not persist in the population and under this condition the endemic equilibrium does not exist. It is unstable for R0>1 and endemic equilibrium exists and the Covid-19 infection is maintained in the population.

3.5. Global stability of Covid-19 free equilibrium

To prove global stability of the Covid-19 free equilibrium E0 of the model Eqs. (1-5), we apply the approach by Castillo-Chavez (2002), by rewriting the model system Eqs. (1-5) as:

17
dK1dt=FK1,K2,
dK2dt=HK1,K2, H(K1,0)=0,

where K1=(X0,G0,R0) andK2=(I0,J0), with the components of K1R3 denoting the uninfected groups and the components of K2R2 denoting the infected groups

The Covid-19 free equilibrium is now denoted as E0=(K1*,0), where K1*=(N0,0).

The two Conditions are:

– L1: dK1dt=F(K1,0).

– L2: H^(K1,K2)=AK2-H(K1,K2).

Theorem 3: If R0<1, the Covid-19 free equilibrium E0 of the model system Eqs. (1-5) is globally asymptotically stable (GAS) in Ω.

Proof: To establish the global stability of Covid-19 free equilibrium E0, the two conditions (L1) and (L2) as in Castillo-Chavez (2002) must hold for R0<1.

Now, for the first condition, that is global asymptotic stability of K10, gives:

18
dK1dt=FK1,0=Λ+ϕG0-ωX0-μX0ωX0-ϕG0-μG0-μR0.

Solving linear differential equations of Eq. (18), we have:

19
X0t=Λ+ϕG0ω+μ-Λ+ϕG0ω+μe-ω+μt+X00e-ω+μt,
G0t=ωX0ϕ+μ-ωX0ϕ+μe-ϕ+μt+G00e-ϕ+μt,
R0t=R00e-μt.

Now, clearly Eq. (10) gives [X0(t)+,G0(t)+,R0(t)]N0(t) as t regardless of the value of X0(t), G0(t) and R0(t). Thus, K1*=(N0,0) is globally asymptotically stable.

Next, for the second condition, that is H^(K1,K2)=AK2-H(K1,K2), gives:

20
A=eαX0N0-(δ1+μ+μI)0δ1+(1-ϕ)αG0N0-(γ+τ+μ+μJ),
HK1,K2=eαI0X0N0-δ1+μ+μII0δ1I0+1-ϕαI0G0N0-γ+τ+μ+μJJ0.

Therefore:

21
H^(K1,K2)=AK2-H(K1,K2)
=eαX0N0-(δ1+μ+μI)0δ1+(1-ϕ)αG0N0-(γ+τ+μ+μJ)I0J0
-eαI0X0N0-(δ1+μ+μI)I0δ1I0+(1-ϕ)αI0G0N0-(γ+τ+μ+μJ)J0
=eαX0I0N0-(δ1+μ+μI)I00δ1+(1-ϕ)αI0G0N0-(γ+τ+μ+μJ)J0
-eαI0X0N0-(δ1+μ+μI)I0δ1I0+(1-ϕ)αI0G0N0-(γ+τ+μ+μJ)J0
=00H^K1,K2=AK2-HK1,K2=00,

i.e H^(K1,K2)=0,0T.

It is thus obvious that H^(K1,K2)=0. Hence, the proof is complete. Therefore, Covid-19 free equilibrium is globally asymptomatically Stable in the feasible region Ω if R0<1.

4. Numerical simulations

According to Nigeria Centre for Disease Control (NCDC) reported on 17th June, 2020 on Covid-19 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 = Xt+It+Jt+Gt+Rt, where, Xt= 200,000,000 – [7,735 + 10,000] – 200,000 – (5,967 + 469) = 199,775,829.

Recovery rate (γ) = number of recovery / number of cases.

Infection or attack rate (α) = number of new cases in population at risk / number of persons at risk in the population.

Disease death rate = number of death from Covid-19 infection/ total number of cases death rate due to infection (μI).

Contact tracing rate (ω) = number of tracing or sample tested / total population of the area.

Uninfected contact tracing rate (ϕω) = number of contact tracing tested but uninfected / total number tested for Covid-19.

Table 1Values for variables and parameters of the model

Parameters
Definitions
Values
Λ
Recruitment rate
271.23/day
σ
Usage of face mask
{0, 1}
ε
Practicing social and physical distancing
{0, 1}
ϕ
Uninfected contact tracing rate
0.83/day
ω
Contact tracing rate
0.00052/day
α
Infection or attack rate
0.001/day
γ
Recovery rate
0.336/day
τ
Booster
{0,0.5}
δ
Movement rate from I to J
Assumed (0.6)
μ
Natural death rate
0.0017/day
μI
Death rate due to infection
0.0212/day
μJ
Death rate at isolation centre
0.0053/day
N
Total population number
200,000,000

The data for numerical simulations of R0 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. 2-3.

Fig. 2 shows the effect of face masks usage in reducing the spread of Covid-19 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 Covid-19 infection without isolation, we found that as death rate due to Covid-19 infection without isolation increases, it decreased the reproduction number.

Fig. 2Effect of face masks usage on reproduction number on reproduction number (R0)

Effect of face masks usage on reproduction number on reproduction number (R0)

Fig. 3Effect of uninfected contact tracing on reproduction number (R0)

Effect of uninfected contact tracing on reproduction number (R0)

Fig. 4Effect of infection rate on reproduction number (R0)

Effect of infection rate  on reproduction number (R0)

Fig. 5Effect of rate of I to J on reproduction number (R0)

Effect of rate of I to J  on reproduction number (R0)

Fig. 6Variation of proportion of contact tracing for different values of ω, ε

Variation of proportion of contact tracing for different values of ω, ε

a)

Variation of proportion of contact tracing for different values of ω, ε

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. 7Population of isolation for different values of treatment,τ

Population of isolation  for different values of treatment,τ

Fig. 8Recovery individual population for different values of γ

Recovery individual population  for different values of γ

Fig. 9Variation of proportion of total population in different classes (infected, contact tracing and isolation)

Variation of proportion of total population in different classes (infected, contact tracing and isolation)

Fig. 10Variation of total population in different classes (infected, recovered and isolation)

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 Covid-19 infection dynamics in a population. The invariant region of the model was established. The Covid-19 free equilibrium state was obtained and the stability of Covid-19 free equilibrium was investigated. The model showed that the Covid-19 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 Covid-19 free equilibrium state was established, by the construction of Lyapunov method, it was showed that Covid-19 free equilibrium is globally asymptomatically stable in the feasible region Ω if R0<1. Based on the results of the study and simulations, we conclude that the most effective way to reduce the transmission of Covid-19 Infection is to educate individuals on the use of face masks and sustain physical distancing, because use of face mask reduced the spread of Covid-19 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 Covid-19 infection.

References

  • Y. Li, et al., “Mathematical modeling and epidemic prediction of COVID-19 and its significance to epidemic prevention and control measures,” Annals of Infectious Disease and Epidemiology, Vol. 5, No. 1, 2020, p. 1052.
  • 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/s15010-020-01401-y
  • Q. Li, et al., “Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia,” New England Journal of Medicine, Vol. 382, No. 13, pp. 1199-1207, 2020, https://doi.org/10.1056/NEJMoa2001316
  • C. Rothe, et al. “Transmission of 2019-nCoV 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
  • A. Sahin, et al. “2019 novel coronavirus (COVID-19) 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
  • “WHO statement regarding cluster of pneumonia cases in Wuhan, China,” 2020, https://www.who.int/china/news/detail/09-01-2020-\who-statement-regarding-cluster-of-pneumonia-cases-in-wuhan-china
  • J. M. Read, J. R. E. Bridgen, D. A. T. Cummings, A. Ho, and C. P. Jewell, “Novel coronavirus 2019-nCoV: 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
  • B. Tang, et al., “Estimation of the transmission risk of 2019-nCoV and its implication for public health interventions”, Journal of Clinical Medicine, Vol. 9, No. 2, 2020, p. 462, https://doi.org/10.3390/jcm9020462
  • N. Imai, et al., “Report 3: Transmissibility of 2019-nCoV”, Reference Source, 2020, https://www.imperial.ac.uk/mrc-global-infectious-disease-analysis/news--wuhan-coronavirus
  • 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.
  • “WHO novel coronavirus (2019-nCoV) situation reports”, https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports
  • “Centers for Disease Control and Prevention: 2019 novel coronavirus”, https://www.cdc.gov/coronavirus/2019-ncov
  • 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
  • L. E. Gralinski, and V. D. Menachery, “Return of the coronavirus: 2019-nCoV,” Viruses, Vol. 12, No. 12, 2020, p. 135, https://doi.org/10.3390/v12020135
  • 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
  • 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
  • J.-T. Wu, K. Leung, and G. M Leung, “Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study”, Lancet, Vol. 395, pp. 689–697, 2020, https://doi.org/10.1016/S0140-6736(20)30260-9
  • 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
  • N. O. Lasisi, “Effect of public awareness, behaviours and treatment on infection-age-structured 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
  • N. O. Lasisi, N. I. Akinwande, and F. A. Oguntolu, “Development and exploration of a Mathematical Model for Transmission of Monkey-Pox in Humans”, Journal of Mathematical Models in Engineering, Vol. 6, No. 1, pp. 23–33, 2020, https://doi.org/10.21595/mme.2019.21234
  • 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
  • H. T. Alemneh and G. T. Tilahun, “Mathematical modeling and optimal control analysis of COVID-19 in Ethiopia,” medRxiv, July 24, 2020, https://doi.org/10.1101/2020.07.23.20160473

About this article

Received
04 April 2021
Accepted
04 May 2021
Published
11 June 2021
Keywords
distancing
reproduction number
wearing of mask
local and global stability
modelling
transmission