Next Article in Journal
Analysis of the Mathematical Modelling of COVID-19 by Using Mild Solution with Delay Caputo Operator
Next Article in Special Issue
A Secure and Lightweight Multi-Party Private Intersection-Sum Scheme over a Symmetric Cryptosystem
Previous Article in Journal
Enumerating Subtrees of Flower and Sunflower Networks
Previous Article in Special Issue
Edge Computing-Based VANETs’ Anonymous Message Authentication
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Extinction and Ergodic Stationary Distribution of COVID-19 Epidemic Model with Vaccination Effects

1
School of Mathematics, Liaocheng University, Liaocheng 252000, China
2
School of Mathematical Sciences, Suzhou University of Science and Technology, Suzhou 215000, China
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(2), 285; https://doi.org/10.3390/sym15020285
Submission received: 8 October 2022 / Revised: 9 November 2022 / Accepted: 24 December 2022 / Published: 19 January 2023

Abstract

:
Human society always wants a safe environment from pollution and infectious diseases, such as COVID-19, etc. To control COVID-19, we have started the big effort for the discovery of a vaccination of COVID-19. Several biological problems have the aspects of symmetry, and this theory has many applications in explaining the dynamics of biological models. In this research article, we developed the stochastic COVID-19 mathematical model, along with the inclusion of a vaccination term, and studied the dynamics of the disease through the theory of symmetric dynamics and ergodic stationary distribution. The basic reproduction number is evaluated using the equilibrium points of the proposed model. For well-posedness, we also test the given problem for the existence and uniqueness of a non-negative solution. The necessary conditions for eradicating the disease are also analyzed along with the stationary distribution of the proposed model. For the verification of the obtained result, simulations of the model are performed.

1. Introduction

For the control and extinction of COVID-19, several effective measurements have been taken from many researchers and policymakers of each country of the globe. Many of them are very economical and have produced a disturbance in the routine activities of social media and psychological effects. Controlling of the said disease viruses for all people showing signs of infection or not is the significant and key work of the higher-ups and managerial staff. Still now in the process of developing safety against COVID-19, the most important tool is vaccination [1,2]. This infection may spread through human-to-human contact in the society, like (MERS-CoV), which is transmitted from civet cats into the human environment. Various efforts have been made against it, which can be studied in [3,4].
The mathematical formulation may well present the dynamical behavior of the infectious diseases of different epidemics. Therefore, different approaches of mathematical modeling for the biological dynamics have been used in the past history. Among them is deterministic mathematical modeling, which mostly discusses the idealistic situation of real-world problems. It depends on input and output data. The stochastic modeling is another approach which mostly describes the real situation. It presents both the randomness of various compositions with division and also the deterministic version. Because of the noise terms as an input source, the said division may have dynamics of uncertainty [5,6,7,8].
So, for vaccination, the conversion of real global problems into mathematical terms is a very significant tool for the analysis of epidemics. Epidemiological infections are transferred through mathematical formulation to various types of differential equations and systems of it. Policy makers and many other scholars also provide some ideas and parameters which can be used for reducing infections [9]. From the beginning of the twenty-first century, the conversion of real-world problems to mathematical terms made the information very easy and made the situation very simple regarding future predictions [10,11]. Mostly, the infections are written in mathematical formulation. Different mathematical formulation are applied to mostly infectious diseases. The techniques of mathematical modeling are divided into different mechanisms. Among them are the deterministic, stochastic and difference equations. The stochastic-type approach is much more realistic as compared to other mechanisms of mathematical formulations. Such a scheme type gives the result in the form of up and down for unknown quantities along with the distribution of an approximate solution [10,11,12,13]. For such dynamics, one can study different articles, such as that given in reference [14,15].
The stochastic analysis of real-world phenomena is very good, as it involves both the input and output of the data of that responsible for the transmission of the disease. Due to this, it answers all questions related to infectious disease modeling. It discusses how to construct a mathematical model for infectious disease. It also answers the question of how to formulate the intermediate frequently transmitted disease. The rate of disease and its suitable position in the model is also discussed by the stochastic approach. So for the need of these points, we formulate the COVID-19 problem by the mathematical term of stochastic approach. Changing the densities of different model quantities, we construct the model for the said infection with a long period of time investigation.
Several biological problems have the aspects of symmetry, and this theory has many applications in explaining the dynamics of biological models. The systematic theory of symmetry has several applications in the development and study of models of biological shapes and phenomena. In engineering and physics, symmetries are frequently accurate or nearly so, whereas in the biological sciences, perfect symmetry is uncommon. Because symmetries are only approximations in most biological systems, their application in models is an idealization. This type of idealization is beneficial because it simplifies the mathematical analysis and also because systems with approximated symmetries frequently reflect idealized symmetric models more precisely than ordinary asymmetric ones. The symmetrical methods are especially suited in biological environments where somewhat consistent patterns are found.
The analysis of symmetry features of ODEs and PDEs is a popular and well-established topic (see [16,17]) and provides a valuable tool for better comprehension of the solution’s qualitative behavior. In contrast, a similar theory of SDE symmetries was just recently constructed; the readers are suggested to see [18,19,20]. In this work, we studied a stochastic COVID-19 model by using a very general and elegant approach for solving a SDE via symmetries in a dynamical perspective. Unfortunately, the scale and complexity of many models in mathematical biology renders a brute force application of symmetry methods impractical. In other words, the analysis of SDEs using symmetry methods is non-standard in mathematical biology; thus, we will focus mainly on other available techniques for SDEs and will give less attention to the underlying symmetry methods.
The remaining part of this research paper is organized as follows: In Section 2, we develop the COVID-19 mathematical formulation with stochastic fluctuation in the rate of infection. In Section 3, we present the qualitative analysis of a non-negative solution for a long duration of time to the proposed stochastic model. We give some important conditions for the COVID-19 infection to vanish it from the society in Section 4. In Section 6, the condition for the stationary distribution existence is provided. The scheme of the analytical solution and their graphical representation are given in Section 7. For the interested reader, we conclude the article with some future remarks in Section 8.

2. Models Formulation

