Next Article in Journal
Canonical Coordinates and Natural Equation for Lorentz Surfaces in R13
Next Article in Special Issue
The Stability of Functional Equations with a New Direct Method
Previous Article in Journal
Global and Local Analysis for a Cournot Duopoly Game with Two Different Objective Functions
Previous Article in Special Issue
Existence and Uniqueness of Nontrivial Periodic Solutions to a Discrete Switching Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of the Stochastic Modeling of COVID-19 with Environmental Noise from the Analytical and Numerical Point of View

1
Faculty of Informatics and Computing, Besut Campus, Universiti Sultan Zainal Abidin (UniSZA), Besut 22200, Terengganu, Malaysia
2
Department of Mathematics, Shaheed Benazir Bhutto Univeresity, Dir Upper, Sheringal 18050, Khyber Pakhtunkhwa, Pakistan
3
Department of Mathematics, Azarbaijan Shahid Madani University, Tabriz 53751-71379, Iran
4
Department of Medical Research, China Medical University Hospital, China Medical University, Taichung 40402, Taiwan
5
Mathematics Department, Faculty of Science and Technology, Suan Dusit University, Bangkok 10300, Thailand
6
Department of Mathematics, Faculty of Applied Science, King Mongkut’s University of Technology North Bangkok, Bangkok 10800, Thailand
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2021, 9(23), 3122; https://doi.org/10.3390/math9233122
Submission received: 7 November 2021 / Revised: 26 November 2021 / Accepted: 30 November 2021 / Published: 3 December 2021
(This article belongs to the Special Issue Functional Differential Equations and Epidemiological Modelling)

Abstract

:
In this article, we propose a novel mathematical model for the spread of COVID-19 involving environmental white noise. The new stochastic model was studied for the existence and persistence of the disease, as well as the extinction of the disease. We noticed that the existence and extinction of the disease are dependent on R 0 (the reproduction number). Then, a numerical scheme was developed for the computational analysis of the model; with the existing values of the parameters in the literature, we obtained the related simulations, which gave us more realistic numerical data for the future prediction. The mentioned stochastic model was analyzed for different values of σ 1 , σ 2 and β 1 , β 2 , and both the stochastic and the deterministic models were compared for the future prediction of the spread of COVID-19.

1. Introduction and Backgrounds

