# Sensitivity and Optimal Control Analysis of an Extended SEIR COVID-19 Mathematical Model

Authors: C. M. Wachira; G. O. Lawi; L. O. Omondi
Journal: Journal of Mathematics
 (2022)

---

## Abstract

In this paper, a mathematical model based on a system of ordinary differential equations is developed with vaccination as an intervention for the transmission dynamics of coronavirus 2019 (COVID-19). The model solutions are shown to be well posed. The vaccine reproduction number is computed by using the next-generation matrix approach. The sensitivity analysis carried out on this model showed that the vaccination rate and vaccine efficacy are among the most sensitive parameters of the vaccine reproduction number,RV. The optimal control problem is solved with the rate of vaccination and the transition rate from the vaccinated class to the infected class as control variables. Finally, the numerical simulations showed that the control intervention should aim to increase the vaccination rate with a high-efficacy vaccine.


---

## Body

## 1. Introduction

COVID-19 is a highly contagious respiratory disease caused by the novel severe acute respiratory syndrome coronavirus–2 (SARS-CoV-2). The increasing prevalence rate of the virus nearly devastated healthcare systems across the world [1, 2]. On March 11, 2020, the World Health Organization (WHO) declared the disease a global health pandemic. It is transmitted through inhalation of respiratory droplets from an infectious person emitted through sneezing, coughing, and having a close conversation, as well as contact with contaminated surfaces [3]. The infection prominently affects the respiratory tract. The initial COVID-19 symptoms are fever, dry cough, and fatigue. Other symptoms include body aches and pains, a sore throat, diarrhoea, headache, loss of taste and smell, breathing difficulties, skin rash, and discoloration of the fingers and toes [3].COVID-19 morbidity and mortality can be compounded by pre-existing diseases such as cardiovascular disease, respiratory disease, cancer, infectious diseases, and substance abuse. Furthermore, pre-existing conditions such as environmental, demographic, and socioeconomic factors may influence the COVID-19 incidence rate [4]. To prevent and control the spread of COVID-19, health providers and governments around the world implemented containment measures such as lockdowns, travel bans, cessation of movement, social distancing, proper hygiene, and proper hand washing, among others. Despite these safeguards, the virus continues to spread, although at a slow rate.Vaccination has been a major public health tool in modern medicine used to minimize the impact of many infectious diseases on humans [5, 6]. A vaccine is any biologically derived substance that elicits a protective immune response when administered to a susceptible host. Vaccines help the body to prepare for disease by taking the advantage of the fact that immunity knows how to defend against infectious organisms, which are typically viruses, bacteria, or toxins [7]. Given the emergence of highly transmissible new variants, decreased vaccine effectiveness, and unequal vaccine availability, there is a growing concern that vaccination may not result in herd immunity. Vaccines reduce the risk of severe disease and death; therefore, a country’s vaccine coverage must be high in order to protect the healthcare system from an infection surge [2]. Since late 2020, several SARS-CoV-2 variants of concern that are more transmissible and/or vaccine-resistant have been identified [8].In modern medicine, vaccination has been an important public health technique for reducing the burden of many human infectious diseases [5, 6, 9]. Several COVID-19 vaccine candidates have been developed and are already in use [10]. Vaccines such as Pfizer, Moderna, AstraZeneca, Sputnik, Sinopharm, Johnson & Johnson, and many others are currently being administered around the world to date. It is a social responsibility choice of an individual to get vaccinated. However, offsetting the fear and adverse side effects of the vaccine has to be carried out[11]. The shortage of supply of COVID-19 vaccines in many developing countries has resulted in the implementation of different immunization strategies, with varying results. For example, in Kenya, prioritization is given to the adult population [2].A mathematical model provides insight into infectious disease transmission and control [12, 13]. Many mathematical models on the dynamics of COVID-19 with control strategies have been developed and studied, for example, in studies carried out in [9, 10]. Optimal control is regarded as a powerful mathematical tool for optimizing control problems in various fields [14–18]. In the context of vaccination, the dynamics of disease transmission and immunity in a population are typically modelled by using network compartmental approaches [19]. A mathematical model for the dynamics of COVID-19 is presented in [18]. An optimal control function is added to the model in order to effectively control the outbreak. The main controls are isolation, quarantine, and hospitalization. The result shows that adopting the available control measures to their full potential will greatly reduce infectious populations.A new COVID-19 model with an optimal control analysis is presented in [20]. The authors considered four different controls to develop an optimal control model, including prevention, vaccination control, rapid screening of people who are exposed, and people who are infected without screening. The forward-backward Runge–Kutta method is used to resolve the model with and without control. The findings suggest that control can be useful for reducing infected individuals and improving population health. A novel fractional-order mathematical model that explains the behaviour of COVID-19 in Ethiopia is developed and analyzed in [21]. An inexact Newton iterative method is used to solve the model system. The model considers the impact of various control techniques on disease transmission.Zamir et al. in [22] developed a mathematical model to study the transmission dynamics of COVID-19. The study used a nonclinical approach to investigate optimal control. A COVID-19 mathematical model that takes into account the susceptible, exposed, infected, asymptotic, quarantine/isolation, and recovered classes is developed and analyzed in [23]. Elasticity and sensitivity analysis shows that the model is more sensitive to transmission rates from infected to exposed classes than transmission rates from susceptible to exposed classes.Alqarni et al. [24] developed a new mathematical model for the transmission dynamics of the coronavirus (COVID-19) using cases reported in the Kingdom of Saudi Arabia from March 2, 2020, to July 31, 2020. The model’s stability results are shown when the basic reproduction number is R0<1, and the model is locally asymptotically stable. Using the PRCC method, they show some important parameters that are more sensitive to the basic reproduction number R0. The sensitive parameters that act as control parameters and can reduce and control infection in the population are depicted graphically. The study results suggested that if controls are followed, the infection rate in the Kingdom of Saudi Arabia can be reduced significantly. In this study, an extended SEIR COVID-19 dynamics model with vaccination is proposed to investigate the impact of vaccination as a control intervention on the pandemic spread.