The virus of COVID-19 has a proper incubation period and will result in confinement in the beginning from the community and after the detection of disease. Healthy individuals will have an incubation duration before disease onset. We construct the COVID-19 model with five ordinary equations of derivative. The quantities of the model are composed of Susceptible ( S ) , Immunize or Vaccinated ( V ) , Asymptomatic or Exposed ( E ) , Symptomatic or Infectious ( I ) , hospitalized ( H ) , and recover individuals ( R ) , i.e., S ( t ) + V ( t ) + E ( t ) + I ( t ) + H ( t ) + R ( t ) ) = N ( t ) represents the whole density size, which is varying with time t. The finalized quantities epidemic model for the COVID-19 dynamical behavior with the vaccination effect is given as follows:
d S ( t ) d t = Π η S ( t ) I ( t ) N ( κ + δ ) S ( t ) , d V ( t ) d t = κ S ( t ) ( 1 τ ) η V ( t ) I ( t ) N δ V ( t ) , d E ( t ) d t = η S ( t ) I ( t ) N + ( 1 τ ) η V ( t ) I ( t ) N ( t ) ( σ + δ ) E ( t ) , d I ( t ) d t = σ E ( t ) ( ρ 1 + ρ 2 + ρ 3 + δ ) I ( t ) , d H ( t ) d t = ρ 1 I ( t ) ( μ 1 + μ 2 + δ ) H ( t ) , d R ( t ) d t = ρ 2 I ( t ) + μ 1 H ( t ) δ R ( t ) .
The used parameters in the proposed model are given in Table 1.
The vaccination parameter is assumed to be effective in the analyzed article, i.e., it has no cure rate of the COVID-19 virus. By this, powerful immunized density is required when coming into contact with infectious people. It should be noted that 0 < τ < 1 ( τ = 1 shows a big effective vaccine, while τ = 0 represents a vaccine that has no immunity).
The spreading dynamical analysis of the aforesaid COVID-19 problem is fully investigated by the basic reproductive value, which is given as
R 0 D = ρ F V 1 = α η Π [ ( 1 τ ) ρ + d ] d ( d + ρ ) ( α + d ) ( δ + d ) .
Further, if R 0 D < 1 , (DFE) X 0 = S 0 , V 0 , E 0 , I 0 , H 0 , R 0 = Π ρ + d , ρ Π d ( ρ + d ) , 0 , 0 , 0 , 0 is local asymptotic stable if R 0 D < 1 , and it is global asymptotic stable in the same region. Next, if R 0 D > 1 , then DFE is unstable. In such a situation, Equation (1) contains another equilibrium point called the endemic equilibrium (EE) X * = S * , V * , E * , I * , H * , H * which is local and global asymptotic stable if R 0 D < 1 and not stable if R 0 D > 1 .
It is noted that the uncertainty of the duration of incubation and the variation of detection along with population movements with stochastically given dynamics of the proposed system (1) incorporate random perturbation by transferring into a mathematical investigation. In this article, we take the random perturbation, which is directly related to the changing of Healthy ( S ) , Immunized or Vaccinated ( V ) , Asymptomatic or Exposed ( E ) , Symptomatic or Infected ( I ) , Hospitalized ( H ) , and Recover ( R ) under the effect of white noise as given in the stochastic model
d S ( t ) = Π η S ( t ) I ( t ) N ( κ + δ ) S ( t ) d t + Φ 1 S ( t ) d B 1 ( t ) , d V ( t ) = κ S ( t ) ( 1 τ ) η V ( t ) I ( t ) N δ V ( t ) d t + Φ 2 V ( t ) d B 2 ( t ) , d E ( t ) = η S ( t ) I ( t ) N + ( 1 τ ) η V ( t ) I ( t ) N ( σ + δ ) E ( t ) d t + Φ 3 E d B 3 ( t ) , d I ( t ) = σ E ( t ) ( ρ 1 + ρ 2 + ρ 3 + δ ) I ( t ) d t + Φ 4 I ( t ) d B 4 ( t ) , d H ( t ) = ρ 1 I ( t ) ( μ 1 + μ 2 + δ ) H ( t ) d t + Φ 5 H ( t ) d B 5 ( t ) , d R ( t ) = ρ 2 I ( t ) + μ 1 H ( t ) δ R ( t ) d t + Φ 6 R ( t ) d B 6 ( t ) .
here, B 1 ( t ) , B 2 ( t ) , B 3 ( t ) , B 4 ( t ) , B 5 ( t ) , B 6 ( t ) are the parameters for the Brownian motions, and Φ 1 , Φ 2 , Φ 3 , Φ 4 , Φ 5 , Φ 6 are the frequency intensity of the Gaussian white noise.

3. The Qualitative Analysis for the Positive Solution

In this section, we investigate the existence of a positive solution along with the uniqueness of the said solution of the proposed stochastic model (2).
Theorem 1.
We consider ( S t , V t , E t , I t , H t , R t ) as a positive solution of the given stochastic epidemic model (2) is one on t 0 with initial approximations ( S 0 , V 0 , E 0 , I 0 , R 0 ) R + 6 . Further, the root will always lie in R + 6 having one probability, as ( S 0 , V 0 , E 0 , I 0 , R 0 ) R + 6 t 0 mostly sure (a.s).
Proof. 
As the initial conditions, ( S 0 , V 0 , E 0 , I 0 , H 0 , R 0 ) R + 6 are coefficients that are pre-defined and local Lipschitz. Therefore, we have unique root ( S t , V t , E t , I t , H t , R t ) of the system over t [ 0 , τ e ) . For extra analysis of more time τ e , we give the reference [21,22]. For computation, the globally nature of the root, we have to derive that τ e = almost surely. Let us have a large positive value k 0 , all of the initial conditions on the state variable defined in [ 1 k 0 , k 0 ] . Consider that for all positive integer k k 0 , the time for vanishing is
τ k = t [ 0 , τ e ) : m i n { S t , V t , E t , I t , H t , R t } 1 k or m a x { S t , V t , E t , I t , H t , R t } .
In the whole article, the taking of inf ϕ = is chosen. Here, ϕ is the void set. From the output of τ k , we know that this increases as k goes to . Set τ = lim k and τ e τ , almost surely deriving τ = . The verification about the τ e = is noted and so, ( S t , V t , E t , I t , H t , R t ) will be in R + 5 almost surely t 0 . So, this will be sufficient to derive that τ e = almost surely. This shows that in the given case, there are two positive constant ϵ , from ( 0 , 1 ) and T:
P { T τ } > ϵ .
So, let k 1 k 0 and
P { T τ k } ϵ , k 1 k .
Next, we define a C 2 mapping H : R + 5 R + as follows:
G ( S , V , E , I , H , R ) = S + V + E + I + H + R 6 ( log S + log V + log E + log I + log H + log R ) .
It is validated that the H is a positive operator, and may be conformed from the statement 0 y logy 1 , 0 < y . Consider that k 0 K and 0 < T are any arbitrary fixed values. Using I t o ^ s for Equation (5) provides
d G ( S , V , E , I , H , R ) = L G ( S , V , E , I , H , R ) + Φ 1 ( S 1 ) d B 1 ( t ) + x i 2 ( V 1 ) d B 2 ( t ) + Φ 3 ( E 1 ) d B 3 ( t ) + Φ 4 ( I 1 ) d B 4 ( t ) + Φ 5 ( H 1 ) d B 5 ( t ) + Φ 6 ( R 1 ) d B 6 ( t ) .
In Equation (6), L G : R + 6 R + , whose definition is given by the given equation:
L G ( S , V , E , I , H , R ) = 1 1 S Π η S I N ( κ + δ ) S + Φ 1 2 2 + 1 1 V κ S ( 1 τ ) η V I N δ V ( t ) + Φ 2 2 2 + 1 1 E η S I N + ( 1 τ ) η V I N ( σ + δ ) E + Φ 3 2 2 + 1 1 I σ E ( ρ 1 + ρ 2 + ρ 3 + δ ) I + Φ 4 2 2 + 1 1 H ρ 1 I ( μ 1 + μ 2 + δ ) H + Φ 5 2 2 + 1 1 R ρ 2 I + μ 1 H δ R + Φ 6 2 2 .
L G Π + η + ( 1 τ ) η + 5 δ + κ + σ + ρ 1 + ρ 2 + ρ 3 + μ 1 + μ 2 + Φ 1 2 + Φ 2 2 + Φ 3 2 + Φ 4 2 + Φ 5 2 + Φ 6 2 2 : = K .
The remaining proof is almost the same as that presented in Theorem 2.1 of [22]. Hence, we must omit it here, and this completes the proof of the theorem. □

4. Extinction