Coronavirus disease 2019 (COVID-19) is a communicable respiratory disease. The disease is caused by a newly discovered virus strain, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) [1]. COVID-19 was first identified in December 2019 in Wuhan, China, and spread quickly over four months. In a short period, more than 2.9 million people in 185 nations around the world were infected and 206 thousand people had died [2]. On March 11, 2020, The World Health Organization (WHO) confirmed that the spread of this disease constituted a pandemic [3]. This disease can be spread primarily from droplets produced when coughing or sneezing, by person-to-person contact, or even through conversation. By contacting contaminated surfaces, susceptible individuals can also become infected. The most prevalent signs of this disease are fever, nausea, dry cough, fatigue, and shortness of breath. All these signs constitute COVID-19 [4]. Some patients can have joint pain, nasal congestion, diarrhea, runny nose, or sore throat. The symptoms are typically mild, but can slowly worsen. In order to prevent infection, frequent hand washing, covering the nose and mouth when sneezing or coughing, avoiding touching the nose, mouth, or eyes, and other preventive steps are advised, including social distancing.
Due to the seriousness of the COVID-19 pandemic, many countries have taken extreme measures to curb its spread. In addition, they have evaluated and adapted their healthcare systems accordingly. Hence, they have canceled public events, closed public venues, schools, and borders, placed restrictions on travel, implemented lockdown measures, etc. While those measures have been helpful, lockdowns have led to negative socioeconomic impacts such as the bankruptcy of businesses, the loss of employment, etc. Furthermore, lockdowns have disrupted supply chains and decreased productivity. The shutdown of China’s drug-producing plants, which constitute the second largest pharmaceutical product exporters, has delayed the delivery of generic drugs [5]. The air transport, tourism, and oil sectors have been visibly affected. It is also expected that there will be unforeseen impacts irrespective of the duration of the pandemic. According to The International Monetary Fund, the worldwide economy was expected to shrink by 3% in 2020 [6].
Governments must prevent the failure of the economy, taking measures to relax the lockdowns while preventing the spread of the disease. Some developed countries intend to or have already issued immunity identity cards. However, The World Health Organization has disapproved of this technique, since there is a lack of adequate scientific proof regarding reinfection, possibly rendering such measures ineffective. A risk-balancing strategy has been adopted by the South African government to lift the lockout restrictions progressively. We refer the readers to [7,8] for some scientific works performed on infectious diseases and, in particular, for several mathematical models related to abstract and real-world boundary problems via fractional operators, to [9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24].
Brownian motion was discovered by the scholar Robert Brown in 1827. While Brown was studying pollen particles drifting in water in a microscope, he noticed that the particles exhibited an unsteady motion. After replicating the experiment, he concluded that the motion was due to the particles being alive; however, the initiation of the motion remained unexplained. The first to provide a theory of Brownian motion was Louis Bachelier in 1900 in his Ph.D. thesis. However, it was only in 1905 that Albert Einstein, using a probabilistic model, could adequately explain Brownian motion. He determined that the kinetic energy of fluids was the cause, in this case water molecules moving randomly. Thus, a small particle would receive a random number of impacts of random strength and from random directions in a short period of time. This random bombardment by the molecules of the fluid would cause a sufficiently small particle to move exactly as Brown described [25].
In probability theory and associated fields, a stochastic or random process is a mathematical object, usually characterized as a group of random variables. In the 1930s, the first mathematical definition of a stochastic process as a family of random variables indexed by means of the change of the real line was given by Aleksandr Khinchin. Khinchin, Andrey Kolmogorov, Joseph Doob, and William Feller furthered the probability concept and stochastic techniques. In mathematics, the theory of stochastic processes is a significant aspect of probability theory and a subject of research for both theory and applications. The word stochastic is used to portray different terms and objects in mathematics. Few example, this includes a stochastic matrix, which describes a stochastic process known as a Markov process, and stochastic calculus, which includes differential equations and integrals depending on stochastic processes such as the Wiener process, also called the Brownian motion (white noise) process.
Stochastic processes are broadly used as mathematical models of systems and phenomena that appear to vary in a random manner, for example the growth of a bacterial population, an electrical flow fluctuating because of thermal noise, or the moment of gas particles [26]. Stochastic processes have applications in several fields such as biology, chemistry, physics, image processing, computer science, and cryptography. Stochastic processes consist of the Wiener process or Brownian motion technique, utilized by Louis Bachelier to study the value changes on the Paris Bourse, and the Poisson technique, used by A.K. Erlang to study the number of phone calls happening in a specific time frame. Stochastic processes can be arranged into different classes, which include random walks, martingales, and Markov processes [27,28]. Mathematical knowledge on probability, calculus, linear algebra, set theory, and topology, as well as branches of mathematical analysis such as real analysis and measure theory are used in the study of stochastic techniques. The theory of stochastic processes is viewed as a significant contribution to mathematics, and it continues to be an active topic of research for both theoretical reasons and applications. The Wiener process is a stochastic process with stationary and independent addition, which are typically dispersed depending on the size of the increments. The Wiener process is named after Norbert Wiener, who demonstrated its mathematical existence, but the process is also called the Brownian motion (white noises) process or just Brownian motion because of its historical association as a model for Brownian motion in liquids. Brownian motion (white noises) refers to either the physical phenomenon by which minute particles immersed in a fluid move around randomly or the mathematical models used to describe those random movements [29].
Compartmental models are a type of general modeling tool that are frequently used in infectious diseases’ mathematical modeling. The population is divided into sections denoted by the letters S , I , or R (susceptible, infectious, or recovered). The flow patterns between the sections are commonly indicated by the order of the abbreviations; for example, SEIS stands for susceptible, exposed, infected, and susceptible. Mathematical modeling begin with the well-known contributions of Ross [30] in 1916, Ross and Hudson in 1917 [31,32], and Kermack and McKendrick in 1927 [33]. The work of Kermack and McKendrick, published in 1927, on inoculation against smallpox, had a major impact on the modeling of the system. Their SIR model is still used to model epidemics of infectious diseases. The SIR models are based on ordinary differential equations (which are deterministic), yet can likewise be utilized with a stochastic (random) structure, which is more sensible, yet substantially more challenging to investigate. The SIR model [34,35,36] is the basic approach for a dynamical study, which can then be extended with further subcompartments.
The S I R model comprises three classes: S are the susceptible people. At the point when a susceptible and an infectious person come into “contact”, the susceptible becomes infected and advances to the infected class. I is the class of infectious people. These are people who are infected and are able to infect uninfected people. R stands for those who were infected, but are now immune or are now deceased, thus entering the removed compartment. It is expected that the quantity of deceased individuals will be negligible concerning the absolute population. This compartment may likewise be classified as “recovered” or “safe”.
The S I R model predicts infections that are communicated from one person to another. These factors address the quantity of individuals in every compartment at a specific time. The quantity of susceptible, infected, and recovered people might differ over time (regardless of whether the absolute population size stays consistent), so we make the exact numbers a function of t (time): S ( t ) , I ( t ) , and R ( t ) . For a particular disease in a particular area, these compartments might be determined to anticipate possible outbreaks and bring them under control [37].
As suggested by the variable function of t, the model is dynamic in that the numbers in every class might change over time. The significance of this unique perspective is generally clear in an endemic disease with a short infectious period, for example measles in the U.K. before the introduction of immunization in 1968. Such sicknesses will often happen in patterns of episodes due to the variation in the number of susceptibles over time. During an outbreak, the quantity of susceptible people falls quickly as a greater amount of them are infected and consequently enter the infectious and removed compartments. An outbreak cannot make an appearance again until the susceptible population has increased again, for example because of a future generation being naturally introduced to the susceptible compartment. Such models are still under consideration by the research community, and one can find some related work in [38,39,40]. These aspects highlight some of the innovations in and contributions to the study of such models.
The notion of stochastic modeling can be investigated from different dimensions of engineering and science. For instance, stochastic models of financial markets are designed to analyze commodities, relative stock prices, and interest rates; in biological stochastic systems, we see that adding stochastic “noise” to internal links provides a balance and other noise interactions optimize function and behavior.It is believed that manufacturing processes have random structures. In fact, regarding both continuous and batch manufacturing processes, such a hypothesis is essentially valid. A process control chart represents a given parameter of process control in terms of time and is plotted to record and analyze the testing and monitoring of the relevant processes. Stochastic social science theory deals with subconscious processes. Since the number of involved variables is large, every event creates its own set of possibilities, which results in every event being unpredictable. The theory of stochastic models of different dynamical epidemic problems is of great interest among scientists. In the published literature, different aspects of theoretical analysis, including existence, extinction, persistence, and numerical simulations, can be observed in stochastic epidemic models. To better understand the subject, we recommend interested readers see [41,42,43,44,45] and the references therein.
In this article, we considered that the fluctuations of β 1 and β 2 are random and β 1 = β 1 + σ 1 d B 1 and β 2 = β 2 + σ 2 d B 2 . With these considerations, we have the following new stochastic model for the spread of COVID-19:
d S 1 = Λ 1 β 1 I S 1 α S 1 d t σ 1 I S 1 d B 1 , d S 2 = Λ 2 β 2 I S 2 α S 2 d t σ 2 I S 2 d B 2 , d I = I β 1 S 1 + β 2 S 2 ( μ + α ) I d t + I σ 1 S 1 d B 1 + σ 2 S 2 d B 2 , d H = μ I ( ρ + α ) H d t , d R = ρ H α R d t .
Here, S 1 represents the class of uninfected susceptible people, S 2 represents the class of older people or those having other diseases, I represents the class of people infected with COVID-19, H represents the class of people who are undergoing treatment, and R represents the class of recovered individuals. S 1 , S 2 are continuously increased with a constant ratio denoted by Λ 1 and Λ 2 , respectively. Furthermore, the parameters β 1 and β 2 are the ratios of the infection of S 1 and S 2 upon their contact with infected individuals I . The infected individuals shift to the treated class with a ratio of μ , and the recovery rate is ρ . σ 1 , σ 2 represent the intensity of environmental white noise, and B 1 , B 2 are the independent Brownian motions. For the case of σ 1 = σ 2 = 0 , we have the following deterministic model [46]:
d d t S 1 = Λ 1 β 1 I S 1 α S 1 , d d t S 2 = Λ 2 β 2 I S 2 α S 2 , d d t I = I β 1 S 1 + β 2 S 2 ( μ + α ) I , d d t H = μ I ( ρ + α ) H , d d t R = ρ H α R .
Our objectives were to study the analytical and computational approaches to the stochastic epidemic model (1). We compared our results with the deterministic model (2) numerically for different parametric values for the prediction of the future. Thus, the study can be helpful to the community regarding the awareness about the spread of the disease with a more realistic interpretation. Furthermore:
d N d t = ( Λ 1 + Λ 2 ) μ N ,
in which N ( t ) = S 1 ( t ) + H ( t ) + S 2 ( t ) + R ( t ) + I ( t ) shows the total constant population for ( Λ 1 + Λ 2 ) μ N and N ( 0 ) = S 1 ( 0 ) + S 2 ( 0 ) + I ( 0 ) + H ( 0 ) + R ( 0 ) .
The above equation has an exact solution, that is:
N ( t ) = e α t Ņ ( 0 ) + ( Λ 1 + Λ 2 ) α e α t ,
which implies that N = ( Λ 1 + Λ 2 ) α .
Furthermore, we have S 1 ( 0 ) 0 , I ( 0 ) 0 , S 2 ( 0 ) 0 , R ( 0 ) 0 , H ( 0 ) 0 , implying that S 1 ( t ) 0 , I ( t ) 0 S 2 ( t ) 0 , R ( t ) 0 , H ( t ) 0 . Therefore, the solution has a positive property.