## 2. Model Formulation

A susceptible-vaccinated-exposed-infected-recovered (SVEIR) model with a standard incidence rate is developed and analyzed. Despite the fact that there are underlying conditions that predispose someone to a higher risk of contracting COVID-19 disease, the population under study is assumed to have an equal level of susceptibility. Individuals in this category are known as susceptible individuals denoted bySt. Several COVID-19 vaccine candidates have been developed and are already in use [10]. Individuals who have been vaccinated are hereby categorised as the vaccinated class, denoted by Vt. When susceptible individuals are exposed to the virus, they undergo an incubation period, which is the time between virus exposure and the onset of symptoms. This is on average 5–6 days, but it can be as long as 14 days [25]. The incubating individuals in this study are categorised as the exposed class denoted by Et. Following the end of the exposure period, the incubating individual’s transit into the infection class is denoted by It. The infectious potential of COVID-19 is significantly greater just before or within the first five days of symptom onset [26]. The infected class It, in this study, includes both presymptomatic and asymptomatic cases. Depending on the severity of symptoms and the intervention strategies in place, the infected individuals may die or recover from the disease. Thus, the class Rt denotes the number of individuals who have recovered from the infection. The study assumes that infection with COVID-19 confers temporary immunity upon recovery.The recruitment of individuals into the susceptible class is through births at the per capita rateΛ. Following an outbreak of the disease, several containment measures are implemented to control the spread. The susceptible humans acquire infection from an infectious individual via a force of infection βI, where β is the transmission coefficient. The saturated incidence rates given by βSI/1+ηI and β1VI/1+ηI tend to the saturation level when I grows large. The constant η measures the saturation effects caused by infective individuals in the presence of containment measures. The proportion of vaccination is given by γ, 0<γ≤1. The transmission coefficient of the vaccinees is denoted by β1:=1−ϑβ, where the parameter ϑ, 0<ϑ≤1 measures the vaccine efficacy. The vaccinees are considered to have acquired partial immunity against COVID-19 subject to the type of vaccine (efficacy level) administered to an individual. The transmission coefficient of the vaccinees, β1<β since the vaccinees are assumed to have acquired a vaccine-induced immunity [27, 28]. The rate at which exposed individuals become infectious is denoted by ε. The rates of natural and disease-induced mortality are μ and δ, respectively. The recovery rate from the infection is taken as λ.The above model description translates into the following schematic flow diagram in Figure1.Figure 1
A schematic flow diagram of the COVID-19 dynamics model with vaccination.From the description above, the corresponding mathematical model is represented by the following set of ordinary differential equations:(1)dSdt=Λ−βSI1+ηI−μ+γS,dVdt=γS−β1VI1+ηI−μV,dEdt=βS+β1VI1+ηI−μ+εE,dIdt=εE−μ+λ+δI,dRdt=λI−μR,where β1<β, S=St,V=Vt,E=Et,I=It,R=Rt, S,V,E,I,R∈ℝ+5.