This section is for the exploration of the parameters’ values for the extinction of infection given in model (2). We derive the significant result of our article in the form of a lemma as follows.
Lemma 1.
For any given initial value ( S 0 , V 0 , E 0 , I 0 , H 0 , R 0 ) R + 6 , the solution ( S t , V t , E t , I t , H t , R t ) of system (2) has the following properties:
lim t S ( t ) t = 0 , lim t V ( t ) t = 0 , lim t E ( t ) t = 0 , lim t I ( t ) t = 0 , lim t H ( t ) t = 0 , lim t R ( t ) t = 0 , a . s .
Furthermore, when δ > 1 2 ( Φ 1 2 Φ 2 2 Φ 3 2 Φ 4 2 Φ 5 2 Φ 6 2 ) holds, then
lim t 1 t 0 t S ( r ) d B 1 ( r ) = 0 , lim t 1 t 0 t V ( r ) d B 2 ( r ) = 0 , lim t 1 t 0 t E ( r ) d B 3 ( r ) = 0 , lim t 1 t 0 t I ( r ) d B 4 ( r ) = 0 , lim t 1 t 0 t H ( r ) d B 5 ( r ) = 0 , lim t 1 t 0 t R ( r ) d B 6 ( r ) = 0 , a . s .
Proof. 
The derivation of Lemma 1 is similar to [21], thus we skip it. □
Lemma 2
([11,12] (High number strong Law)). Consider that M = { M } t 0 is a continuous function of the real output, having a non-global Martingale and finishing at t = 0 , then
lim t M , M t = , a . s . , implies that lim t M t M , M t = 0 , a . s , and also lim t s u p M , M t t < 0 , a . s . , implies that lim t M t t = 0 , a . s .

5. Extinction

Define the following threshold quantity:
R 0 e = 2 σ η ( 2 τ ) ( σ + δ ) ρ 1 + ρ 2 + ρ 3 + δ + Φ 4 2 2 ( σ + δ ) 2 σ 2 Φ 3 2 2 .
Now, we construct a sufficient condition to ensure the extinction of COVID-19 in the stochastic formulation (2).
Theorem 2.
Let ( S t , V t , E t , I t , H t , R t ) be a unique and positive solution of the stochastic Formulation (2) with a positive starting data ( S 0 , V 0 , E 0 , I 0 , H 0 , R 0 ) . If the stochastic threshold R 0 e is strictly less than one, then E ( t ) and I ( t ) classes will go extinct almost surely, i.e.,
lim t E ( t ) = 0 , a n d lim t I ( t ) = 0 a . s .
Meanwhile,
lim t S ( t ) = Π ( κ + δ ) , a n d lim t V ( t ) = κ Π δ ( κ + δ ) a . s .
Proof. 
Firstly, we define the following combination linear of E and I :
ε ( E , I ) = a 1 E ( t ) + a 2 I ( t ) ,
where a 1 = σ and a 2 = σ + δ . By using Itô’s formula to ln ε , we obtain
d ln ε ( E , I ) = { σ η S I N ( σ E + ( σ + δ ) I ) + σ η ( 1 τ ) V I N ( σ E + ( σ + δ ) I ) ( ρ 1 + ρ 2 + ρ 3 + δ ) ( σ + δ ) I σ E + ( σ + δ ) I σ 2 Φ 3 2 E 2 + ( σ + δ ) Φ 4 2 I 2 2 ( σ E + ( σ + δ ) I ) 2 } d t + σ Φ 3 E σ E + ( σ + δ ) I d B 3 ( t ) = M 1 ( t ) + ( σ + δ ) Φ 4 I σ E + ( σ + δ ) I d B 4 ( t ) = M 2 ( t ) .
Then, we obtain that
d ln ε ( E , I ) σ η ( 2 τ ) ( σ + δ ) ( ρ 1 + ρ 2 + ρ 3 + δ + Φ 4 2 2 ) ( σ + δ ) 2 I 2 + σ 2 Φ 3 2 2 E 2 [ σ E + ( σ + δ ) I ] 2 d t + M 1 ( t ) + M 2 ( t ) , σ η ( 2 τ ) ( σ + δ ) ( ρ 1 + ρ 2 + ρ 3 + δ + Φ 4 2 2 ) ( σ + δ ) 2 σ 2 Φ 3 2 2 2 ( σ + δ ) 2 d t + M 1 ( t ) + M 2 ( t ) .
Integrating from 0 to t and then dividing by t on both sides leads to
ln ε ( E ( t ) , I ( t ) ) t σ η ( 2 τ ) ( σ + δ ) ( ρ 1 + ρ 2 + ρ 3 + δ + Φ 4 2 2 ) ( σ + δ ) 2 ( σ 2 Φ 3 2 2 ) 2 ( σ + δ ) 2 ln ε ( E ( 0 ) , I ( 0 ) ) t + σ Φ 3 t 0 t E ( r ) σ E ( r ) + ( σ + δ ) I ( r ) d B 3 ( r ) = M 1 i ( t ) + ( σ + δ ) Φ 4 t 0 t I ( r ) σ E ( r ) + ( σ + δ ) I ( r ) d B 4 ( r ) = M 2 i ( t ) ,
where M 1 i ( t ) and M 2 i ( t ) are two local continuous martingales. Then based on the strong law of large numbers, it implies that lim t t 1 M 1 i ( t ) = 0 and lim t t 1 M 2 i ( t ) = 0 . Provided that R 0 e < 1 , taking the superior limit of both sides leads to
lim t ln ε ( E ( t ) , I ( t ) ) t ( ρ 1 + ρ 2 + ρ 3 + δ + Φ 4 2 2 ) ( σ + δ ) 2 σ 2 Φ 3 2 2 2 ( σ + δ ) 2 R 0 e 1 < 0 .
Due to the positivity of the solution, we obtain
lim t E ( t ) = 0 , a n d lim t I ( t ) = 0 a . s .
Now, we will conclude the result presented in (12). From the first equation of the perturbed model (2), we obtain
S ( t ) S ( 0 ) t = Π η t 0 t S ( r ) I ( r ) N ( r ) d r ( κ + μ ) t 0 t S ( r ) d r + Φ 1 t 0 t S ( r ) d B 1 ( r ) .
Under the hypothesis R 0 e < 1 , we have
lim t S ( t ) = Π ( κ + δ ) a . s .
In the same way, we obtain from the second equation of (2) that
lim t V ( t ) = κ Π δ ( κ + δ ) a . s .
The derivation of the resulting Theorem 2 is completed. □

6. The Stationary Distribution of the Disease

In this section, we assume that the stochastically considered system has no endemic equilibrium point. So, the investigation of stability cannot be used as a scheme for studying epidemic infection. One may say that either the root is a Lie group or not, and if it is, then it must be so with the help of the theory of stationary distribution as the key task, and we may search for the Lie group of the endemic.With this task, we give the following along with the citation of a well-known theorem from Hasminskii [23]. Let
X ( t ) = 1 t 0 t X ( r ) d r .

Stationary Distribution