2. Existence of the Solution

We made the following assumptions:
  • Set R + d = { ( h ˙ ¯ i R d , h ˙ ¯ i > 0 , 1 d } ;
  • Assume a complete probability space ( Ω ˙ , ˙ , { ˙ } t t 0 ) .
For the existence of the solution to our problem (1), we considered a general 5D stochastic differential equation, which takes the form:
d h ˙ ¯ ( t ) = h ( h ˙ ¯ ( t ) , t ) d t + g ( h ˙ ¯ ( t ) , t ) d B ( t ) f o r t t 0 ,
via the initial value h ˙ ¯ ( 0 ) = h ˙ ¯ 0 R + 5 . Furthermore, an operator L ˙ is defined as:
L ˙ = t + = 1 d h ( h ˙ ¯ , t ) h ˙ ¯ + 1 2 , j ̣ = 1 5 [ g T ( h ˙ ¯ , t ) g ( h ˙ ¯ , t ) ] j ̣ 2 h ˙ ¯ h ˙ ¯ j ̣ .
By operating L ˙ on a function V ˙ = ( R d × R ˙ + ; R ˙ + ) , then:
L ˙ V ˙ ( h ˙ ¯ , t ) = V ˙ t ( h ˙ ¯ , t ) + V ˙ h ˙ ¯ ( h ˙ ¯ , t ) f ( h ˙ ¯ , t ) + 1 2 t r a c e [ g T ( h ˙ ¯ , t ) V ˙ h ˙ ¯ h ˙ ¯ ( h ˙ ¯ , t ) g ( h ˙ ¯ , t ) ] .
The existence of the solution is ensured by the following lemma.
Lemma 1.
For ( S 1 ( 0 ) , H ( 0 ) , S 2 ( 0 ) , I ( 0 ) , R ( 0 ) ) R + 5 , a positive solution of (1) exists uniquely on t 0 in R + 5 with probability equal to one.
Proof. 
Since (1) fulfills the local Lipschitz criterion, in this case for:
( S 1 ( 0 ) , H ( 0 ) , S 2 ( 0 ) , I ( 0 ) , R ( 0 ) ) R + 5 ,
we have ( S 1 , H , S 2 , I , R ) as a unique local solution on [ 0 , τ e ) in which τ e is the explosion time. Now, our aim is to show τ e = for the global solution of (1). Let k 0 0 be very large so that ( S 1 ( 0 ) , S 2 ( 0 ) , I ( 0 ) , H ( 0 ) , R ( 0 ) ) lies in [ 1 k 0 , k 0 ] . For k k 0 , define:
τ e = inf t [ 0 , τ e ) : min { S 1 ( t ) , H ( t ) , S 2 ( t ) , I ( t ) , R ( t ) } 1 k or max { S 1 ( t ) , S 2 ( t ) , H ( t ) , I ( t ) , R ( t ) } k .
By assuming as the empty set, inf = . Since τ k is increasing for k , set τ = l i m k τ k . Then, we have τ τ e a.s.
Next, we aim to show τ = a.s. If this claim is not true, then T > 0 and ϵ ( 0 , 1 ) such that P { τ T } > ϵ . As a result, we have k 1 k 0 s.t. P { τ T } ϵ , k k 1 and for t τ k .
Now, consider the C 2 -function V ˙ : R + 5 R ˙ + by the following formulation:
V ˙ ( S 1 , S 1 , I 1 , H , R ) = ( S 1 c c ln S 1 c ) + ( S 2 1 ln S 2 ) + ( I 1 ln I ) + ( H 1 ln H ) + ( R 1 ln R ) .
By Ito’s formula, we have:
d V ˙ ( S 1 , S 2 , I 1 , H , R ) = ( 1 c S 1 ) d S 1 + 1 2 S 1 2 ( d S 1 ) 2 + ( 1 1 S 2 ) d S 2 + 1 2 S 2 2 ( d S 2 ) 2 + ( 1 1 I ) d I + 1 2 I 2 ( d I ) 2 + ( 1 1 H ) d H + ( 1 1 R ) d R = L ˙ V ˙ + σ 1 ( I S 1 ) d B 1 + σ 2 ( I S 2 ) d B 2
where L ˙ V ˙ : R + 5 R ˙ + acts as:
L ˙ V ˙ = 1 c S 1 Λ 1 β 1 I S 1 α S 1 + 1 2 σ 1 2 I 2 + 1 1 S 2 Λ 2 β 2 I S 2 α S 2 + 1 2 σ 2 2 I 2 + 1 1 I I β 1 S 1 + β 2 S 2 ( μ + α ) I ) + 1 2 σ 1 2 S 1 2 + 1 2 σ 2 2 S 2 2 + 1 1 H μ I ( ρ + α ) H + 1 1 R ρ H α R = Λ 1 β 1 I S 1 α S 1 c Λ 1 S 1 + c β 1 I + c α + 1 2 σ 1 2 I 2 + λ 2 β 2 I S 2 α S 2 λ 2 S 2 + β 2 I + α + 1 2 σ 2 2 I 2 + I ( β 1 S 1 + β 2 S 2 ) ( μ + α ) I ( β 1 S 1 + β 2 S 2 ) + ( μ + α ) + 1 2 σ 1 2 S 1 2 + 1 2 σ 2 2 S 2 2 + μ I ( ρ + α ) H μ I H + ( ρ + α ) + ρ H α R ρ H R + α Λ 1 + Λ 2 + β 2 I + 3 α + μ + ρ + 1 2 σ 1 2 I 2 + 1 2 σ 2 2 I 2 + 1 2 σ 1 2 S 1 2 + 1 2 σ 2 2 S 2 2 = Λ 1 + Λ 2 + β 2 I + 3 α + μ + ρ + σ 1 2 I 2 + σ 2 2 I 2 + σ 1 2 S 1 2 + σ 2 2 S 2 2 2 .
Set the constant : K = Λ 1 + Λ 2 + β 2 I + 3 α + μ + ρ + σ 1 2 I 2 + σ 2 2 I 2 + σ 1 2 S 1 2 + σ 2 2 S 2 2 2 .
Then:
d V ˙ K d t + σ 1 ( I S 1 ) d B 1 + σ 2 ( I S 2 ) d B 2 .
Now, integrate both sides of (11) between 0 t o σ T * and calculate the expectation; then, we obtain the inequality:
E [ V ˙ ( S 1 ( σ T * ) , S 2 ( σ T * ) , I ( σ T * ) , H ( σ T * ) , R ( σ T * ) ) V ˙ ( 0 ) + E 0 σ T * K d s V ˙ ( 0 ) + K E ( σ T * ) V ˙ ( 0 ) + K T * ,
where:
V ˙ ( 0 ) = V ˙ ( S 1 ( 0 ) , S 2 ( 0 ) , I ( 0 ) , H ( 0 ) , R ( 0 ) ) .
We take Ω = { σ T * } , 1 , which implies for every P ( Ω ) ξ . We remark that for every ω Ω , there exist at least S 1 ( σ , ω ) , S 2 ( σ , ω ) , I ( σ , ω ) , H ( σ , ω ) , R ( σ , ω ) , equaling the number or 1 , such that:
V ˙ [ S 1 ( σ T * , ω ) , S 2 ( σ T * , ω ) , I ( σ T * , ω ) , H ( σ T * , ω ) , R ( σ T * , ω ) ] ( 1 c c l n ( c ) ) ( 1 l n ) .
From equation(12), we obtain the following relation:
V ˙ ( 0 ) + K T * [ 1 Ω ( S 1 ( σ T * ) , S 2 ( σ T * ) , I ( σ T * ) , H ( σ T * ) , R ( σ T * ) ] ξ ( 1 c c l n ( c ) ) ( 1 l n ) ,
in which 1 Ω stands for the indicator function. Furthermore, it is seen that when K , then:
> V ˙ ( 0 ) + K T * =
implies that σ = a.s., which means that the solution ( S 1 , S 2 , I , H , R ) is globally positive a.s. □
Theorem 1.
Let ( S 1 t , S 2 t , I t , H t , R t ) be the solution of the stochastic epidemic model (1). Then, for R 0 = β 1 Λ 1 + β 2 Λ 2 α   ( μ + α ) < 1 , we have:
lim t 1 t log ( I t ) < 0 , lim t H t = 0 , lim t 0 t S 2 ( t ) d t = Λ 2 ( μ + α ) β 2 Λ 2 + α ( μ + α ) , lim t 0 t S 1 ( t ) d t = Λ 1 ( μ + α ) β 1 Λ 1 + α ( μ + α ) , lim t 1 t 0 t R ( t ) d t = 0 .
Proof. 
With the help of (1), we have:
d log ( I ) = β 1 S 1 + β 2 S 2 ( μ + α ) d t + σ 1 S 1 d B 1 ( t ) + σ 2 S 2 d B 2 ( t ) .
This implies:
log ( I t ) log ( I 0 ) t = 1 t 0 t β 1 S 1 + β 2 S 2 ( μ + α ) d s + 1 t 0 t σ 1 S 1 d B 1 ( s ) + 1 t 0 t σ 2 S 2 d B 2 ( s ) β 1 Λ 1 α + β 2 Λ 2 α ( μ + α ) + Λ 1 α 1 t 0 t σ 1 d B 1 ( s ) + 1 t Λ 2 α 0 t σ 2 d B 2 ( s ) .
By the help of the large number approach in local martingales, studied in [47], we have:
lim t Λ 1 α 1 t 0 t σ 1 d B 1 ( s ) + 1 t Λ 2 α 0 t σ 2 d B 2 ( s ) 0 .
Thus, for R 0 < 1 , (16) yields:
log ( I t ) log ( I 0 ) t β 1 Λ 1 α + β 2 Λ 2 α ( μ + α ) = ( μ + α ) β 1 Λ 1 + β 2 Λ 2 α ( μ + α ) 1 = ( μ + α ) R 0 1 < 0 .
We define Ω 1 = { Ω : lim sup t I ( , t ) = 0 } . By (18), we have P ( Ω 1 ) = 1 . This proves that the infection is bounded and controlled with ξ 1 .
Next, with the help of (1), for the infection, we have:
d H ( , t ) = ( μ I ( , t ) ( ρ + α ) H ( , t ) ) d t .
This yields the following:
H ( , t ) μ ξ 1 ( α + ρ ) .
By the arbitrariness of the ξ 1 , we have H ( , t ) 0 . Furthermore, on the other side,
lim t H ( , t ) 0 .
Thus, there exists some constant, say ξ 2 , such that H ( t ) < ξ 2 .
Next for S 2 ( t ) , we have:
β 2 Λ 2 α + μ + α 0 t S 2 ( t ) d t = Λ 2 S 2 t S 2 0 t ,
and so:
0 t S 2 ( t ) d t = Λ 2 α + μ β 2 Λ 2 + α ( α + μ ) α + μ β 2 Λ 2 + α ( α + μ ) S 2 t S 2 0 t .
This implies:
lim t 0 t S 2 ( t ) d t = Λ 2 ( α + μ ) β 2 Λ 2 + α ( α + μ ) .
Similarly, for S 1 ( t ) , we have:
lim t 0 t S 1 ( t ) d t = Λ 1 ( α + μ ) β 1 Λ 1 + α ( α + μ ) ,
and for R ( t ) , we obtain:
lim t 0 t R ( t ) d t = ρ ξ 2 α .
However, due to the arbitrariness of the ξ 2 and t , we have:
lim t 1 t 0 t R ( t ) d t 0 .
This completes the proof. □
Theorem 2.
Let ( S 1 t , S 2 t , I t , H t , R t ) be the solution of the stochastic epidemic model (1), via the initial values given by ( S 1 t ( 0 ) , S 2 t ( 0 ) , I t ( 0 ) , H t ( 0 ) , R t ( 0 ) ) . Then, for R 0 = β 1 Λ 1 + β 2 Λ 2 α ( μ + α ) > 1 , we have:
lim t inf S 1 Λ 1 α , lim t inf S 2 Λ 2 α , lim t inf H μ ( Λ 1 + Λ 2 ) ( ρ + α ) ( α + μ ) , lim t inf I ( Λ 1 + Λ 2 ) ( α + μ ) , lim t inf R ρ μ ( Λ 1 + Λ 2 ) α ( ρ + α ) ( α + μ ) .
Proof. 
In view of the first three equations of the proposed stochastic system (1), we obtain:
d S 1 + d S 2 + d I = Λ 1 + Λ 2 α S 1 α S 2 ( μ + α ) I .
Integrating (19), we have:
S 1 ( 0 ) S 2 ( 0 ) + S 1 ( t ) + S 2 ( t ) + I ( t ) I ( 0 ) = ( Λ 1 + Λ 2 ) t α 0 t S 1 ( t ) d t α 0 t S 2 ( t ) d t ( μ + α ) 0 t I ( t ) d t .
This implies:
( μ + α ) t 0 t I ( t ) d t = ( Λ 1 + Λ 2 ) α t 0 t S 1 ( t ) d t α t 0 t S 2 ( t ) d t S 1 ( 0 ) S 2 ( 0 ) + S 1 ( t ) + S 2 ( t ) + I ( t ) I ( 0 ) t .
This further gives:
lim t 1 t 0 t I ( t ) d t = ( Λ 1 + Λ 2 ) ( μ + α ) lim t α t 0 t S 1 ( t ) d t lim t α t 0 t S 2 ( t ) d t lim t S 1 ( t ) S 1 ( 0 ) S 2 ( 0 ) + S 2 ( t ) + I ( t ) I ( 0 ) t ( μ + α ) .
Ultimately, we have:
lim t inf 1 t 0 t I ( t ) d t ( Λ 1 + Λ 2 ) ( μ + α ) > 0 .
Next, for S 1 ( t ) , we have:
d S 1 = ( Λ 1 β 1 I S 1 α S 1 ) d t σ 1 I S 1 d B 1 ( t ) .
Integrating (21), we have:
S 1 ( t ) S 1 ( 0 ) t = Λ 1 1 t β 1 0 t I S 1 d s + α 0 t S 1 d s + 0 t σ 1 I S 1 d B 1 ( s ) .
Thus:
lim t α t 0 t S 1 ( s ) d s = Λ 1 lim t 1 t S 1 ( t ) S 1 ( 0 ) t + β 1 0 I t S 1 d s + 0 t σ 1 I S 1 d B 1 ( s ) .
This further gives:
lim t inf 1 t 0 t S 1 ( s ) d s Λ 1 α .
Similarly, we have:
lim t inf 1 t 0 t S 2 ( s ) d s Λ 2 α .
For H , integrating the fourth equation in the system (1), we obtain:
H ( 0 ) + H ( t ) t = μ t 0 t I ( s ) d s ( ρ + α ) t 0 t H ( s ) d s .
This implies that:
( ρ + α ) t 0 t H ( s ) d s = μ t 0 t I ( s ) d s ( H ( 0 ) + H ( t ) t ) .
With the help of (20), we further have:
lim t 1 t 0 t H ( s ) d s 1 ( ρ + α ) μ ( Λ 1 + Λ 2 ) ( α + μ ) lim t H ( t ) H ( 0 ) t ,
and:
lim t inf 1 t 0 t H ( s ) d s μ ( Λ 1 + Λ 2 ) ( ρ + α ) ( α + μ ) .
Integrating the fifth equation of the system (1), we obtain:
R ( 0 ) + R ( t ) t = ρ t 0 t H ( s ) d s α 0 t R d s .
Therefore:
α 0 t R d s = ρ t 0 t H ( s ) d s R ( t ) R 0 t .
With the help of (22), we reach:
lim t inf 0 t R d s μ ρ ( Λ 1 + Λ 2 ) α ( ρ + α ) ( α + μ ) .
This ends the proof. □

3. Numerical Scheme and Simulations

In this section, we provide a numerical scheme for the suggested stochastic COVID-19 epidemic model (1). We have the parametric values Λ 1 = 0.2 , Λ 1 = 0.05 , α = 0.25 , μ = 0.1 , ρ = 0.3 , S 1 ( 0 ) = 0.45 , S 2 ( 0 ) = 0.15 , I ( 0 ) = 0.1 , H ( 0 ) = 0.45 , H ( 0 ) = 0.1 . For the numerical values, we refer the readers to [48,49]. For the numerical results, we discretize the suggested stochastic mathematical model (1). The following discretization is based on the work given in [50,51] where the notion of the numerical scheme and its stability were discussed.
S 1 k + 1 = S 1 k + Λ 1 β 1 I k S 1 k α S 1 k Δ t σ 1 I k S 1 k Δ t τ k σ 1 2 2 I k S 1 k ( τ k 2 1 ) Δ t , S 2 k + 1 = S 2 k + Λ 1 β 2 I S 2 α S 2 Δ t σ 2 I k S 2 k Δ t τ k σ 2 2 2 I k S 2 k ( τ k 2 1 ) Δ t , I k + 1 = I k + I k β 1 S 1 k + β 2 S 2 k ( μ + α ) I k Δ t + σ 1 S 1 k I k τ k Δ t + σ 1 2 2 I k S 1 k ( τ k 2 1 ) Δ t + σ 2 S 2 k I k τ k Δ t + σ 2 2 2 I k S 2 k ( τ k 2 1 ) Δ t , H k + 1 = H k + μ I k ( ρ + α ) H k Δ t , R k + 1 = R k + ρ H k α R k Δ t .
The graphical explanation of the model (1) is given onward. In Figure 1, a numerical solution of the model (1) is given for the zero noise, i.e., σ 1 = σ 2 = 0.0 . In Figure 2, we give a numerical solution to the aforementioned epidemic stochastic model (1) for the zero noise, i.e., σ 1 = σ 2 = 0.30 . For the zero noise, we have a numerical solution given by Figure 3 of the model (1), i.e., σ 1 = σ 2 = 0.40 . In Figure 4, we give a numerical solution to (1) for the zero noise, i.e., σ 1 = σ 2 = 0.50 . In Figure 5, another numerical solution is plotted to the model (1) for the zero noise, i.e., σ 1 = σ 2 = 0.60 . The class S 1 is compared for different values of σ 1 = σ 2 = 0.20 , 0.40 , 0.50 , 0.60 , 0.0 in Figure 6. Figure 7 represents a numerical analyses for the class S 2 keeping σ 1 = σ 2 = 0.20 , 0.40 , 0.50 , 0.60 , 0.0 . These analyses show that there is a vital role of the noise in the spread of COVID-19. Next, we give the graphical presentation for the role of β 1 and β 2 for fixed σ 1 and σ 2 .
We give the numerical results with the help of graphs for the analysis of the stochastic model (1) based on β 1 and β 2 . In Figure 8, we give the numerical presentation of S 1 , S 2 , I for the role of β 1 and β 2 such that β 1 = 0.25 , β 2 = 0.2 keeping σ 1 = σ 2 = 0.3 . In Figure 9, we have numerical data S 1 , S 2 , I for the role of β 1 and β 2 with β 1 = 0.25 , β 2 = 0.45 and σ 1 = σ 2 = 0.3 . In Figure 10, we have numerical values S 1 , S 2 , I for the role of β 1 and β 2 with β 1 = 0.30 , β 2 = 0.50 and σ 1 = σ 2 = 0.3 . In Figure 11, we have a numerical presentation of S 1 , S 2 , I for the role of β 1 and β 2 with β 1 = 0.35 ,   β 2 = 0.55 and σ 1 = σ 2 = 0.3 . In Figure 12, we have a numerical presentation of the role of β 1 and β 2 with β 1 = 0.25 , 0.3 , 0.35 , 0.4 , 0.6 , and β 2 = 0.45 , 0.5 , 0.55 , 0.6 , 0.2 and σ 1 = σ 2 = 0.3 for the comparative analysis of the S 1 class. In Figure 13, we have a numerical presentation of the role of β 1 and β 2 with β 1 = 0.25 , 0.3 , 0.35 , 0.4 , 0.6 , and β 2 = 0.45 , 0.5 , 0.55 , 0.6 , 0.2 and σ 1 = σ 2 = 0.3 for the comparative analysis of the S 2 class.

4. Conclusions

In this article, we designed a new stochastic mathematical model for the spread of COVID-19 along with the white noise in the environment. The model was investigated mathematically in terms of the existence of the solution, the persistence of the disease, and the extinction of the disease. After the mathematical analysis of the model, we gave the numerical scheme for the computational analysis of the model, and the scheme was then tested for the available parametric and initial data in the literature. The computational results presented more realistic data for us. The numerical analysis was presented with the help of several graphs with respect to different noises σ 1 , σ 2 and the parameters β 1 , β 2 . These parametric values have important roles in the stochastic model (1). Our results showed that the disease dies out in the case of R 0 < 1 and persists if R 0 > 1 . In future works, we will aim to implement these criteria for different stochastic epidemic models of diseases. With the help of Figure 2, Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7, we observed that by increasing the values of σ 1 and σ 2 , an increase of the susceptible compartment occurs. Our presumed model is an S I R model where the susceptible class is categorized into two classes: S 1 , S 2 . In the future, the readers can follow the same technique to consider the epidemic models studied in [52,53,54,55] with a similar procedure as a continuation of the work and a further generalization of the study.

Author Contributions

Conceptualization, S.H., E.N.M. and H.K.; formal analysis, S.H., E.N.M., H.K., S.E., S.R., T.S. and N.P.; funding acquisition, N.P. and T.S.; methodology, S.H., E.N.M., H.K., S.E., S.R., T.S. and N.P.; software, S.H., H.K. and S.E.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by King Mongkut’s University of Technology North Bangkok, Contract Number KMUTNB-64-KNOW-43.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable for this article as no datasets were generated nor analyzed during the current study.

Acknowledgments

The fourth and fifth authors would like to thank Azarbaijan Shahid Madani University. Furthermore, the authors would like to thank the reviewers for their valuable and constructive comments, which improved the quality of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gorbalenya, A.E.; Baker, S.C.; Baric, R.S.; de Groot, R.J.; Drosten, C.; Gulyaeva, A.A.; Haagmans, B.L.; Lauber, C.; Leontovich, A.M.; Neuman, B.W. Coronaviridae study group of the international committee on taxonomy of viruses. The species severe acute respiratory syndrome-related coronavirus: Classifying 2019-nCoV and naming it SARS-CoV-2. Nat. Microbiol. 2020, 5, 536–544. [Google Scholar]
  2. Morrison, S.A.; Jurak, G.; Starc, G. Responding to a global pandemic: Republic of Slovenia on maintaining physical activity during self-isolation. Scand. J. Med. Sci. Sports 2020, 30, 1546–1548. [Google Scholar] [CrossRef] [PubMed]
  3. Kokudo, N.; Sugiyama, H. Call for international cooperation and collaboration to effectively tackle the COVID-19 pandemic. Global Health Med. 2020, 2, 60–62. [Google Scholar] [CrossRef]
  4. Omer, S.; Ali, S.; ud Din Babar, Z. Preventive measures and management of COVID-19 in pregnancy. Drugs Ther. Perspect. 2020, 36, 246–249. [Google Scholar] [CrossRef]
  5. Owusu, K.F.; Goufo, E.F.D.; Mugisha, S. Modelling intracellular delay and therapy interruptions within Ghanaian HIV population. Adv. Differ. Equ. 2020, 401, 2020. [Google Scholar] [CrossRef]
  6. Anderson, M.; Mckee, M.; Mossialos, E. Developing a sustainable exit strategy for COVID-19: Health, economic and public policy implications. J. Roy. Soc. Med. 2020, 113, 176–178. [Google Scholar] [CrossRef]
  7. Goufo, E.F.D.; Maritz, R. A note on Ebola’s outbreak and human migration dynamic. J. Hum. Ecol. 2015, 51, 257–263. [Google Scholar] [CrossRef]
  8. Djomegni, P.M.T.; Tekle, A.; Dawed, M.Y. Pre-exposure prophylaxis HIV/AIDS mathematical model with non classical isolation. Jpn. J. Ind. Appl. Math. 2020, 37, 781–801. [Google Scholar] [CrossRef]
  9. Khan, Z.A.; Jarad, F.; Khan, A.; Khan, H. Nonlinear discrete fractional sum inequalities related to the theory of discrete fractional calculus with applications. Adv. Differ. Equ. 2021, 100, 2021. [Google Scholar] [CrossRef]
  10. Khan, A.; Alshehri, H.M.; Abdeljawad, T.; Al-Mdallal, Q.M.; Khan, H. Stability analysis of fractional nabla difference COVID-10 model. Res. Phys. 2021, 22, 103888. [Google Scholar]
  11. Shah, A.; Khan, R.A.; Khan, A.; Khan, H.; Gomez-Aguilar, J.F. Investigation of a system of nonlinear fractional order hybrid differential equations under usual boundary conditions for existence of solution. Math. Methods Appl. Sci. 2021, 44, 1628–1638. [Google Scholar] [CrossRef]
  12. Shah, K.; Alqudah, M.A.; Jarad, F.; Abdeljawad, T. Semi-analytical study of pine wilt disease model with convex rate under Caputo-Febrizio fractional order derivative. Chaos Solitons Fractals 2020, 135, 109754. [Google Scholar] [CrossRef]
  13. Boutiara, A.; Etemad, S.; Hussain, A.; Rezapour, S. The generalized U-H and U-H stability and existence analysis of a coupled hybrid system of integro-differential IVPs involving φ-Caputo fractional operators. Adv. Differ. Equ. 2021, 95, 2021. [Google Scholar] [CrossRef]
  14. Amara, A.; Etemad, S.; Rezapour, S. Topological degree theory and Caputo-Hadamard fractional boundary value problems. Adv. Differ. Equ. 2020, 369, 2020. [Google Scholar] [CrossRef]
  15. Sulaiman, I.M.; Mamat, M.; Umar, A.O.; kamfa, K.; Madi, E.N. Some three-term conjugate gradient algorithms with descent condition for unconstrained optimization models. J. Adv. Res. Dyn. Control Syst. 2020, 12, 2494–2501. [Google Scholar]
  16. Madi, E.N.; Naim, S.; Yaafar, A.; Yaakob, A.M.; Yusoff, B. Agreement matrix based on fuzzy decision-making to rank ship Berthing criteria. Int. J. Eng. Trends Tech. 2020, 68, 31–36. [Google Scholar] [CrossRef]
  17. Etemad, S.; Rezapour, S.; Samei, M.E. α-ψ-contractions and solutions of a q-fractional differential inclusion with three-point boundary value conditions via computational results. Adv. Differ. Equ. 2020, 218, 2020. [Google Scholar] [CrossRef]
  18. Zhu, L.; Hu, H. A stochastic SIR epidemic model with density dependent birth rate. Adv. Differ. Equ. 2015, 330, 2015. [Google Scholar] [CrossRef] [Green Version]
  19. Mohammadi, H.; Kumar, S.; Rezapour, S.; Etemad, S. A theoretical study of the Caputo-Fabrizio fractional modeling for hearing loss due to Mumps virus with optimal control. Chaos Solitons Fractals 2021, 144, 110668. [Google Scholar] [CrossRef]
  20. Rezapour, S.; Ntouyas, S.K.; Iqbal, M.Q.; Hussain, A.; Etemad, S.; Tariboon, J. An analytical survey on the solutions of the generalized double-order ϕ-integrodifferential equation. J. Funct. Spaces 2021, 6667757, 2021. [Google Scholar]
  21. Rezapour, S.; Etemad, S.; Mohammadi, H. A mathematical analysis of a system of Caputo-Fabrizio fractional differential equations for the anthrax disease model in animals. Adv. Differ. Equ. 2020, 481, 2020. [Google Scholar] [CrossRef]
  22. Kongson, J.; Thaiprayoon, C.; Sudsutad, W. Analysis of a fractional model for HIV CD4+ T-cell with treatment under generalized Caputo fractional derivative. AIMS Math. 2021, 6, 7285–7304. [Google Scholar] [CrossRef]
  23. Matar, M.M.; Abbas, M.I.; Alzabut, J.; Kaabar, M.K.A.; Etemad, S.; Rezapour, S. Investigation of the p-Laplacian nonperiodic nonlinear boundary value problem via generalized Caputo fractional derivatives. Adv. Differ. Equ. 2021, 68, 2021. [Google Scholar] [CrossRef]
  24. Baleanu, D.; Etemad, S.; Rezapour, S. On a fractional hybrid integro-differential equation with mixed hybrid integral boundary value conditions by using three operators. Alex. Eng. J. 2020, 59, 3019–3027. [Google Scholar] [CrossRef]
  25. Brown, R. XXVII. A brief account of microscopical observations made in the months of June, July and August 1827, on the particles contained in the pollen of plants; and on the general existence of active molecules in organic and inorganic bodies. Philos. Mag. 1828, 4, 161–173. [Google Scholar] [CrossRef] [Green Version]
  26. Doob, J.L. Stochastic Processes; Wiley: New York, NY, USA, 1991; 664p. [Google Scholar]
  27. Williams, D. Probability with Martingales; Cambridge University Press: Cambridge, UK, 1991. [Google Scholar]
  28. Rogers, L.C.G.; Williams, D. Diffusions, Markov Processes and Martingales: Volume 1, Foundations, 2nd ed.; John Wiley and Sons, Ltd.: Chichester, UK, 1994. [Google Scholar]
  29. Burdzy, K. Brownian Motion and its Applications to Mathematical Analysis; Springer: Cham, Switzerland, 2014. [Google Scholar]
  30. Ross, R. An application of the theory of probabilities to the study of a priori pathometry, Part I. Proc. R. Soc. Lond. Ser. A 1916, 92, 204–230. [Google Scholar]
  31. Ross, R.; Hudson, H.P. An application of the theory of probabilities to the study of a priori pathometry, Part II. Proc. R. Soc. Lond. Ser. A 1917, 93, 212–225. [Google Scholar]
  32. Ross, R.; Hudson, H.P. An application of the theory of probabilities to the study of a priori pathometry, Part III. Proc. R. Soc. Lond. Ser. A 1917, 93, 225–240. [Google Scholar]
  33. Kermack, W.O.; McKendrick, A.G. A contribution to the mathematical theory of epidemics. Proc. R. Soc. Lond. Ser. A 1927, 115, 700–721. [Google Scholar]
  34. Harko, T.; Lobo, F.S.N.; Mak, M.K. Exact analytical solutions of the susceptible-infected-recovered (SIR) epidemic model and of the SIR model with equal death and birth rates. Appl. Math. Comput. 2014, 236, 184–194. [Google Scholar] [CrossRef] [Green Version]
  35. Kroger, M.; Schlickeiser, R. Analytical solution of the SIR-model for the temporal evolution of epidemics. Part A: Time-independent reproduction factor. J. Phys. A Math. Theor. 2020, 53, 505601. [Google Scholar] [CrossRef]
  36. Schlickeiser, R.; Kroger, M. Analytical solution of the SIR-model for the temporal evolution of epidemics: Part B. Semi-time case. J. Phys. A Math. Theor. 2021, 54, 175601. [Google Scholar] [CrossRef]
  37. Yang, W.; Zhang, D.; Peng, L.; Zhuge, C.; Hong, L. Rational evaluation of various epidemic models based on the COVID-19 data of China. Epidemics 2021, 37, 100501. [Google Scholar] [CrossRef] [PubMed]
  38. Ajbar, A.; Alqahtani, R.T.; Boumaza, M. Dynamics of an SIR-based COVID-19 model with linear incidence rate, nonlinear removal rate, and public awareness. Front. Phys. 2021, 9, 13. [Google Scholar] [CrossRef]
  39. Alanazi, S.A.; Kamruzzaman, M.M.; Alruwaili, M.; Alshammari, N.; Alqahtani, S.A.; Karime, A. Measuring and preventing COVID-19 using the SIR model and machine learning in smart health care. J. Healthc. Eng. 2020, 2020, 8857346. [Google Scholar] [CrossRef] [PubMed]
  40. Moein, S.; Nickaeen, N.; Roointan, A.; Borhani, N.; Heidary, Z.; Javanmard, S.H.; Ghaisari, J.; Gheisari, Y. Inefficiency of SIR models in forecasting COVID-19 epidemic: A case study of Isfahan. Sci. Rep. 2021, 11, 4725. [Google Scholar] [CrossRef]
  41. Sweilam, N.H.; Al-Mekhlafi, S.M.; Baleanu, D. A hybrid stochastic fractional order Coronavirus (2019-nCov) mathematical model. Chaos Solitons Fractals 2021, 145, 110762. [Google Scholar] [CrossRef]
  42. Babaei, A.; Jafari, H.; Banihashemi, S.; Ahmadi, M. Mathematical analysis of a stochastic model for spread of Coronavirus. Chaos Solitons Fractals 2021, 145, 110788. [Google Scholar] [CrossRef]
  43. Danane, J.; Allali, K.; Hammouch, Z.; Nisar, K.S. Mathematical analysis and simulation of a stochastic COVID-19 Levy jump model with isolation strategy. Res. Phys. 2021, 23, 103994. [Google Scholar] [CrossRef]
  44. Moualkia, S.; Xu, Y. On the existence and uniqueness of solutions for multidimensional fractional stochastic differential equations with variable order. Mathematics 2021, 9, 2106. [Google Scholar] [CrossRef]
  45. Sakthivel, R.; Revathi, P.; Ren, Y. Existence of solutions for nonlinear fractional stochastic differential equations. Nonlilear Anal. 2013, 81, 70–86. [Google Scholar] [CrossRef]
  46. Rafiq, M.; Ali, J.; Riaz, M.B.; Awrejcewicz, J. Numerical analysis of a bi-modal COVID-19 sitr model. Alex. Eng. J. 2022, 61, 227–235. [Google Scholar] [CrossRef]
  47. Mao, X. Stochastic Differential Equations and Applications, 2nd ed.; Woodhead Publishing: Chichester, UK, 2008. [Google Scholar]
  48. Sanchez, Y.G.; Sabir, Z.; Guirao, J.L.G. Design of a nonlinear SITR fractal model based on the dynamics of a novel Coronavirus (COVID-19). Fractals 2020, 28, 2040026. [Google Scholar] [CrossRef]
  49. Umar, M.; Sabir, Z.; Raja, M.A.Z.; Shoaib, M.; Gupta, M.; Sanchez, Y.G. A stochastic intelligent computing with neuro-evolution heuristics for nonlinear SITR system of novel COVID-19 dynamics. Symmetry 2020, 12, 1628. [Google Scholar] [CrossRef]
  50. Higham, D.J. An algorithmic introduction to numerical simulation of stochastic differential equations. SIAM Rev. 2001, 43, 525–546. [Google Scholar] [CrossRef]
  51. Zhou, Y.; Zhang, W.; Yuan, S. Survival and stationary distribution of a SIR epidemic model with stochastic perturbations. Appl. Math. Comput. 2014, 244, 118–131. [Google Scholar] [CrossRef]
  52. Kabir, K.M.A.; Risa, T.; Tanimoto, J. Prosocial behavior of wearing a mask during an epidemic: An evolutionary explanation. Sci. Rep. 2021, 11, 12621. [Google Scholar] [CrossRef]
  53. Chowdhury, A.; Kabir, K.M.A.; Tanimoto, J. How quarantine and social distancing policy can suppress the outbreak of novel coronavirus in developing or under poverty level countries: A mathematical and statistical analysis. Biom. Biostat. Int. J. 2021, 10, 145–152. [Google Scholar]
  54. Tanimoto, J. Sociophysics Approach to Epidemics; Springer: Singapore, 2021; Volume 23. [Google Scholar]
  55. Madi, E.N.; Yusoff, B. Modelling perceptive-based information (words) for decision support system. Int. J. Recent Technol. Eng. 2019, 7, 665–671. [Google Scholar]
Figure 1. Numerical solution of the model (1) for σ 1 = σ 2 = 0.0 .
Figure 1. Numerical solution of the model (1) for σ 1 = σ 2 = 0.0 .
Mathematics 09 03122 g001
Figure 2. Numerical solution of the model (1) for σ 1 = σ 2 = 0.3 .
Figure 2. Numerical solution of the model (1) for σ 1 = σ 2 = 0.3 .
Mathematics 09 03122 g002
Figure 3. Numerical solution of the model (1) for σ 1 = σ 2 = 0.4 .
Figure 3. Numerical solution of the model (1) for σ 1 = σ 2 = 0.4 .
Mathematics 09 03122 g003
Figure 4. Numerical solution of the model (1) for σ 1 = σ 2 = 0.5 .
Figure 4. Numerical solution of the model (1) for σ 1 = σ 2 = 0.5 .
Mathematics 09 03122 g004
Figure 5. Numerical solution of the model (1) for σ 1 = σ 2 = 0.6 .
Figure 5. Numerical solution of the model (1) for σ 1 = σ 2 = 0.6 .
Mathematics 09 03122 g005
Figure 6. Comparative analysis of S 1 ( t ) for σ 1 = σ 2 = 0.2 , 0.4 , 0.5 , 0.6 , 0.0 .
Figure 6. Comparative analysis of S 1 ( t ) for σ 1 = σ 2 = 0.2 , 0.4 , 0.5 , 0.6 , 0.0 .
Mathematics 09 03122 g006
Figure 7. Comparative analysis of S 2 ( t ) for σ 1 = σ 2 = 0.2 , 0.4 , 0.5 , 0.6 , 0.0 .
Figure 7. Comparative analysis of S 2 ( t ) for σ 1 = σ 2 = 0.2 , 0.4 , 0.5 , 0.6 , 0.0 .
Mathematics 09 03122 g007
Figure 8. Comparative analysis of S 1 ( t ) , S 2 ( t ) , I ( t ) for β 1 = 0.25 , β 2 = 0.2 .
Figure 8. Comparative analysis of S 1 ( t ) , S 2 ( t ) , I ( t ) for β 1 = 0.25 , β 2 = 0.2 .
Mathematics 09 03122 g008
Figure 9. Comparative analysis of S 1 ( t ) , S 2 ( t ) , I ( t ) for β 1 = 0.25 , β 2 = 0.45 .
Figure 9. Comparative analysis of S 1 ( t ) , S 2 ( t ) , I ( t ) for β 1 = 0.25 , β 2 = 0.45 .
Mathematics 09 03122 g009
Figure 10. Comparative analysis of S 1 ( t ) , S 2 ( t ) , I ( t ) for β 1 = 0.30 , β 2 = 0.50 .
Figure 10. Comparative analysis of S 1 ( t ) , S 2 ( t ) , I ( t ) for β 1 = 0.30 , β 2 = 0.50 .
Mathematics 09 03122 g010
Figure 11. Comparative analysis of S 1 ( t ) , S 2 ( t ) , I ( t ) for β 1 = 0.35 , β 2 = 0.55 .
Figure 11. Comparative analysis of S 1 ( t ) , S 2 ( t ) , I ( t ) for β 1 = 0.35 , β 2 = 0.55 .
Mathematics 09 03122 g011
Figure 12. Comparative analysis of S 1 ( t ) .
Figure 12. Comparative analysis of S 1 ( t ) .
Mathematics 09 03122 g012
Figure 13. Comparative analysis of S 2 ( t ) .
Figure 13. Comparative analysis of S 2 ( t ) .
Mathematics 09 03122 g013
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hussain, S.; Madi, E.N.; Khan, H.; Etemad, S.; Rezapour, S.; Sitthiwirattham, T.; Patanarapeelert, N. Investigation of the Stochastic Modeling of COVID-19 with Environmental Noise from the Analytical and Numerical Point of View. Mathematics 2021, 9, 3122. https://doi.org/10.3390/math9233122

AMA Style

Hussain S, Madi EN, Khan H, Etemad S, Rezapour S, Sitthiwirattham T, Patanarapeelert N. Investigation of the Stochastic Modeling of COVID-19 with Environmental Noise from the Analytical and Numerical Point of View. Mathematics. 2021; 9(23):3122. https://doi.org/10.3390/math9233122

Chicago/Turabian Style

Hussain, Shah, Elissa Nadia Madi, Hasib Khan, Sina Etemad, Shahram Rezapour, Thanin Sitthiwirattham, and Nichaphat Patanarapeelert. 2021. "Investigation of the Stochastic Modeling of COVID-19 with Environmental Noise from the Analytical and Numerical Point of View" Mathematics 9, no. 23: 3122. https://doi.org/10.3390/math9233122

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