## 3. Well-Posedness of the Model

In this section, well-posedness of the model solutions is discussed. Model (1) describes the human population, and therefore, its solutions, as shown below, are positive and bounded for all time t≥0.


### 3.1. Boundedness of Solutions

Model (1) is analyzed in a suitable feasible region(2)Ω=S,V,E,I,Rt∈ℝ+5:St+Vt+Et+It+Rt≤Λμ,where St+Vt+Et+It+Rt=Nt gives the total human population.Using Proposition1 below, the model solutions are shown to be bounded for all t≥0 in the region Ω.Proposition 1.
For all timet≥0, the solutions of model (1) are invariant in the region Ω.Proof.
The solutions of model (1) are positively invariant of Ω, i.e., all solutions start in Ω and remain in the region Ω for all t≥0. The rate of change of human population Nt is given by(3)dNdt=dSdt+dVdt+dEdt+dIdt+dRdt.
which implies that(4)dNdt≤Λ−μN.
By variation-of-constant formula, it follows that(5)limsupt⟶+∞Nt≤Λμ.
Thus,Nt≤Λ/μ. This implies that the solution set St,Vt,Et,It,Rt is bounded in the feasible region Ω, i.e.,(6)Ω=St,Vt,Et,It,Rt|St+Vt+Et+It+Rt≤Λμ,S,V,E,I,R0≥0.
Hence, all solutions of model (1) are bounded in the region Ω.


### 3.2. Positivity of Solutions

Proposition 2.
All solution setsSt,Vt,Et,It,Rt of model (1) with non-negative initial conditions are positive ∀t>0.Proof.
Considering the first equation in model (1), that is,(7)dSdt=Λ−βSI1+ηI−μ+γS,where(8)dSdt≥−βI1+ηI+μ+γS.
integration by variable separation to get(9)∫dSS≥∫−βI1+ηI+μ+γdt.
thus(10)St≥S0exp−∫0tKdτ+μ+γt,where K=βI/1+ηI. This implies that(11)St≥0∀t≥0.
In a similar way, all the other variables can be shown to be positive∀t≥0. Hence, all solutions of model (1) are positive in the region Ω.
Clearly, from Propositions1 and 2, all solutions of model (1) are shown to be positively invariant in the region Ω. Thus, model (1) is mathematically and epidemiologically well posed in a biological feasible region Ω.


## 3.1. Boundedness of Solutions

Model (1) is analyzed in a suitable feasible region(2)Ω=S,V,E,I,Rt∈ℝ+5:St+Vt+Et+It+Rt≤Λμ,where St+Vt+Et+It+Rt=Nt gives the total human population.Using Proposition1 below, the model solutions are shown to be bounded for all t≥0 in the region Ω.Proposition 1.
For all timet≥0, the solutions of model (1) are invariant in the region Ω.Proof.
The solutions of model (1) are positively invariant of Ω, i.e., all solutions start in Ω and remain in the region Ω for all t≥0. The rate of change of human population Nt is given by(3)dNdt=dSdt+dVdt+dEdt+dIdt+dRdt.
which implies that(4)dNdt≤Λ−μN.
By variation-of-constant formula, it follows that(5)limsupt⟶+∞Nt≤Λμ.
Thus,Nt≤Λ/μ. This implies that the solution set St,Vt,Et,It,Rt is bounded in the feasible region Ω, i.e.,(6)Ω=St,Vt,Et,It,Rt|St+Vt+Et+It+Rt≤Λμ,S,V,E,I,R0≥0.
Hence, all solutions of model (1) are bounded in the region Ω.


## 3.2. Positivity of Solutions