Let X ( t ) be defined on its domain and obey the Markov procedure (homogeneity of time) in R + n with different dynamics as follows:
d X ( t ) = b ( X ) d t + r k σ r d B r ( t ) .
The defuse matrix is
A ( X ) = [ a i j ( x ) ] , a i j ( x ) = r = 1 k σ r i ( x ) σ j r ( x ) .
Lemma 3
([11,12]). Technique X ( t ) has one distribution for its stationary condition m ( . ) if their lie group, the input data, has bounds with continuous boundary. U , U ¯ R d U ¯ closure U ¯ R d , have the following properties:
1. 
In both sides, open input U and in its neighbor, the smallest eigenvalue of A(t) has bounds that are separate.
2. 
If x R d U , the average time τ (at which a curve starts from x going to the set U) is of finiteness, and S u p x k E x τ < for every compact subset K R n . Next, if f ( . ) is an integrating function having measurement π, then
P lim T 1 T 0 T f ( X x ( t ) ) d t = R d f ( x ) π ( d x ) = 1
for all x R d .
Define a parameter
R 0 s = δ η σ η + κ + δ + Φ 1 2 2 σ + δ + Φ 3 2 2 ρ 1 + ρ 2 + ρ 3 + δ + Φ 4 2 2 .
Theorem 3.
The root ( S ( t ) , V ( t ) , E ( t ) , I ( t ) , H ( t ) , R ( t ) ) of system (2) is ergodic and has one distribution for its stationary point π ( . ) whenever R 0 S > 1 .
Proof. 
For the validation of Equation (2) of Lemma (3), we take a positive C 2 -operator
V 1 = S + V + E + I + H + R c 1 ln S c 2 ln E c 3 ln I ,
where c 1 , c 2 and c 3 are the positive constants to be calculated later in the rest of the sections. Using the Itô’s result and the take problem (2), we obtain
L ( S + V + E + I + H + R ) = Π δ ( S ( t ) + V ( t ) + E ( t ) + I ( t ) + H ( t ) + R ( t ) ) ρ 3 I μ 2 H , L ( ln S ) = Π S + η I N + ( κ + δ ) + Φ 1 2 2 , L ( ln V ) = κ S V + ( 1 τ ) η I N + δ + Φ 2 2 2 , L ( ln E ) = η S I N E ( 1 τ ) η V I N E + ( σ + δ ) + Φ 3 2 2 , L ( ln I ) = σ E I + ( ρ 1 + ρ 2 + ρ 3 + δ ) + Φ 4 2 2 , L ( ln H ) = ρ 1 I H + μ 1 + μ 2 + δ + Φ 5 2 2 , L ( ln R ) = ρ 2 I R μ 1 H R + δ + Φ 6 2 2 .
So, we have
L V 1 = Π δ ( S ( t ) + V ( t ) + E ( t ) + I ( t ) + H ( t ) + R ( t ) ) c 1 A S + c 1 η I N + c 1 κ + δ + Φ 1 2 2 c 2 η S I N E c 2 ( 1 τ ) η V I N E + c 2 σ + δ + Φ 3 2 2 c 3 σ E I + c 3 ρ 1 + ρ 2 + ρ 3 + δ + Φ 4 2 2 .
This implies that
L V 1 4 δ ( S ( t ) + V ( t ) + E ( t ) + I ( t ) + H ( t ) + R ( t ) ) × c 1 Π S × c 2 η S I N E × σ E I 1 4 + c 1 ( η + κ + δ + Φ 1 2 2 ) + c 2 σ + δ + Φ 3 2 2 + c 3 ρ 1 + ρ 2 + ρ 3 + δ + Φ 4 2 2 + Π .
Let
c 1 ( η + κ + δ + Φ 1 2 2 ) = c 2 σ + δ + Φ 3 2 2 = c 3 ρ 1 + ρ 2 + ρ 3 + δ + Φ 4 2 2 = Π .
Namely,
c 1 = Π ( η + κ + δ + Φ 1 2 2 ) , c 2 = Π σ + δ + Φ 3 2 2 , c 3 = Π ρ 1 + ρ 2 + ρ 3 + δ + Φ 4 2 2 .
Consequently,
L V 1 4 Π 4 δ η σ ( η + κ + δ + Φ 1 2 2 ) σ + δ + Φ 3 2 2 ρ 1 + ρ 2 + ρ 3 + δ + Φ 4 2 2 1 4 Π c 2 ( 1 τ ) η V I N E .
L V 1 4 Π ( R 0 S ) 1 / 4 1
In addition, we obtain
V 2 = c 4 ( S + V + E + I + H + R c 1 ln S c 2 ln E c 3 ln I ) ln S ln V ln H ln R + S ( t ) + V ( t ) + E ( t ) + I ( t ) + H ( t ) + R ( t ) , = ( c 4 + 1 ) ( S + V + E + I + H + R ) ( c 1 c 4 + 1 ) ln S c 2 c 4 ln E c 3 c 4 ln I ln V ln H ln R ,
here, the fixed c 4 > 0 will be computed later on. It very difficult to show that
lim inf ( S , V , E , I , H , R ) R + 6 U k V 2 ( S , V , E , I , H , R ) = + , as k ,
here U k = ( 1 k , k ) × ( 1 k , k ) × ( 1 k , k ) . Further for the new steps we show that V 2 ( S , V , E , I , H , R ) has one and at least one value V 2 ( S 0 , V 0 , E 0 , I 0 , H 0 , R 0 ) .
Partial order differentiation of V 2 ( S , V , E , I , H , R ) with respect to t ( S , V , E , I , H , R ) is as follows
V 2 ( S , V , E , I , H , R ) S = 1 + c 4 1 + c 1 c 4 S , V 2 ( S , V , E , I , H , R ) V = 1 + c 4 1 V , V 2 ( S , V , E , I , H , R ) E = 1 + c 4 c 2 c 4 E , V 2 ( S , V , E , I , H , R ) I = 1 + c 4 c 3 c 4 I , V 2 ( S , V , E , I , H , R ) H = 1 + c 4 1 H , V 2 ( S , V , E , I , H , R ) R = 1 + c 4 1 R .
It is very easy to obtain that V 2 has a one point of stagnation
( S ( 0 ) , V ( 0 ) , E ( 0 ) , I ( 0 ) , H ( 0 ) , R ( 0 ) ) = 1 + c 1 c 4 1 + c 4 , 1 1 + c 4 , c 2 c 4 1 + c 4 , c 3 c 4 1 + c 4 , 1 1 + c 4 , 1 1 + c 4 .
Next, the hessian matrix of V 2 ( S , V , E , I , H , R ) at ( S ( 0 ) , V ( 0 ) , E ( 0 ) , I ( 0 ) , H ( 0 ) , R ( 0 ) ) is
B = 1 + c 1 c 4 S 2 ( 0 ) 0 0 0 0 0 0 1 V 2 ( 0 ) 0 0 0 0 0 0 c 2 c 4 E 2 ( 0 ) 0 0 0 0 0 0 c 3 c 4 I 2 ( 0 ) 0 0 0 0 0 0 1 H 2 ( 0 ) 0 0 0 0 0 0 1 R 2 ( 0 ) .
Definitely, the matrix of the Hessian is definitely positive. Hence, V 2 ( S , V , E , I , H , R ) must have a small value V 2 ( S ( 0 ) , V ( 0 ) , E ( 0 ) , I ( 0 ) , H ( 0 ) , R ( 0 ) ) . By Equation (23) and by the assumption of the continuity of V 2 ( S , V , E , I , H , R ) , we write that V 2 ( S , V , E , I , H , R ) has one and at least one constant point V 2 ( S ( 0 ) , V ( 0 ) , E ( 0 ) , I ( 0 ) , H ( 0 ) , R ( 0 ) ) in R + 6 .
Next, we give the formula for a non-negative C 2 operator V : R + 6 R + as under
V ( S , V , E , I , H , R ) = V 2 ( S , V , E , I , H , R ) V 2 ( S ( 0 ) , V ( 0 ) , E ( 0 ) , I ( 0 ) , H ( 0 ) , R ( 0 ) ) .
By the application of I t o s formula to the said problem, we obtain
L ( V ) c 4 4 Π ( R ˜ 0 S ) 1 / 4 1 Π S + η I N + ( κ + δ ) + Φ 1 2 2 κ S V + ( 1 τ ) η I N + δ + Φ 2 2 2 ρ 1 I H + μ 1 + μ 2 + δ + Φ 5 2 2 ρ 2 I R μ 1 H R + δ + Φ 6 2 2 ρ 3 I μ 2 H + Π δ ( S ( t ) + V ( t ) + E ( t ) + I ( t ) + H ( t ) + R ( t ) ) .
leading to the assertion as follows:
L ( V ) c 4 c 5 Π S + η I N + ( κ + δ ) + Φ 1 2 2 κ S V + ( 1 τ ) η I N + δ + Φ 2 2 2 ρ 1 I H + μ 1 + μ 2 + δ + Φ 5 2 2 ρ 2 I R μ 1 H R + δ + Φ 6 2 2 ρ 3 I μ 2 H + Π δ N .
where
C 5 = 4 Π ( R 0 S ) 1 / 4 1 > 0 .
Next we formulate
D = { ϵ 1 < S < 1 ϵ 2 , ϵ 1 < V < 1 ϵ 2 , ϵ 1 < E < 1 ϵ 2 , ϵ 1 < I < 1 ϵ 2 ϵ 1 < H < 1 ϵ 2 ϵ 1 < R < 1 ϵ 2 } ,
where ϵ i > 0 for ( i = 1 , 2 , , 12 ) are minimum constant values to be computed later on. For simplification, we can write R + 6 D in the domain as follows:
D 1 = ( S , V , E , I , H , R ) R + 5 , 0 < S ϵ 1 , D 2 = ( S , V , E , I , H , R ) R + 5 , 0 < V ϵ 1 , S > ϵ 2 , D 3 = ( S , V , E , I , H , R ) R + 5 , 0 < E ϵ 1 , V > ϵ 2 , D 4 = ( S , V , E , I , H , R ) R + 5 , 0 < I ϵ 1 , E > ϵ 2 , D 5 = ( S , V , E , I , H , R ) R + 5 , 0 < H ϵ 1 , I > ϵ 2 , D 6 = ( S , V , E , I , H , R ) R + 5 , 0 < R ϵ 1 , H > ϵ 2 , D 7 = ( S , V , E , I , H , R ) R + 5 , S 1 ϵ 2 , D 8 = ( S , V , E , I , H , R ) R + 5 , V 1 ϵ 2 , D 9 = ( S , V , E , I , H , R ) R + 5 , E 1 ϵ 2 , D 10 = ( S , V , E , I , H , R ) R + 5 , I 1 ϵ 2 , D 11 = ( S , V , E , I , H , R ) R + 5 , H 1 ϵ 2 , D 12 = ( S , V , E , I , H , R ) R + 5 , R 1 ϵ 2 .
Going ahead, we have to derive L V ( S , V , E , I , H , R ) < 0 on R + 6 D , which is similar to how it is shown on the above-cited eight areas.
Case 1. If ( S , V , E , I , H , R ) D 1 , then by Equation (27), we obtain
L V c 4 c 5 Π S + η I N + ( κ + δ ) + Φ 1 2 2 κ S V + ( 1 τ ) η I N + δ + Φ 2 2 2 ρ 1 I H + μ 1 + μ 2 + δ + Φ 5 2 2 ρ 2 I R μ 1 H R + δ + Φ 6 2 2 ρ 3 I μ 2 H + Π δ N , c 4 c 5 + η I N + ( κ + δ ) + ( 1 τ ) η I N + δ + μ 1 + μ 2 + δ + δ + Π + Φ 1 2 + Φ 2 2 + Φ 5 2 + Φ 6 2 2 δ ϵ 1 .
Setting ϵ 1 > 0 , gives L V < 0 for every ( S , V , E , I , H , R ) D 1 .
Case 2. If ( S , V , E , I , H , R ) D 2 , then from Equation (27), we can get
L V c 4 c 5 Π S + η I N + ( κ + δ ) + Φ 1 2 2 κ S V + ( 1 τ ) η I N + δ + Φ 2 2 2 ρ 1 I H + μ 1 + μ 2 + δ + Φ 5 2 2 ρ 2 I R μ 1 H R + δ + Φ 6 2 2 ρ 3 I μ 2 H + Π δ N , c 4 c 5 + η I N + ( κ + δ ) + ( 1 τ ) η I N + δ + μ 1 + μ 2 + δ + δ + Π + Φ 1 2 + Φ 2 2 + Φ 5 2 + Φ 6 2 2 Π ϵ 1 .
We choose maximally high c 4 > 0 and maximally low ϵ 1 > 0 , so we can obtain L V < 0 for any ( S , V , E , I , H , R ) D 2 .
Case 3. If ( S , V , E , I , H , R ) D 3 , then from Equation (27), we obtain
L V c 4 c 5 Π S + η I N + ( κ + δ ) + Φ 1 2 2 κ S V + ( 1 τ ) η I N + δ + Φ 2 2 2 ρ 1 I H + μ 1 + μ 2 + δ + Φ 5 2 2 ρ 2 I R μ 1 H R + δ + Φ 6 2 2 ρ 3 I μ 2 H + Π δ N , c 4 c 5 + η I N + ( κ + δ ) + ( 1 τ ) η I N + δ + μ 1 + μ 2 + δ + δ + Π + Φ 1 2 + Φ 2 2 + Φ 5 2 + Φ 6 2 2 S ϵ 1 .
Selecting small ϵ 1 > 0 , thus, we have L V < 0 for every ( S , V , E , I , H , R ) D 3 .
Case 4. I f ( S , V , E , I , H , R ) D 4 , from Equation (27), we obtain
L V c 4 c 5 Π S + η I N + ( κ + δ ) + Φ 1 2 2 κ S V + ( 1 τ ) η I N + δ + Φ 2 2 2 ρ 1 I H + μ 1 + μ 2 + δ + Φ 5 2 2 ρ 2 I R μ 1 H R + δ + Φ 6 2 2 ρ 3 I μ 2 H + Π δ N , c 4 c 5 + η I N + ( κ + δ ) + ( 1 τ ) η I N + δ + μ 1 + μ 2 + δ + δ + Π + Φ 1 2 + Φ 2 2 + Φ 5 2 + Φ 6 2 2 S ϵ 1 .
Select small ϵ 2 > 0 to obtain L V < 0 for each ( S , V , E , I , H , R ) D 4 .
Case 5. If ( S , V , E , I , H , R ) D 5 , from Equation (27), we obtain
L V c 4 c 5 Π S + η I N + ( κ + δ ) + Φ 1 2 2 κ S V + ( 1 τ ) η I N + δ + Φ 2 2 2 ρ 1 I H + μ 1 + μ 2 + δ + Φ 5 2 2 ρ 2 I R μ 1 H R + δ + Φ 6 2 2 ρ 3 I μ 2 H + Π δ N , c 4 c 5 + η I N + ( κ + δ ) + ( 1 τ ) η I N + δ + μ 1 + μ 2 + δ + δ + Π + Φ 1 2 + Φ 2 2 + Φ 5 2 + Φ 6 2 2 δ ϵ 2 .
We choose small ϵ 2 > 0 , so we can obtain L V < 0 for any ( S , V , E , I , H , R ) D 5 .
Case 6. If ( S , V , E , I , H , R ) D 6 , from Equation (27), we obtain
L V c 4 c 5 Π S + η I N + ( κ + δ ) + Φ 1 2 2 κ S V + ( 1 τ ) η I N + δ + Φ 2 2 2 ρ 1 I H + μ 1 + μ 2 + δ + Φ 5 2 2 ρ 2 I R μ 1 H R + δ + Φ 6 2 2 ρ 3 I μ 2 H + Π δ N , c 4 c 5 + η I N + ( κ + δ ) + ( 1 τ ) η I N + δ + μ 1 + μ 2 + δ + δ + Π + Φ 1 2 + Φ 2 2 + Φ 5 2 + Φ 6 2 2 δ ϵ 2 .
We can choose sufficiently small ϵ 2 > 0 , so we can obtain L V < 0 for any ( S , V , E , I , H , R ) D 6 .
Case 7. If ( S , V , E , I , H , R ) D 7 , from Equation (27), we obtain
L V c 4 c 5 Π S + η I N + ( κ + δ ) + Φ 1 2 2 κ S V + ( 1 τ ) η I N + δ + Φ 2 2 2 ρ 1 I H + μ 1 + μ 2 + δ + Φ 5 2 2 ρ 2 I R μ 1 H R + δ + Φ 6 2 2 ρ 3 I μ 2 H + Π δ N , c 4 c 5 + η I N + ( κ + δ ) + ( 1 τ ) η I N + δ + μ 1 + μ 2 + δ + δ + Π + Φ 1 2 + Φ 2 2 + Φ 5 2 + Φ 6 2 2 Π ϵ 2 .
By considering the smallest value of ϵ 2 > 0 , we can obtain L V < 0 for any ( S , V , E , I , H , R ) D 7 .
Case 8. If ( S , V , E , I , H , R ) D 8 , from Equation (27), we obtain
L V c 4 c 5 Π S + η I N + ( κ + δ ) + Φ 1 2 2 κ S V + ( 1 τ ) η I N + δ + Φ 2 2 2 ρ 1 I H + μ 1 + μ 2 + δ + Φ 5 2 2 ρ 2 I R μ 1 H R + δ + Φ 6 2 2 ρ 3 I μ 2 H + Π δ N , c 4 c 5 + η I N + ( κ + δ ) + ( 1 τ ) η I N + δ + μ 1 + μ 2 + δ + δ + Π + Φ 1 2 + Φ 2 2 + Φ 5 2 + Φ 6 2 2 κ ϵ 1 ϵ 1 .
Let us select the smallest value of ϵ 1 , ϵ 2 > 0 so that we can obtain L V < 0 for any ( S , V , E , I , H , R ) D 8 .
Case 9. If ( S , V , E , I , H , R ) D 9 , from Equation (27), we obtain
L V c 4 c 5 Π S + η I N + ( κ + δ ) + Φ 1 2 2 κ S V + ( 1 τ ) η I N + δ + Φ 2 2 2 ρ 1 I H + μ 1 + μ 2 + δ + Φ 5 2 2 ρ 2 I R μ 1 H R + δ + Φ 6 2 2 ρ 3 I μ 2 H + Π δ N , c 4 c 5 + η I N + ( κ + δ ) + ( 1 τ ) η I N + δ + μ 1 + μ 2 + δ + δ + Π + Φ 1 2 + Φ 2 2 + Φ 5 2 + Φ 6 2 2 δ ϵ 2 .
If ϵ 2 > 0 , then we can find L V < 0 for each ( S , V , E , I , H , R ) D 9 .
Case 10. If ( S , V , E , I , H , R ) D 10 , from Equation (27), we obtain
L V c 4 c 5 Π S + η I N + ( κ + δ ) + Φ 1 2 2 κ S V + ( 1 τ ) η I N + δ + Φ 2 2 2 ρ 1 I H + μ 1 + μ 2 + δ + Φ 5 2 2 ρ 2 I R μ 1 H R + δ + Φ 6 2 2 ρ 3 I μ 2 H + Π δ N , c 4 c 5 + η I N + ( κ + δ ) + ( 1 τ ) η I N + δ + μ 1 + μ 2 + δ + δ + Π + Φ 1 2 + Φ 2 2 + Φ 5 2 + Φ 6 2 2 δ ϵ 2 .
If ϵ 2 > 0 , then we can find L V < 0 for each ( S , V , E , I , H , R ) D 10 .
Case 11. If ( S , V , E , I , H , R ) D 11 , from Equation (27), we obtain
L V c 4 c 5 Π S + η I N + ( κ + δ ) + Φ 1 2 2 κ S V + ( 1 τ ) η I N + δ + Φ 2 2 2 ρ 1 I H + μ 1 + μ 2 + δ + Φ 5 2 2 ρ 2 I R μ 1 H R + δ + Φ 6 2 2 ρ 3 I μ 2 H + Π δ N , c 4 c 5 + η I N + ( κ + δ ) + ( 1 τ ) η I N + δ + μ 1 + μ 2 + δ + δ + Π + Φ 1 2 + Φ 2 2 + Φ 5 2 + Φ 6 2 2 ρ 1 I ϵ 2 .
If ϵ 2 > 0 , then we can find L V < 0 for each ( S , V , E , I , H , R ) D 11 .
Case 12. If ( S , V , E , I , H , R ) D 12 , from Equation (27), we obtain
L V c 4 c 5 Π S + η I N + ( κ + δ ) + Φ 1 2 2 κ S V + ( 1 τ ) η I N + δ + Φ 2 2 2 ρ 1 I H + μ 1 + μ 2 + δ + Φ 5 2 2 ρ 2 I R μ 1 H R + δ + Φ 6 2 2 ρ 3 I μ 2 H + Π δ N , c 4 c 5 + η I N + ( κ + δ ) + ( 1 τ ) η I N + δ + μ 1 + μ 2 + δ + δ + Π + Φ 1 2 + Φ 2 2 + Φ 5 2 + Φ 6 2 2 ρ 1 ϵ 2 ϵ 1 .
For the smallest value of ϵ 1 , ϵ 2 > 0 , we can obtain L V < 0 for any ( S , V , E , I , H , R ) D 12 .
So, we reach the concluding remarks that here is a fixed W > 0 :
L V ( S , V , E , I , H , R ) < W < 0 for all ( S , V , E , I , H , R ) R + 6 D .
Therefore,
d V ( S , V , E , I , R ) < W d t + [ ( c 4 + 1 ) S ( c 1 c 4 + 1 ) ζ 1 ] d B 1 ( t ) + [ ( c 4 + 1 ) V ζ 2 ] d B 2 ( t ) + [ ( c 4 + 1 ) E c 2 c 4 ζ 3 ] d B 3 ( t ) + [ ( c 4 + 1 ) I c 3 c 4 ζ 4 ] d B 4 ( t ) + [ ( c 4 + 1 ) H ζ 5 ] d B 5 ( t ) + [ ( c 4 + 1 ) R ζ 5 ] d B 6 ( t ) .
Assume that ( S ( 0 ) , V ( 0 ) , E ( 0 ) , I ( 0 ) , H ( 0 ) , R ( 0 ) ) = ( x 1 , x 2 , x 3 , x 4 , x 5 , x 6 ) = x R + 6 D , and τ x is that time at which a path starting from x reaches set D,
τ n = i n f { t : | X ( t ) | = n } and τ ( n ) ( t ) = min { τ x , t , τ n } .
By taking the integral of inequality (28) from zero to τ ( n ) ( t ) , considering the expected result, and by application of Dynkin’s expression, we obtain
E V ( S ( τ ( n ) ( t ) ) , V ( τ ( n ) ( t ) ) , E ( τ ( n ) ( t ) ) , I ( τ ( n ) ( t ) ) , H ( τ ( n ) ( t ) ) , R ( τ ( n ) ( t ) ) ) V ( x ) = E 0 τ ( n ) ( t ) L V ( S ( u ) , V ( u ) , E ( u ) , I ( u ) , H ( u ) , R ( u ) ) d u , E 0 τ ( n ) ( t ) W d u = W E τ ( n ) ( t ) .
As V ( x ) > 0 , so
E τ ( n ) ( t ) V ( x ) W .
for the proof of (3), we obtain P { τ e = } = 1 . Alternatively, we can say that the model (2) is continuous. So, we write that as t and n , then one has τ ( n ) ( t ) τ x almost surely. □
By Fatou’s Lemma,
E τ ( n ) ( t ) V ( x ) W < .
Thus, s u p x K E τ x < , where K is a compact subset of R + 6 . It is proved in a direct way by the result (ii) of Lemma 3.
Furthermore, the matrix of diffusion model (2) is as follows:
B = ξ 1 2 S 2 0 0 0 0 0 0 ξ 2 2 V 2 0 0 0 0 0 0 ξ 3 2 E 2 0 0 0 0 0 0 ξ 4 2 I 2 0 0 0 0 0 0 ξ 5 2 H 2 0 0 0 0 0 0 ξ 6 2 R 2 .
Choosing
M = min ( S , V , E , I , H , R ) D ¯ R + 6 { ξ 1 2 S 2 , ξ 2 2 V 2 , ξ 3 2 E 2 , ξ 4 2 I 2 , ξ 5 2 H 2 , ξ 6 2 R 2 } ,
we obtain
i , j = 1 6 a i j ( S , V , E , I , H , R ) ξ i ξ j = ξ 1 2 S 2 ζ 1 2 + ξ 2 2 V 2 ζ 2 2 + ξ 3 2 E 2 ζ 2 + ξ 4 2 I 2 ζ 4 2 + ξ 5 2 ζ 5 2 + ξ 6 2 R 2 ζ 6 2 M | ξ | 2 , ( S , V , E , I , H , R ) D ¯ ,
here ξ = ( ζ 1 , ζ 2 , ζ 3 , ζ 4 , ζ 5 , ζ 6 ) R + 6 . This implies that Equation (1) of Lemma 3 also fulfilled. From this, we conclude that Lemma 3 is sure for model (2), which is ergodic and on stationary distribution.

7. Numerical Simulations

For the validation of our obtained scheme, we establish the graphical representation using the approximate scheme to the model (2). The graphical representation depends on the qualitative analysis of the discussion, and the numerics of the used parameters are epidemiologically justifiable. Further, for the use of a stochastic approach of the RK method of the 4th order, the discretization of model (2) is
S i + 1 = S i + Π η S i I i N ( κ + δ ) S i t + ξ 1 S i t ζ 1 , i + ξ 1 2 2 S i ( ζ 1 , i 2 1 ) t , V i + 1 = V i + S i ( 1 τ ) η V i I i N δ V i t + ξ 2 V i t ζ 2 , i + ξ 2 2 2 V i ( ζ 2 , i 2 1 ) t , E i + 1 = E i + η S i I i N + ( 1 τ ) η V i I i N ( σ + δ ) E i t + ξ 3 E i t ζ 3 , i + ξ 3 2 2 E i ( ζ 3 , i 2 1 ) t , I i + 1 = I i + σ E i ( ρ 1 + ρ 2 + ρ 3 + δ ) I i t + ξ 4 I i t ζ 4 , i + ξ 4 2 2 I i ( ξ 4 , i 2 1 ) t , H i + 1 = H i + ρ 1 H i ( μ 1 + μ 2 + δ ) H i t + ξ 4 H i t ζ 4 , i + ξ 5 2 2 I i ( ξ 5 , i 2 1 ) t , R i + 1 = R i + ρ 2 I I + μ 1 H I μ R i t + ξ 6 R i t ζ 6 , i + ξ 6 2 2 R i ( ζ 6 , i 2 1 ) t .
ζ i , j ( i = 1 , 2 , 3 , 4 , 5 ) , is normal distribution satisfying the division of N ( 0 , 1 ) , and the difference Δ t . Take ξ i > 0 , ( i = 1 , 2 , 3 , 4 , 5 ) as the white noise intensities.
For the well dynamics, the stability of the stochastically described model, and optimal controlling, we must point out the parameters values for the numerical simulations of (2).
Now here, we give the discussion of the graphical representation and epidemiological feasibility of the proposed problem (2). For this, we use the numerical parameters values of Table 2 ( S 1 ). The initial values for all the agents are also written Table 2 ( S 1 ). Using Theorem 2 gives the conditions for vanishing the epidemic. Few of the achievements are obtained from the analysis of stability in the stochastically investigated model. Theorem 2 is satisfied if the reproduction value is R 0 < 1 . This shows that the epidemic vanishing probability will be one in the stochastic model. Similarly, if R 0 D < 1 , the idealistic problem (2) will be asymptotic and globally stable. For both curves, they converge to the free critical value. This is given in Figure 1, which implies that COVID-19 infection will die out of society, given some important conditions.The simulations of all the agents are provided in Figure 1a–e.
Next, we provide the digital findings for the stationary points of distributions. In Lemma 3, an ergodic condition is used to prove that the system is more realistic in the stochastic approach and has one stationary distribution. The stochastically given problem (2) is proposed for the parameter values given in Table 2  ( S 2 ) and evaluated as R 0 s > 1 , so using Theorem 3, which is fulfilled, because of small intensities of white noises, the infection reflection will lie. We observe in Figure 2 that the infected system (2) will lie or remain in the mean, which may validate the output of Theorem 3 which implies that the system (2) lies on ergodic stationary distributions. The Susceptible, Vaccinated, Expose, Infected, and Recovered individual stationary distribution numerical simulations can be clearly seen in Figure 2a–f, respectively. Theorem 3 implies that model (2) must be ergodic stationary distributions. Figure 3 confirms this.