Proposition 2.
All solution setsSt,Vt,Et,It,Rt of model (1) with non-negative initial conditions are positive ∀t>0.Proof.
Considering the first equation in model (1), that is,(7)dSdt=Λ−βSI1+ηI−μ+γS,where(8)dSdt≥−βI1+ηI+μ+γS.
integration by variable separation to get(9)∫dSS≥∫−βI1+ηI+μ+γdt.
thus(10)St≥S0exp−∫0tKdτ+μ+γt,where K=βI/1+ηI. This implies that(11)St≥0∀t≥0.
In a similar way, all the other variables can be shown to be positive∀t≥0. Hence, all solutions of model (1) are positive in the region Ω.
Clearly, from Propositions1 and 2, all solutions of model (1) are shown to be positively invariant in the region Ω. Thus, model (1) is mathematically and epidemiologically well posed in a biological feasible region Ω.


## 4. The Vaccine Reproduction Number and the Disease-Free Equilibrium

In this section, the disease-free equilibrium and the vaccine reproduction numberRV of the model (1) are computed.


### 4.1. The Disease-Free Equilibrium

The disease-free equilibrium (DFE) point of model (1) is defined as the state in which there is no COVID-19 infection in the population under study.Proposition 3.
There exists a DFE of model (1) given by(12)E0=S0,V0,E0,I0,R0=Λμ+γ,γΛμμ+γ,0,0,0.Proof.
At the disease-free equilibrium point, there is no infection in the population. Therefore, withE=I=R=0, then model (1) is given as follows:(13)Λ−βSI1+ηI−μ+γS=0,γS−β1VI1+ηI−μV=0,βS+β1VI1+ηI−μ+εE=0,εE−μ+λ+δI=0,λI−μR=0.
Solving model (13) S0=Λ/μ+γ and V0=γΛ/μμ+γ. Therefore, model (1) has a disease-free equilibrium given by(14)E0=Λμ+γ,γΛμμ+γ,0,0,0.


### 4.2. Reproduction Number

The basic reproduction number, usually denoted byR0, is defined as the average number of secondary infections due to a single infectious individual introduced into a fully susceptible population during his/her period of infectivity [29]. The basic reproduction number is the spectral radius of the matrix, FV−1, R0=ρFV−1, where F and V are the next-generation matrices [30]. The operator FV−1, the next generation matrix, is formed from matrices of partial derivatives of ℱi (the rate of appearance of new infection in the ith compartment) and Vi=Vi−−Vi+ (the rate of transfer/transition rate in and out of the disease compartment i) with respect to the infected compartments (E and I) evaluated at DFE. The matrices F and V are given by(15)F=∂ℱiE0∂xj,V=∂ViE0∂xj,where the transition matrices F and V evaluated at E0 are given by(16)F=0Λβμ+β1Λγμμ+γ00,and(17)V=μ+ε0−εμ+λ+δ.MatrixV is invertible and(18)V−1=1μ+ε0εμ+εμ+λ+δ1μ+λ+δ.Thus, the matrixFV−1 is(19)FV−1=Λεβμ+1−ϑγμμ+γμ+εμ+λ+δΛβμ+1−ϑγμμ+γμ+εμ+λ+δ00.Since model (1) involves vaccination as an intervention, its associated reproduction number is called the vaccine reproduction number denoted by RV. This is the threshold quantity that can predict the spread of the disease in a given population in the presence of vaccination. The vaccine reproduction number computed by using the next-generation matrix approach above is given by(20)RV=Λεβμ+1−ϑγμμ+εμ+λ+δμ+γ,RV=R0μ+1−ϑγμ+γ.If the vaccine efficacyϑ=0, then Rv=R0 which is the basic reproduction number given by R0=Λεβ/μμ+εμ+λ+δ.The vaccine reproduction number,RV, is the measure of the severity of an epidemic in the presence of vaccination and one of the most important parameters since it determines whether or not the disease will invade a population. Epidemiologically, if RV<1, then by definition, the infection does not spread in the population. On the other hand, if RV>1, then the infection spreads in the population and may result into an epidemic.


## 4.1. The Disease-Free Equilibrium

The disease-free equilibrium (DFE) point of model (1) is defined as the state in which there is no COVID-19 infection in the population under study.Proposition 3.
There exists a DFE of model (1) given by(12)E0=S0,V0,E0,I0,R0=Λμ+γ,γΛμμ+γ,0,0,0.Proof.
At the disease-free equilibrium point, there is no infection in the population. Therefore, withE=I=R=0, then model (1) is given as follows:(13)Λ−βSI1+ηI−μ+γS=0,γS−β1VI1+ηI−μV=0,βS+β1VI1+ηI−μ+εE=0,εE−μ+λ+δI=0,λI−μR=0.
Solving model (13) S0=Λ/μ+γ and V0=γΛ/μμ+γ. Therefore, model (1) has a disease-free equilibrium given by(14)E0=Λμ+γ,γΛμμ+γ,0,0,0.


## 4.2. Reproduction Number

The basic reproduction number, usually denoted byR0, is defined as the average number of secondary infections due to a single infectious individual introduced into a fully susceptible population during his/her period of infectivity [29]. The basic reproduction number is the spectral radius of the matrix, FV−1, R0=ρFV−1, where F and V are the next-generation matrices [30]. The operator FV−1, the next generation matrix, is formed from matrices of partial derivatives of ℱi (the rate of appearance of new infection in the ith compartment) and Vi=Vi−−Vi+ (the rate of transfer/transition rate in and out of the disease compartment i) with respect to the infected compartments (E and I) evaluated at DFE. The matrices F and V are given by(15)F=∂ℱiE0∂xj,V=∂ViE0∂xj,where the transition matrices F and V evaluated at E0 are given by(16)F=0Λβμ+β1Λγμμ+γ00,and(17)V=μ+ε0−εμ+λ+δ.MatrixV is invertible and(18)V−1=1μ+ε0εμ+εμ+λ+δ1μ+λ+δ.Thus, the matrixFV−1 is(19)FV−1=Λεβμ+1−ϑγμμ+γμ+εμ+λ+δΛβμ+1−ϑγμμ+γμ+εμ+λ+δ00.Since model (1) involves vaccination as an intervention, its associated reproduction number is called the vaccine reproduction number denoted by RV. This is the threshold quantity that can predict the spread of the disease in a given population in the presence of vaccination. The vaccine reproduction number computed by using the next-generation matrix approach above is given by(20)RV=Λεβμ+1−ϑγμμ+εμ+λ+δμ+γ,RV=R0μ+1−ϑγμ+γ.If the vaccine efficacyϑ=0, then Rv=R0 which is the basic reproduction number given by R0=Λεβ/μμ+εμ+λ+δ.The vaccine reproduction number,RV, is the measure of the severity of an epidemic in the presence of vaccination and one of the most important parameters since it determines whether or not the disease will invade a population. Epidemiologically, if RV<1, then by definition, the infection does not spread in the population. On the other hand, if RV>1, then the infection spreads in the population and may result into an epidemic.


## 5. Sensitivity Analysis

The degree to which an input parameter influences a model’s output is known as parameter sensitivity. Sensitivity analysis of the basic reproductive number can be used to develop a mitigation strategy that will slow the spread of the pandemic by loweringRV [14]. Sensitive parameters are those that have a significant impact on the transmission dynamics of an infection. Using the normalized forward sensitivity indices [31], the sensitivity indices are computed by(21)ϒMRV=∂RV∂M×MRV,where M is the parameter whose sensitivity index is computed.Table1 gives a summary of the sensitivity indices of RV evaluated at the baseline parameters values given in Table 2.Table 1
Sensitivity indices ofRV with respect to the model parameters.
ParameterDescriptionSensitivity indexΛRecruitment rate+1βTransmission coefficient+1εTransition rate fromE to I−0.999843μNatural death rate+7.81816576×10−5γRate of vaccination−0.8984463ϑVaccine efficacy−0.9998436λRecovery rate−0.999679δDisease mortality rate−0.000008237Table 2
The descriptive summary of the model parameters.
ParameterDescriptionUnit/value unitsSourceΛRecruitment rate3.178×10−5 day−1[28]μNatural death rate3.91×10−5 day−1[10]δDisease mortality rate0.103×10−5 day−1[28]βTransmission coefficient0.02 day−1Estimatedβ1Transition rate from V to E0.05 day−1EstimatedϑVaccine efficacy(0–1.0)VariableγRate of vaccination(0–1.0)VariableεTransition rate fromE to I0.5 day−1EstimatedηHuman saturation constant0.05EstimatedλHuman recovery rate0.125 day−1[32]From the sensitivity analysis as presented in Table1, an increase of the rate of vaccination γ by 1% leads to a decrease of the effective reproduction RV by 0.89844626%. An increase of the vaccine efficacy ϑ by 1% leads to a decrease in the effective reproductive number RV by 0.9990436%. Clearly, the rate of vaccination and the vaccine efficacy are among the sensitive parameters of RV. An increase in the rate of vaccination, with high vaccine efficacy, leads to a decrease in the vaccine reproduction number. Consequently, control strategies should target an increased rate of vaccination and administration of vaccines of high efficacy levels.


## 6. Optimal Analysis with Vaccination Program as a Control Intervention

An optimal control problem is constructed to control the spread of the COVID-19 virus and optimize the vaccination program. Therefore, in this section, optimal control is performed to understand the effects of the sensitive parameters on the optimum vaccination program. To determine the optimum vaccination program for COVID-19, let the variablesγ∗ and β1∗=:1−ϑβ be the control variables. As a result, an optimal control problem is constructed, with the goal of reducing the number of individuals infected with COVID-19. The following objective function is created to accomplish this:(22)J=∫0τP0S+P1V+P2E+P3γ2+P4β12dt.where 0,τ is the entire time horizon over the control applied and P0, P1, P2, P3, P4 are positive weights that balance the relative importance of terms in the objective functional J. An optimal control γ∗, β1∗ is chosen as(23)Jγ∗,β1∗=minJγ,β1.such thatγ,β1 are measurable with 0≤γ≤β1≤1. This is the necessary condition that the optimal control must satisfy from Pontryagin’s maximum principle [33, 34]. The Hamiltonian function is(24)H=P0S+P1V+P2E+P3γ2+P4β12+ΦsΛ−βSI1+ηI−μ+γS+ΦVγS−β1VI1+ηI−μV+ΦEβS+β1VI1+ηI−μ+εE+ΦIεE−μ+λ+δI+ΦRλI−μR,where ΦS, ΦV, ΦE, ΦI, and ΦR are the adjoint variables. To obtain the expression of optimal controls in order to minimize the number of infections in the population and the cost of control strategies, the following proposition is applied.Proposition 4.
For the optimal controlγ,β1 that minimizes Jγ,β1, the adjoint variables ΦS, ΦV, ΦE, ΦI, and ΦR satisfy the following ordinary differential equations:(25)dΦSdt=P0−ΦSβI1+ηI−μ+γΦS+γΦV+ΦEβI1+ηI,dΦVdt=P1−ΦVβI1+ηI−μΦV+ΦEβ1I1+ηI,dΦEdt=P2−μ+εΦE+εΦI,dΦIdt=−ΦSβS1+ηI2−ΦVβ1V1+ηI2+ΦEβS+β1V1+ηI2−ΦIμ+λ+δ+λΦR,dΦRdt=μΦR.with transversality conditions(26)ΦSτ=ΦVτ=ΦEτ=ΦIτ=ΦRτ=0.Thus, the optimal control takes the characterization form, which is given as follows:(27)γ∗=max0,min1,SΦS−ΦV2P3,(28)β1∗=max0,min1,VIΦV−ΦE21+ηIP4.Proof.
The HamiltonianH in equation (25) is differentiated with respect to the state variables, S,V,E,I,R, respectively. Thus, the adjoint of the system can be written as(29)ΦM=−∂H∂M,for(30)M=S,V,E,I,R.
By Pontryagin’s maximum principle [34], H can be maximized with respect to γ and β1, that is,(31)0=∂H∂γ|γ∗=2P3γ−ΦSS−ΦVS,⇒γ∗=SΦS−ΦV2P3.0=∂H∂β1|β1∗=2P4β1−ΦVVI1+ηI+ΦEVI1+ηI,⇒β1∗=VIΦV−ΦE21+ηIP4.
Taking the bounds onγ and β1 into account, the characterization of γ∗ and β1∗ is obtained as shown in (27) and (28), respectively. Now, using the control arguments 0≤γ≤β1≤1, then we obtain(32)γ∗=0ifξ1∗≤0ξ1∗if0<ξ1∗<11ifξ1∗≥1,β1∗=0ifξ2∗≤0ξ2∗if0<ξ2∗<11ifξ2∗≥1,where(33)ξ1∗=SΦS−ΦV2P3,ξ2∗=VIΦV−ΦE21+ηIP4.
Since the optimal control switches at most once, then the control objective function constructed in this study is at optimum.
In recent decades, control theory has been widely applied in many fields. Optimal control, particularly in epidemiology, could be very useful in controlling mathematical models depicting the spread of infectious diseases [35]. The appropriate regulation of disease dynamics is specified in the form of restrictions according to the biological interpretation of the objective function. The control’s objectives must be met exclusively within these limits. Furthermore, lowering illness prevalence requires reducing the total number of infectious cases. This entails increasing the rate of vaccination γ while decreasing the coefficient of contracting the infection β1=1−ϑβ. Increased vaccination rates and the administering of highly efficacious vaccines would help maximize control strategies against COVID-19 transmission dynamics.


## 7. Numerical Simulation and Discussion

The parameters used for simulation are obtained from literature, and others are estimates, as listed in Table2. These parameter values are varied within realistic limits.For purposes of simulation, unless otherwise stated, the initial populations are taken to beS0=3000, V0=2000, E0=1500, I0=1200, and R0=1000.The numerical simulation aims to analyze the change in the state of COVID-19 virus progression over time, as well as the impact of variation in the vaccination rate and vaccine efficacy on COVID-19 transmission dynamics. This is accomplished by varying the parametersϑ and γ while holding the other parameters constant. Simulation analyses of the model (1) are presented in Figures 2 and 3.Figure 2
The effect of variation in vaccine efficacyϑ on infection dynamics.Figure 3
The effect of variation in the vaccination rateγ on infection dynamics.Figure2 shows the effect of variation in vaccine efficacy ϑ on infection dynamics. When the vaccine administered is of high efficacy, say 90% efficacy, it would take less time for the pool of infected individuals to reduce. On the other hand, when the vaccine efficacy is low, it would take a longer time for the number of infected individuals to reduce, and therefore, COVID-19 infection would persist in the population. Figure 3 shows the effect of variation in the vaccination rate γ on infection dynamics. It is shown that when the rate of vaccination is high, the number of infected individuals decreases sharply and vice versa. Figure 4 shows the graph of the vaccine reproduction number RV against the vaccination rate γ. It is observed that with an increase in the rate of vaccination γ, the vaccine reproduction number decreases sharply. Figure 5 shows the 3D plot of the vaccine reproduction number RV against the rate of vaccination γ and vaccine efficacy ϑ. It is observed that with an increase in both the rate of vaccination and vaccine efficacy, the vaccine reproduction number decreases sharply.Figure 4
Plot of the vaccine reproduction numberRV against the rate of vaccination γ.Figure 5
Plot of the vaccine reproduction numberRV against the rate of vaccination γ and vaccine efficacy ϑ.Vaccination offers a very powerful method of COVID-19 disease control. The critical level of vaccination of above50% with an>80% vaccine efficacy level would be required to reduce the severity and eventually eradicate the infection from the population. Herd immunity against COVID-19 can be achieved by immunizing a significant proportion of the susceptible population. To prevent a large outbreak in the future, in addition to vaccination, a variety of other control measures (e.g., COVID-19 protocols such as the use of masks, social distancing, proper hygiene, and so on) are required [14, 36].


## 8. Conclusion

In this paper, the effect of vaccination as an intervention is investigated through an SVEIR model. The well-posedness of the model and the existence of disease-free equilibrium are shown. Sensitivity analysis here suggests that the rate of vaccination,γ, and the vaccine efficacy ϑ are among the most influential parameters of RV. For optimal control, increasing the vaccination rate and administering highly efficacious vaccines would help maximize control strategies against COVID-19 transmission dynamics. This is illustrated by the numerical simulations. This study, like others [19, 37, 38], confirms that vaccination may possibly eradicate the infection in the population. Future research could focus on an age-dependent SVEIR model to determine the suitable vaccination strategy to be applied with respect to age.



---

*Source: 1476607-2022-10-15.xml*