8. Conclusions

In this article, we discussed the asymptotic dynamics of a stochastic COVID-19 epidemic model with a general incident rate and the effect of vaccination. Firstly, we analyzed the considered model for the unique globally non-negative root along with an initial approximation. The stability of the solution of the given model was also computed by the help of the Lyapunov operator. For eliminating the infection from society, we derived the reproduction value R 0 < 1 . The Lyapunov function method proved that an ergodic stationary division for the non-negative root of the given system exists and is unique, which shows that for R 0 s > 1 , the disease may lie in the community. Few of the numerical simulations have been performed by the R K 4 approach for the validation of the obtained results. We also will study in our future research work the effect of Levy noises on the dynamics of a more complex population system.

Author Contributions

Conceptualization, H.B. and Z.S.; Formal analysis, H.B. and W.L.; Funding acquisition, Z.S.; Investigation, H.B.; Methodology, H.B. and W.L.; Software, H.B.; Validation, H.B.; Writ draft, H.B. and W.L.; Writ editing, Z.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grant 11801249; the Nature Science Foundation of Shandong Province China under Grant ZR2020MF040, and in part by the Open Project of Liaocheng University under Grant 319312101-01.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China under Grant 11801249; the Nature Science Foundation of Shandong Province China under Grant ZR2020MF040, and in part by the Open Project of Liaocheng University under Grant 319312101-01.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Huang, C.; Wang, Y.; Li, X.; Ren, L.; Zhao, J.; Hu, Y.; Zhang, L. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 2020, 395, 497–506. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ding, Y.; Jiao, J.; Zhang, Q.; Zhang, Y.; Ren, X. Stationary Distribution and Extinction in a Stochastic SIQR Epidemic Model Incorporating Media Coverage and Markovian Switching. Symmetry 2021, 13, 1122. [Google Scholar] [CrossRef]
  3. Atangana, A. Modelling the spread of COVID-19 with new fractal-fractional operators: Can the lockdown save mankind before vaccination? Chaos Solitons Fractals 2020, 136, 109860. [Google Scholar] [CrossRef] [PubMed]
  4. Tul, A.Q.; Anjum, N.; Zeb, A.; Djilali, S.; Khan, Z.A. On the analysis of Caputo fractional order dynamics of Middle East Lungs Coronavirus (MERS-CoV) model. Alex. Eng. J. 2022, 61, 5123–5131. [Google Scholar]
  5. Ali, I.; Khan, S.U. Threshold of Stochastic SIRS Epidemic Model from Infectious to Susceptible Class with Saturated Incidence Rate Using Spectral Method. Symmetry 2022, 14, 1838. [Google Scholar] [CrossRef]
  6. Anwarud, D.; Li, Y.; Shah, M.A. The complex dynamics of hepatitis B infected individuals with optimal control. J. Syst. Sci. Complex. 2021, 34, 1301–1323. [Google Scholar]
  7. Allen, L.J.S. An introduction to stochastic epidemic models. In Mathematical Epidemiology; Springer: Berlin/Heidelberg, Germany, 2008; pp. 81–130. [Google Scholar]
  8. Lei, Q.; Yang, Z. Dynamical behaviours of a stochastic SIRI epidemic model. Appl. Anal. 2016, 96, 1–13. [Google Scholar]
  9. Özdemir, N.; Agrawal, O.P.; Karadeniz, D.; İskender, B.B. Fractional optimal control problem of an axis-symmetric diffusion-wave propagation. Phys. Scr. 2009, T136, 014024. [Google Scholar] [CrossRef]
  10. Din, A.; Yassine, S. Long-term bifurcation and stochastic optimal control of a triple-delayed Ebola illness model with vaccination and quarantine strategies. Fractal Fract. 2022, 6, 578. [Google Scholar] [CrossRef]
  11. Khan, A.; Yassine, S. Stochastic modeling of the Monkeypox 2022 epidemic with cross-infection hypothesis in a highly disturbed environment. Math. Biosci. Eng. 2022, 19, 13560–13581. [Google Scholar] [CrossRef] [PubMed]
  12. Din, A.; Li, Y.; Tahir, K.; Gul, Z. Mathematical analysis of spread and control of the novel corona virus (COVID-19) in China. Chaos Solitons Fractals 2020, 141, 110286. [Google Scholar] [CrossRef] [PubMed]
  13. Din, A.; Li, Y.; Abdullahi, Y. Delayed hepatitis B epidemic model with stochastic analysis. Chaos Solitons Fractals 2021, 146, 110839. [Google Scholar] [CrossRef]
  14. Yassine, S.; Khan, A.; Din, A.; Kiouach, D.; Rajasekar, S.P. Determining the global threshold of an epidemic model with general interference function and high-order perturbation. AIMS Math. 2022, 7, 19865–19890. [Google Scholar]
  15. Alzahrani, E.O.; Khan, M.A. Modeling the dynamics of Hepatitis E with optimal control. Chaos Solitons Fractals 2018, 116, 287–301. [Google Scholar] [CrossRef]
  16. Olver, P.J. Applications of Lie Groups to Differential Equations; Springer: Berlin/Heidelberg, Germany, 2000. [Google Scholar]
  17. Olver, P.J. Equivalence, Invariants, and Symmetry; Cambridge University Press: Cambridge, UK, 1995. [Google Scholar]
  18. Zhang, T.; Ding, T.; Gao, N.; Song, Y. Dynamical behavior of a stochastic SIRC model for influenza A. Symmetry 2020, 12, 745. [Google Scholar] [CrossRef]
  19. Gaeta, G.; Rodriguez-Quintero, N. Lie-point symmetries and stochastic differential equations. J. Phys. A Math. Gen. 1999, 32, 8485–8505. [Google Scholar] [CrossRef]
  20. Gaeta, G. Lie-point symmetries and stochastic differential equations: II. J. Phys. A Math. Gen. 2000, 33, 4883–4902. [Google Scholar] [CrossRef]
  21. Zhang, X.-B.; Wang, X.-D.; Huo, H. Extinction and stationary distribution of a stochastic SIRS epidemic model with standard incidence rate and partial immunity. Phys. Stat. Mech. Appl. 2019, 531, 121548. [Google Scholar] [CrossRef]
  22. Din, A.; Saida, A.; Amina, A. A stochastically perturbed co-infection epidemic model for COVID-19 and hepatitis B virus. Nonlinear Dyn. 2022, 28, 1–25. [Google Scholar] [CrossRef] [PubMed]
  23. Khasminskii, R. Stochastic Stability of Differential Equations; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2011; Volume 66. [Google Scholar]
Figure 1. Solution paths of the model (2) when the numerical values are taken as shown in the 2nd column of Table 2.
Figure 1. Solution paths of the model (2) when the numerical values are taken as shown in the 2nd column of Table 2.
Symmetry 15 00285 g001
Figure 2. Solutions paths of the model (2) when the numerical values are taken as shown in the third column of Table 2.
Figure 2. Solutions paths of the model (2) when the numerical values are taken as shown in the third column of Table 2.
Symmetry 15 00285 g002
Figure 3. Ergodic stationary distribution of model (2).
Figure 3. Ergodic stationary distribution of model (2).
Symmetry 15 00285 g003
Table 1. Description of the parameters.
Table 1. Description of the parameters.
ParameterDescription
Π Rate of recruitment.
η Rate of infection effectively
κ Vaccinated population in percentage.
τ effect of Vaccination
δ Rate of natural death
σ Rate of sign reported by lab
ρ 2 Recovery rate from I
ρ 3 COVID-19 death rate
ρ 1 Transferred rate from I to to H
μ 2 COVID-19 death rate
μ 1 recovered rate of H
Table 2. Parameter values used in the simulation of model (2).
Table 2. Parameter values used in the simulation of model (2).
Parameters S 1 S 2 Source
Π 1.503.50assumed
τ 0.020.30assumed
η 0.070.03assumed
δ 0.010.03assumed
σ 0.010.05assumed
κ 0.020.04assumed
ρ 1 0.050.10assumed
ρ 2 0.350.20assumed
ρ 3 0.050.30assumed
μ 1 0.550.40assumed
μ 2 0.150.50assumed
S ( 0 ) 50.04.00assumed
V ( 0 ) 20.03.00assumed
E ( 0 ) 30.01.00assumed
I ( 0 ) 40.02.00assumed
H ( 0 ) 40.02.00assumed
R ( 0 ) 10.01.00assumed
Φ 1 1.251.20assumed
Φ 2 1.231.25assumed
Φ 3 1.351.15assumed
Φ 4 1.201.05assumed
Φ 5 1.151.22assumed
Φ 6 1.101.15assumed
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Batool, H.; Li, W.; Sun, Z. Extinction and Ergodic Stationary Distribution of COVID-19 Epidemic Model with Vaccination Effects. Symmetry 2023, 15, 285. https://doi.org/10.3390/sym15020285

AMA Style

Batool H, Li W, Sun Z. Extinction and Ergodic Stationary Distribution of COVID-19 Epidemic Model with Vaccination Effects. Symmetry. 2023; 15(2):285. https://doi.org/10.3390/sym15020285

Chicago/Turabian Style

Batool, Humera, Weiyu Li, and Zhonggui Sun. 2023. "Extinction and Ergodic Stationary Distribution of COVID-19 Epidemic Model with Vaccination Effects" Symmetry 15, no. 2: 285. https://doi.org/10.3390/sym15020285

APA Style

Batool, H., Li, W., & Sun, Z. (2023). Extinction and Ergodic Stationary Distribution of COVID-19 Epidemic Model with Vaccination Effects. Symmetry, 15(2), 285. https://doi.org/10.3390/sym15020285

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop