
Stochastic differential equation models are important and provide more valuable outputs to examine the dynamics of SARS-CoV-2 virus transmission than traditional models. SARS-CoV-2 virus transmission is a contagious respiratory disease that produces asymptomatically and symptomatically infected individuals who are susceptible to multiple infections. This work was purposed to introduce an epidemiological model to represent the temporal dynamics of SARS-CoV-2 virus transmission through the use of stochastic differential equations. First, we formulated the model and derived the well-posedness to show that the proposed epidemiological problem is biologically and mathematically feasible. We then calculated the stochastic reproductive parameters for the proposed stochastic epidemiological model and analyzed the model extinction and persistence. Using the stochastic reproductive parameters, we derived the condition for disease extinction and persistence. Applying these conditions, we have performed large-scale numerical simulations to visualize the asymptotic analysis of the model and show the effectiveness of the results derived.
Citation: Tahir Khan, Fathalla A. Rihan, Muhammad Bilal Riaz, Mohamed Altanji, Abdullah A. Zaagan, Hijaz Ahmad. Stochastic epidemic model for the dynamics of novel coronavirus transmission[J]. AIMS Mathematics, 2024, 9(5): 12433-12457. doi: 10.3934/math.2024608
[1] | Roshan Ara, Saeed Ahmad, Zareen A. Khan, Mostafa Zahri . Threshold dynamics of stochastic cholera epidemic model with direct transmission. AIMS Mathematics, 2023, 8(11): 26863-26881. doi: 10.3934/math.20231375 |
[2] | Shah Hussain, Naveed Iqbal, Elissa Nadia Madi, Thoraya N. Alharthi, Ilyas Khan . Vaccination strategies in a stochastic SIVR epidemic model. AIMS Mathematics, 2025, 10(2): 4441-4456. doi: 10.3934/math.2025204 |
[3] | Yue Wu, Shenglong Chen, Ge Zhang, Zhiming Li . Dynamic analysis of a stochastic vector-borne model with direct transmission and media coverage. AIMS Mathematics, 2024, 9(4): 9128-9151. doi: 10.3934/math.2024444 |
[4] | Yajun Song, Ruyue Hu, Yifan Wu, Xiaohui Ai . Analysis of a stochastic two-species Schoener's competitive model with Lévy jumps and Ornstein–Uhlenbeck process. AIMS Mathematics, 2024, 9(5): 12239-12258. doi: 10.3934/math.2024598 |
[5] | Zhengwen Yin, Yuanshun Tan . Threshold dynamics of stochastic SIRSW infectious disease model with multiparameter perturbation. AIMS Mathematics, 2024, 9(12): 33467-33492. doi: 10.3934/math.20241597 |
[6] | Shah Hussain, Naveed Iqbal, Elissa Nadia Madi, Mohsen Bakouri, Ilyas Khan, Wei Sin Koh . On the stochastic modeling and forecasting of the SVIR epidemic dynamic model under environmental white noise. AIMS Mathematics, 2025, 10(2): 3983-3999. doi: 10.3934/math.2025186 |
[7] | Tingting Wang, Shulin Sun . Dynamics of a stochastic epidemic model with quarantine and non-monotone incidence. AIMS Mathematics, 2023, 8(6): 13241-13256. doi: 10.3934/math.2023669 |
[8] | Chuanfu Chai, Yuanfu Shao, Yaping Wang . Analysis of a Holling-type IV stochastic prey-predator system with anti-predatory behavior and Lévy noise. AIMS Mathematics, 2023, 8(9): 21033-21054. doi: 10.3934/math.20231071 |
[9] | Hassan Ranjbar, Leila Torkzadeh, Dumitru Baleanu, Kazem Nouri . Simulating systems of Itô SDEs with split-step (α,β)-Milstein scheme. AIMS Mathematics, 2023, 8(2): 2576-2590. doi: 10.3934/math.2023133 |
[10] | Jinji Du, Chuangliang Qin, Yuanxian Hui . Optimal control and analysis of a stochastic SEIR epidemic model with nonlinear incidence and treatment. AIMS Mathematics, 2024, 9(12): 33532-33550. doi: 10.3934/math.20241600 |
Stochastic differential equation models are important and provide more valuable outputs to examine the dynamics of SARS-CoV-2 virus transmission than traditional models. SARS-CoV-2 virus transmission is a contagious respiratory disease that produces asymptomatically and symptomatically infected individuals who are susceptible to multiple infections. This work was purposed to introduce an epidemiological model to represent the temporal dynamics of SARS-CoV-2 virus transmission through the use of stochastic differential equations. First, we formulated the model and derived the well-posedness to show that the proposed epidemiological problem is biologically and mathematically feasible. We then calculated the stochastic reproductive parameters for the proposed stochastic epidemiological model and analyzed the model extinction and persistence. Using the stochastic reproductive parameters, we derived the condition for disease extinction and persistence. Applying these conditions, we have performed large-scale numerical simulations to visualize the asymptotic analysis of the model and show the effectiveness of the results derived.
The emergence of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in late 2019 has caused unprecedented problems in health around the globe. The infection of the novel coronavirus has led to a devastating loss of lives and the production of a substantial burden on healthcare systems around the worldwide. Due to its urgency, the scientific community wants to comprehend the complexity of the disease and has hence investigated the spreading of SARS-CoV-2 not only by symptomatic individuals, but also by those who exhibited no apparent signs of illness but still transmitted the disease to others, i.e., asymptomatic individuals [1]. Symptomatic individuals experience cough, fever, and other respiratory symptoms which are recognized as the primary sources of disease transmission. The role of asymptomatic individuals has been gaining increasing amounts of attention because these individuals with no symptoms can transmit the disease to others making it difficult to control the spread of the novel coronavirus. Scholars have previously noted that 47% and 38% of newly infected cases are reported to be due to the contact with asymptomatic and symptomatic individuals, respectively [2]. Therefore, the existence of a small number of asymptomatic populations can produce a major disaster. The identification and understanding of the dynamics of both asymptomatic and symptomatic individuals can help to control the novel infection with the aid of effective control mechanisms.
Modeling the epidemiological process is a useful tool and has been vastly utilized to test various theories related to the dynamics and control of communicable diseases [3,4,5,6]. Different researchers have used these models to investigate the temporal dynamics of infectious diseases [7,8,9,10,11,12,13]. Dynamical systems representing the SARS-CoV-2 virus transmission are also abundant in literature, and numerous epidemiological models have been frequently exercised to investigate the predictability of the disease spread; see for instance [14,15,16,17]. Indeed, the models reported to represent the temporal dynamics of the novel coronavirus yielded interesting results; however, the authors mostly used a deterministic approach to model the interacting components of the novel model of coronavirus spread. Disease transmission is different everywhere due to the environment of the area and the immune system of the individual as well as the death rate and implementation of the vaccination program varying from place to place. Thus, due to the complex nature of SARS-CoV-2, it is observed that the novel disease has a stochastic nature. Recently, a model utilizing a stochastic approach has been studied to investigate the dynamics of coronavirus 2, as detailed in [18], where the authors did not consider various important factors that greatly influence the disease dynamics, e.g., the classification of infected individuals as a symptomatic and asymptomatic individuals has not been considered. However, models with an appropriate structure provide accurate dynamics and forecast the long-term spread for the disease which is very helpful for public health officials. Thus, we have formulated an epidemiological model by taking the extended version of the model reported in [18] to represent the accurate dynamics of the novel coronavirus and realize the following contributions:
● Infected individuals are included as asymptomatic and symptomatic individuals because they are very important in the pandemic trend of the SARS-CoV-2 virus transmission. In fact, the existence of asymptomatic individuals with a small amount of the virus can create a major outbreak because without experiencing any symptoms, they produce environmental reservoirs and transmit the disease to others.
● We have investigated the multi-infection pattern of the novel of coronavirus by applying various phases of asymptomatic and symptomatic, and we show that successful interaction between the susceptible population and infected individuals either leads to symptomatic or asymptomatic individuals; therefore, have applied probability-based transmission while applying the outflow of the susceptible population to either infected group of individuals.
● The randomness is taken in each group of the model compartments because every parameter involved in the epidemic process of the novel coronavirus disease has a stochastic nature.
● We have developed an algorithm for the model simulations by using Milstein's higher order procedure to check the validity of the analytical work and represent the long-term with the aid of numerical simulations.
Thus, the goal of this study was to explore the dynamics of novel coronavirus transmission by including the above features and enhancing the model by taking an extended version of the model studied in [18]. To this end, we chose to classify the various infected sources of disease transmission as symptomatic, asymptomatic, and environmental reservoirs, and to take the randomness in every group of individuals by using different Brownian motions according to the disease characteristics. Keeping in view these hypotheses, we chose to first formulate the model and discuss its feasibility in the form of existence analysis and positivity. For this, we will use a combination of the Itô formula and the Lyapunov function. We then calculate the various stochastic reproductive parameters to discuss the conditions for the persistence and extinction of the novel coronavirus. In addition, we discuss the detailed extinction and persistence analysis of the model and derived sufficient conditions for it in terms of the stochastic reproductive parameters. Finally, to interpret the obtained results based on graphical visualization, we chose to develop an algorithm by using the Milstein higher-order method and present the large-scale numerical simulations of the model.
To provide a comprehensive and systematic analysis of the paper, we present the organization of the paper as follows. Following the introduction and recent literature, we formulate the model with a complete description in Section 2. To show the feasibility of the epidemic problem, we provide some basic concepts and notations and discuss the model analysis in Section 3. Subsequently, to examine the model extinction and persistence, and to derive sufficient conditions for it, we investigate the extinction and persistence analysis in Section 4. In addition, to perform the large-scale numerical simulations of the model, we have developed the algorithm by using of Milstein's Higher Order method, as presented in Section 5. Particularly, the discretization of the model is reported in Section 5.1, while the numerical experiments are presented in Section 5.2. Lastly, we conclude the summary of our work with the main findings in Section 6.
In this part, we formulated the model for the dynamics of the novel coronavirus transmission keeping in view the characteristics of latent/asymptomatic and symptomatic individuals, as well as the environmental reservoirs. We assume that (Θ,QT,(Qt)t∈[0,T],P) is a space of filtered probability that includes a 5D-Brownian motion W:=(W(t))t≥0, where W(t):=(W1(t),W2(t),W3(t),W4(t),W5(t)). The novel disease transmits via environmental reservoirs and human interaction, differing everywhere. We, therefore, assume that groups of the model population will have a stochastic nature that is driven by different randomness sources i.e., Wi(t), i=1,…,5. The proposed epidemiological model incorporates the variations in all groups of individuals associated with various sources of information denoted by Q=(Qt)t≥0, where Qt:=σ(W(t)) represents the σ−algebra derived from W(t); thus, the fluctuations are considered to be taken in all groups of individuals in the model. Further, the groups of humans are categorized into susceptible, latent or asymptomatic, symptomatic, and recovered groups, which are respectively denoted by S(t), L(t), B(t), and R(t), while the concentration of environmental reservoirs is symbolized by C(t). In addition, we provided the flow diagram that represents the transfer mechanism for the model population in Figure 1.
Thus, the below dynamical system represents the novel coronavirus dynamics:
dS(t)={Π−γL(t)S(t)−βS(t)B(t)−ξS(t)C(t)−(α0+v)S(t)}dt+ϑ1S(t)dW1(t),dL(t)={q(γL(t)S(t)+βS(t)B(t)+ξS(t)C(t))−(α0+γ1)L(t)}dt+ϑ2L(t)dW2(t),dB(t)={(1−q)(γL(t)S(t)+βS(t)B(t)+ξS(t)C(t))+pγ1L(t)−(α1+γ2+α0)B(t)}dt+ϑ3B(t)dW3(t),dR(t)={γ1(1−p)L(t)+γ2B(t)+vS(t)−α0R(t)}dt+ϑ4R(t)dW4(t),dC(t)={η1L(t)+η2B(t)−αC(t)}dt+ϑ5W(t)dW5(t). | (2.1) |
In the model (2.1), the rate of births is denoted by Π, while the parameters γ, β, and ξ are used to represent the disease transmission from various sources. Vaccination for the novel coronavirus plays a significant role and it denoted by v, while α0 denotes the natural death rate. Moreover, the death rate attributable to the novel disease is described by α1, and p is the probability of direct recovery of latent/asymptomatic individuals, as some individuals directly recover without producing any symptoms. We also use q as the probability at which the susceptible individuals enter the asymptomatic/latent population, and γ1 and γ2 are the portions of individuals who recover from the contagious infection of the novel coronavirus disease. In addition, Wi(t) and ϑi, i=1,…,5 are the standard Brownian motions and white noise intensities, respectively, satisfying that Wi(0)=0 and ϑ2i>0 for all i=1,2,3,4,5.
In this part, we show the model feasibility to prove the well-posedness of the proposed epidemic model. Particularly, we discuss the existence analysis of the model and its positivity. However, before presenting the results, first, we briefly recall some basic concepts that are helpful for deriving our conclusion.
Lemma 3.1. [19] Let u=(u1,u2,…,uk) and v=(v1,v2,…,vk) represent two square-integrable adopted processes with k-dimensions. We assume the process Y=(y1,y2,…,yk) for n∈{1,2,3,…,k}; then, Yn is expressed as the following differential equation:
dYn=un(t)dt+vn(t)dW(t),Yn(0)∈R. | (3.1) |
If a function G from Rk to R is twice continuously differentiable, then
dG=k∑n=1∂G∂xndYn(t)+k∑n,m=112∂2G∂xn∂xm⟨dYn,dYm⟩, | (3.2) |
where ⟨,⟩ represents the quadratic variation; therefore, ⟨dW(t),dW(t)⟩=dt and ⟨dt,dW(t)⟩=⟨dW(t),dt⟩=⟨dt,dt⟩=0.
In addition, to examine the extinction and persistence with the aid of the proposed model, we first describe the definitions that will be used in the model extinction and persistence analysis. Let
⟨P(t)⟩=1t∫t0P(a)da. | (3.3) |
Then, the disease persists if lim(inf⟨L(t)⟩), lim(inf⟨B(t)⟩) and lim(inf⟨C(t)⟩) are positive whenever t increases without bound or equivalently, we can say that the epidemiological model (2.1) persists if
limt→∞inf∫t0L(a)da,limt→∞inf∫t0B(a)da,limt→∞inf∫t0C(a)da, | (3.4) |
are positive.
Moreover, we assume that RSa, RSb, RSc, and RSd are the components of the stochastic reproductive parameter, and they are respectively defined by
RSa=q{γ+β+ξ}Π{α0+v}{α0+γ1+ϑ222},RSb=(1−q){γ+β+ξ}Π{α0+v}{α0+α1+γ2+ϑ232},RSc=qγΠ{α0+v}{α0+γ1+ϑ222},RSd=(1−q)βΠ{α0+v}{α0+α1+γ2+ϑ232}. | (3.5) |
Then, the stochastic reproductive parameter for the epidemiological model (2.1) becomes RS0=RSa+RSb+RSc+RSd.
In order to show the model existence analysis, we use a combination of the Itô formula and Lyapunov theory. Define
Φ={(S,L,B,R,C)∈R5+S,R>0,A,B,C≥0,S+L+B+R+C≤1}. | (3.6) |
Thus, regarding the model existence the result will be shown in the following theorem.
Theorem 3.2. For the initial conditions (S(0),L(0),B(0),R(0),C(0)) contained in R5+, there exists a solution (S(t),L(t),B(t),R(t),C(t)) of the epidemiological model (2.1) which is unique and remains in R5+ such that
p{(S,L,B,R,C)∈Φ,∀t≥0}=1. |
Proof. Noted that for the model that is under consideration, the local Lipschitz continuity condition holds by following the methodology as adopted by Lei and Yang in [20]. Let X=(S,L,B,R,C) be the model solution for all t in [0,te) (where te is the explosion time), while X0 symbolizes the initial values (S(0),L(0),B(0),R(0),C(0)) for the model contained in R5+; then, the property of Lipschitz continuity implies that the solution of the model (2.1) is unique and local. To demonstrate that the solution is global, it is enough to show that te=∞. To proceed further, we assume a constant ϱ0≥0 that is sufficiently large apply 1ϱ0<X0<ϱ0 and define the stopping time as follows:
tϱ=inf{t∈[0,te):minX≤1ϱormaxX}. | (3.7) |
Note that infφ=∞, where φ represents the null set; also, since tρ varies with changing ϱ, that is if ϱ increases without bound, tϱ also increases. For t→∞, we set lim=t∞ and t∞=∞. Thus, we need to prove that te=∞ only. Let 0<ζ<1, T>0, and
P{t∞≤T}>ζ. | (3.8) |
For ϱ1≥ϱ0
P{tϱ≤T}≥ζ,for everyϱ≥ϱ1. | (3.9) |
Let H contained in C2 such that H:R5+→R+ be defined by
H=L−log(L)−1+B−log(B)−1+S−log(S)−1+R−1−log(R)+C−log(C)−1. | (3.10) |
Note that H≥0; thus, for ϱ0≤ϱ and T≥0, the implementation of the Itô formula yields
dH=UHdt+ϑ1(S−1)dW1+(L−1)ϑ2dW2+(B−1)ϑ3dW3+ϑ4(R−1)dW4+(C−1)ϑ5dW5, | (3.11) |
where
UH={1−1S}{Π−γLS−βSB−ξSC−(α0+v)S}+12ϑ21 +{1−1L}{q(γLS+βSB+ξSC)−(α0+γ1)L}+12ϑ22 +{1−1B}{(1−q)(γLS+βSB+ξSC)+pγ1L−(α0+α1+γ2)B}+12ϑ23 +{1−1R}{(1−p)γ1L+γ2B+vS−α0R}+12ϑ24 +{1−1C}{η1L+η2B−αC}+12ϑ25. | (3.12) |
Applying some algebraic manipulation and simplification, we arrive at
UH≤Π+4α0+v+γ1+α1+γ2+(γ+γ1+η1)L+(β+γ2+η2)B+ξC+α+12(ϑ21+ϑ22+ϑ23+ϑ24+ϑ25). | (3.13) |
Let K1=max{γ+γ1+η1,β+γ2+η2,ξ}; then, by following S+L+B+R+C≤1 from Eq (3.6), the last inequality appears as follows:
UH≤Π+K1+4α0+v+γ1+α1+γ2+12(ϑ21+ϑ22+ϑ23+ϑ24+ϑ25):=Φ. | (3.14) |
Plugging Eq (3.14) into Eq (3.11), we obtain
dH≤Φdt+(S−1)ϑ1dW1+(L−1)ϑ2dW2+(B−1)ϑ3dW3+(R−1)ϑ4dW4+(C−1)ϑ5dW5. |
Applying the integration on both sides gives
∫tϱ∧T0dH≤∫tϱ∧T0Φdt+∫tϱ∧T0(S−1)ϑ1dW1+∫tϱ∧T0(L−1)ϑ2dW2+∫tϱ∧T0(B−1)ϑ3dW3+∫tϱ∧T0(R−1)ϑ4dW4+∫tϱ∧T0(C−1)ϑ5dW5. | (3.15) |
The implementation of expectation yields that
E[H(S(tϱ∧T),L(tϱ∧T),B(tϱ∧T),R(tϱ∧T),C(tϱ∧T))]≤H(x(0))+E[∫tϱ∧T0Φdt], |
which implies that
E[H(S(tϱ∧T),L(tϱ∧T),B(tϱ∧T),R(tϱ∧T))]≤H(x(0))+TΦ. | (3.16) |
Let Φϱ=T≥tϱ for all ϱ≥ϱ1; then, P(Φϱ)≥ζ. For every ϵ∈Φϱ there exists at least one S(ϵ,tϱ), L(ϵ,tϱ), B(ϵ,tϱ), R(ϵ,tϱ) or C(ϵ,tϱ) that is either equal to ϱ or 1/ϱ, thus,
H(S(ϵ,tϱ),L(ϵ,tϱ),B(ϵ,tϱ),R(ϵ,tϱ),C(ϵ,tϱ))≥(1ϱ−1+logϱ)∧(ϱ−logϱ−1). | (3.17) |
Following Eqs (3.8) and (3.16), we may write
H(X0)+TΦ≥E[1Φϱ(ϵ)G(S(tϱ∧T),L(tϱ∧T),B(tϱ∧T),R(tϱ∧T))],=E[1Φϱ(ϵ)G(S(tϱ,T),L(tϱ,T),B(tϱ,T),R(tϱ,T))]≥E[1Φϱ(ϵ)(logϱ−1+1ϱ) ∧(−logϱ−1+ϱ)]=(logϱ−1+1ϱ)∧(−logϱ−1+ϱ)E[1Φϱ(ϵ)]. |
This gives that
H(X(0))+TΦ≥ζ(logϱ+1ϱ−1)∧(ϱ−logϱ−1). |
In the above equation 1Φϱ(ϵ) represents the indicator function for Φϱ(ϵ). For an increase in n without bound, ∞>H(X0)+ΦT=∞ yields a contradiction; so, t∞=∞.
Remark 1. The existence analysis for the considered model reveals that there exists a unique global solution X∈R5+ for any initial groups of population X(0)∈R5+.
Theorem 3.3. The proposed epidemiological model (2.1) possesses positive solutions for all initial groups of populations in R5+.
Proof. To begin the proof, we follow the work investigated in [21], and assume that [0,+∞) is the interval of solution. Since every differential equation in the model (2.1) corresponds to a geometric Brownian motion, we assume the 5 processes to be (ηS(t))t∈R+, (ηL(t))t∈R+, (ηB(t))t∈R+, (ηR(t))t∈R+ and (ηC(t))t∈R+, respectively denoted
dηS(t)=−{α0+v+γL(t)+βB(t)+ξC(t)}ηS(t)dt+ϑ1ηS(t)dW1(t),ηS(0)=1,dηL(t)={qγL(t)−(α0+γ1)}ηL(t)dt+ϑ2ηL(t)dW2(t),ηL(0)=1,dηB(t)={(1−q)βS(t)−(α0+α1+γ2)}ηB(t)dt+ϑ3ηB(t)dW3(t),ηB(0)=1,dηR(t)=−α0ηR(t)dt+ϑ4ηR(t)dW4(t),ηR(0)=1,dηC(t)=−α0ηC(t)dt+ϑ4ηC(t)dW5(t),ηC(0)=1. |
By applying the Itô formula to the function lnη., the solution of the above system can be easily derived, as follows:
ηS(t)=exp[∫t0(γL(x)+βB(x))dx−(α0+v+ϑ212)t+ϑ1W1(t)],ηL(t)=exp[qβ∫t0S(x)dx−(α0+γ1+ϑ222)t+ϑ2W2(t)],ηB(t)=exp[(1−q)β∫t0S(x)dx−(α0+α1+γ2+ϑ232)t+ϑ3W3(t)],ηR(t)=exp[−{α0+ϑ242}t+ϑ4W4(t)],ηC(t)=exp[−{α0+ϑ242}t+ϑ5W5(t)]. |
The solutions of the model process are derived based on the assumption that each solution can be expressed as a product of the stochastic process and associated geometric Brownian motion, i.e., S(t)=ηS(t)YS(t), and YS(t), respectively [21]. So, the stochastic integration based on the implementation of the integral to S(t)η−1S(t) gives the assertion for YS(t). Hence, the model solutions can be represented as follows:
S(t)=ηS(t){Π∫t0η−1S(x)dx+S(0)},L(t)=ηL(t){q∫t0βB(x)+ξX(x)S(x)η−1L(x)dx+L(0)},B(t)=ηB(t){B(0)+∫t0((1−q)(γS(x)+ξS(x)+pγ1)B(u))η−1B(x)dx},R(t)=ηR(t){R(0)+∫t0(γ1(1−p)L(x)+γ2C(x)+vS(x))η−1R(x)dx},C(t)=ηC(t){C(0)+∫t0(η1L(x)+η2C(x))η−1C(x)dx}. |
Clearly, all the solutions are positive because η. and η−1. are in exponential form while the initial data and all other epidemic parameters of the epidemiological model have non-negative values.
In this section, we discuss the disease extinction and persistence to derive the conditions for it, which will be in the form of an expression containing white noise intensities and model epidemic parameters. Thus, the extinction of the epidemiological model (2.1) is illustrated by giving the subsequent result.
Theorem 4.1. If
a. RSa<1 and RSb<1,
b. (1−q)(γ+β+ξ)(α0+α1+γ2)≥p(v+α0),
then the infection of the novel coronavirus will decay exponentially, that is
limt→∞suplogL(t)t≤{α0+γ1+12ϑ22}(RSa−1)<0, |
and
limt→∞suplogB(t)t≤{α0+α1+γ2+12ϑ23}(RSb−1)<0. |
In addition
limt→∞S(t)=Πv+α0,limt→∞L(t)=limt→∞B(t)=limt→∞C(t)=0,limt→∞R(t)=vΠα0(α0+v). |
Proof. To begin, we integrate both sides of the system (2.1), which gives the integral system as follows
∫t0dS(y)=Πt−∫t0{γL(y)+βB(y)+ξC(y)+(α0+v)}S(y)dy+ϑ1∫t0S(y)dW1(y),∫t0dL(x)=∫t0{q(γL(y)+βB(y)+ξC(y))S(y)−(α0+γ1)L(y)}dy+ϑ2∫t0L(y)dW2(y),∫t0dB(x)=∫t0{(1−q){γL(y)+βB(y)+ξC(y)}S(y)+pγ1L(y)−(α0+α1+γ2)C(y)}dy+ϑ3∫t0C(y)dW3(y),∫t0dR(y)=∫t0{(1−p)γ1L(y)+γ2C(y)+vS(y)−α0R(y)}dy+ϑ4∫t0R(y)dW4(y),∫t0dC(y)=∫t0{η1L(y)+η2C(y)−αC(y)}dy+ϑ5∫t0C(y)dW5(y). |
Considering the above integral system with some algebraic manipulation, we have
S(t)−S(0)t=Π−γ⟨L(t)S(t)⟩−β⟨B(t)S(t)⟩−ξ⟨C(t)s(t)⟩−(α0+v)⟨S(t)⟩+ϑ1t∫t0S(y)dW1(x),L(t)−L(0)t=q{γ⟨L(t)S(t)⟩+β⟨B(t)S(t)⟩+ξ⟨C(t)s(t)⟩}−(α0+γ1)⟨L(t)⟩+ϑ2t∫t0L(y)dW2(y),B(t)−B(0)t=(1−q){γ⟨L(t)S(t)⟩+β⟨B(t)S(t)⟩+ξ⟨C(t)s(t)⟩}+pγ1⟨L(t)⟩−(α0+α1+γ2)⟨B(t)⟩+ϑ3t∫t0C(y)dW3(y),R(t)−R(0)t=v⟨S(t)⟩+γ1(1−p)⟨L(t)⟩+γ2⟨C(t)⟩−α0⟨R(t)⟩+ϑ4t∫t0R(y)dW4(y),C(t)−C(0)t=η1⟨L(t)⟩+η2⟨C(t)⟩−α⟨C(t)⟩+ϑ5t∫t0C(y)dW5(y). |
The addition of the 1st three equations gives
S(t)−S(0)t+L(t)−L(0)t+B(t)−B(0)t= Π−(α0+v)⟨S(t)⟩−(α0+γ1)⟨L(t)⟩ +pγ1⟨L(t)⟩−(α0+α1+γ2)⟨C(t)⟩ +ϑ1t∫t0SdW1(y)+ϑ2t∫t0LdW2(y)+ϑ3t∫t0CdW3(y), | (4.1) |
which implies that
⟨S(t)⟩=Π(α0+v)−{α0+γ1(1−p)}(α0+v)⟨L(t)⟩−(α0+α1+γ2)(α0+v)⟨B(t)⟩+Υ(t), | (4.2) |
where
Υ(t)=1(α0+v){ϑ1t∫t0S(y)dW1(y)+ϑ2t∫t0L(y)dW2(y)+ϑ3t∫t0B(y)dW3(y)−(S(t)−S(0)t+L(t)−L(0)t+B(t)−B(0)t)}. | (4.3) |
Obviously Υ(t)=0 as t increases without bounds to ∞
limt→∞Υ(t)=0. | (4.4) |
By the application of the Itô lemma to the model (2.1), the second equation yields
dlogL(t)=q{γS(t)+βS(t)B(t)L(t)+ξS(t)C(t)L(t)}−(α0+γ1)−ϑ222+ϑ2dW2(t). | (4.5) |
Integrating the above equation and then dividing both sides by t, we get the following assertion:
1tlogL(t)|t0=q{γ⟨S(t)⟩+β⟨S(t)B(t)L(t)⟩+ξ⟨S(t)C(t)L(t)⟩}−(α0+γ1)−ϑ222+ϑ2W2(t)t. | (4.6) |
Since, ⟨S(t)B(t)L(t)⟩≤⟨S(t)B(t)⟩≤⟨S(t)⟩ and ⟨S(t)C(t)L(t)⟩≤⟨S(t)C(t)⟩≤⟨S(t)⟩, the above equation can be re-expressed follows
1tlogL(t)|t0≤q{(γ+β+ξ)⟨S(t)⟩}−(α0+γ1)−ϑ222+ϑ2W2(t)t. | (4.7) |
Plugging Eq (4.2) into last equation, we arrive at
1tlogL(t)|t0≤{α0+γ+ϑ222}{q(γ+β+ξ)Π(α0+v)(α0+γ1+ϑ222)−1}−q(γ+β+ξ){α0+(1−p)γ1}(α0+v)⟨L(t)⟩−q(γ+β+ξ)(α0+α1+γ2)(α0+v)⟨B(t)⟩ +q(γ+β+ξ)Υ(t)+ϑ2W2(t)t. | (4.8) |
Using the stochastic reproductive parameters and following the Strong Law of Large Numbers [22], the final inequality may take the form given by
limt→∞suplogL(t)t≤{α0+γ1+ϑ222}(RSa−1)<0a.s. | (4.9) |
Thus, limL(t)=0 as t increases without bound and RS1<1.
To proceed further by once again using the Itô lemma, the third equation of the epidemiological model (2.1) appears as follows:
dlogB(t)=(1−q){γL(t)S(t)+βB(t)S(t)+ξC(t)S(t)B(t)} +pγ1L(t)B(t)−(α0+α1+γ2)−ϑ232+ϑ3dW3(t), | (4.10) |
or equivalently, we can write
dlogB(t)≤{γS(t)L(t)+βB(t)S(t)+ξC(t)S(t)}(1−q)+pγ1L(t)−(α0+α1+γ2)−ϑ232+ϑ3dW3(t). |
In a similar fashion as above, once again using stochastic integration, we arrive at
1tlogB(t)|t0≤(1−q){γ⟨L(t)S(t)⟩+β⟨S(t)B(t)⟩+ξ⟨S(t)C(t)⟩} +pγ1⟨L(t)⟩−(α0+α1+γ2)−ϑ232+ϑ3W3(t)t. | (4.11) |
Using the same steps as performed for the previous case, the last inequality can be re-casted as follows:
1tlogB(t)|t0≤(1−q){γ⟨s(t)⟩+β⟨s(t)⟩+ξ⟨s(t)⟩} +pγ1⟨L(t)⟩−(α0+α1+γ2)−ϑ232+ϑ3W3(t)t. | (4.12) |
Plugging in the value of ⟨S(t)⟩, Eq (4.12) can be transformed to obtain
1tlogB(t)|t0≤{α0+α1+γ2+ϑ232}{(1−q)(γ+β+ξ)Π(α0+v){α0+α1+γ2+ϑ232}−1} −{(1−q)(γ+β+ξ){α0+α1+γ2}α0+v−pγ1}⟨L(t)⟩ −(1−q)(γ+β+ξ)(α0+α1+γ2)(α0+v)⟨B(t)⟩ +(1−q)(γ+β+ξ)Υ(t)+ϑ3W3(t)t. | (4.13) |
If (1−q)(γ+β+ξ)(α0+α1+γ2)≥pγ1(α0+v), then, through the use of the stochastic reproductive parameters and the implementation of the strong law of large numbers with little rearrangement and simplification, the last inequality is given by
limt→∞suplogB(t)t≤{α0+α1+γ2+ϑ232}(RSb−1)<0. | (4.14) |
Hence, limB(t)=0 as t approaches ∞ whenever RSb<1. In addition, we take the limiting system of the first, fourth, and fifth equations of the proposed model, whose solution yields that limS(t)=Π/(α0+v), limR(t)=vΠα0(v+α0) and limC(t)=0 if t increases without bound. Thus, the extinction is subject to the values of the stochastic reproductive parameters RSa and RSb, and if they are less than unity, the disease will go extinct.
Theorem 4.2. For the initial groups of the population (S(0),L(0),B(0),R(0),C(0)) in Φ, the solutions of the epidemiological model (2.1) satisfy
L2≤limt→∞inf⟨L(t)⟩≤sup⟨L(t)⟩≤L1, |
B2≤limt→∞inf⟨B(t)⟩≤sup⟨B(t)⟩≤B1, |
and
C2≤limt→∞inf⟨C(t)⟩≤sup⟨C(t)⟩≤C1, |
where
L1=(α0+v){α0+γ+ϑ222}(RSa−1)q(γ+β+ξ){α0+(1−p)γ1},L2=(v+α0){α0+γ1+ϑ222}(RSc−1)qγ{α0+γ1(1−p)},B1=(α0+v){α0+α1+γ2+ϑ232}(RSb−1)(1−q)(γ+β+ξ)(α0+α1+γ2),B2={α0+α1+γ2+ϑ232}(RSd−1),C1=η1αL1+η2αB1,C2=η1αL2+η2αB2. |
Proof. The direct use of Eq (4.8) ultimately gives that
1tlogL(t)|t0≤{α0+γ+ϑ222}{RSa−1} −q(γ+β+ξ){α0+(1−p)γ1}(α0+v)⟨L(t)⟩+q(γ+β+ξ)Υ(t)+ϑ2W2(t)t, | (4.15) |
which implies that
⟨L(t)⟩≤(α0+v){α0+γ+ϑ222}(RSa−1)q(γ+β+ξ){α0+(1−p)γ1} +(α0+v)q(γ+β+ξ){α0+(1−p)γ1}{q(γ+β+ξ)Υ(t)+ϑ2W2(t)t−1tlogL(t)|t0}. |
Using the sup property and applying lim as t approaches ∞, we derive the following:
limt→∞sup⟨L(t)⟩≤(α0+v){α0+γ+ϑ222}(RSa−1)q(γ+β+ξ){α0+(1−p)γ1}=L1. | (4.16) |
Similar to Eq (4.6), we can write
1tlogL(t)|t0≥{α0+γ1+ϑ222}(RSc−1)−qγ{α0+(1−p)γ1}(v+α0)⟨L(t)⟩ −qγ(α0+α1+γ2)(v+α0)⟨B(t)⟩+qγΥ(t)+ϑ2W2(t)t. |
Applying some algebraic manipulation with liminf yields
limt→∞inf⟨L(t)⟩≥(v+α0){α0+γ1+ϑ222}(RSc−1)qγ{α0+γ1(1−p)}=L2. | (4.17) |
Similarly, the direct use of Eq (4.13) gives that
1tlogB(t)|t0≤{α0+α1+γ2+ϑ232}(RSb−1) −(1−q)(γ+β+ξ)(α0+α1+γ2)(α0+v)⟨B(t)⟩ +(1−q)(γ+β+ξ)Υ(t)+ϑ3W3(t)t. |
Re-arrangement with some algebraic manipulation yields
⟨B(t)⟩≤(α0+v){α0+α1+γ2+ϑ232}(RSb−1)(1−q)(γ+β+ξ)(α0+α1+γ2) +(α0+v)(1−q)(γ+β+ξ)(α0+α1+γ2){(1−q)(γ+β+ξ)Υ(t)+ϑ3W3(t)t−1tlogB(t)|t0}. |
Through the implementation of lim with sup for an increasing t without bound, the last inequality may take the following form:
limt→∞sup⟨B(t)⟩≤(α0+v){α0+α1+γ2+ϑ232}(RSb−1)(1−q)(γ+β+ξ)(α0+α1+γ2)=B1. | (4.18) |
Likewise, from Eq (4.10), we obtain that
1tlogB(t)|t0≥{α0+α1+γ2+ϑ232}(RSd−1) −(1−q)β(α0+α1+γ2)v+α0⟨B(t)⟩+(1−q)βΥ(t)+ϑ3W3(t)t. |
As t approaches ∞, taking lim and inf, we can obtain the following assertion from the last inequality
limt→∞inf⟨B(t)⟩≥{α0+α1+γ2+ϑ232}(RSd−1)=B2. | (4.19) |
In addition, from the last equation of the integral system of the model (2.1), we can write
⟨C(t)⟩=η1α⟨L(t)⟩+η2α⟨B(t)⟩−1α{C(t)−C(0)t−ϑ5t∫t0C(y)dy}. |
By taking lim as t increases without bound and applying sup and inf respectively to the above equation, we obtain the following assertions
limt→∞sup⟨C(t)⟩≤η1αL1+η2αB1=C1, | (4.20) |
and
limt→∞inf⟨C(t)⟩≥η1αL2+η2αB2=C2. | (4.21) |
Thus, Eqs (4.16)–(4.21) imply that L2≤liminf⟨L(t)⟩≤limsup⟨L(t)⟩≤L1, B2≤liminf⟨B(t)⟩≤limsup⟨B(t)⟩≤B1 and C2≤liminf⟨C(t)⟩≤limsup⟨C(t)⟩≤C1, if t tends to ∞ and RSa>1, RSb>1, RSc>1 and RSd>1, which proves the conclusion.
To perform the numerical assessment of the stochastic epidemiological model (2.1), we chose to use the numerical Milstein higher-order method; for details, see [23]. Here, we will first demonstrate the model discretization that is used for the numerical investigation.
For, this, we assume that if K≥0 and the interval of time is [0,T] with the time step Δt=T/K, then li=iΔt is the point of the discretized interval. We also use (Si,Li,Bi,Ri,Ci) while writing (S(li),L(li),B(li),R(li),C(li)) and Wj(li)=Wji,j=1,2,3,4,5 for simplicity. Thus, the truncation of the Itˆo-Taylor expansion gives the procedure of Milstein's scheme for the reported model (2.1):
Si=Si−1+{Π−γSj−1Li−1−βSj−1Bi−1−ξSj−1Ci−1−(v+α0)}Δt+ϑ1Si−1(W1i−W1i−1)+12ϑ21Si−1{(W1i−W1i−1)2−Δt},Li=Li−1+{q(γSj−1Li−1+βSj−1Bi−1+ξSj−1Ci−1)−(α0+γ1)Li−1}Δt+ϑ2Li−1(W2i−W2i−1)+12ϑ22Li−1{(W2i−W2i−1)2−Δt},Bi=Bi−1+{(1−q)(γSj−1Li−1+βSj−1Bi−1+ξSj−1Ci−1)−(α0+α1+γ2)Bi−1}Δt+ϑ3Bj−1(W3i−W3i−1)+12ϑ23Bi−1{(W3i−W3i−1)2−Δt},Ri=Ri−1+{vSi−1−α0Ri−1+qγ1Li−1+γ2Bi−1}Δt+ϑ4Ri−1(W3i−W3i−1)+12ϑ24Ri−1{(W4i−W4i−1)2−Δt},Ci=Ci−1+{η1Li−1−αCi−1+η2Bi−1}Δt+ϑ5Cj−1(W5i−W5i−1)+12ϑ25Ci−1{(W5i−W5i−1)2−Δt}. |
To calculate the discretized Brownian paths that will be applied for Wj(li)−Wj(li−1) for the model (2.1) discretization, we assume the step size Δt to be an integral multiple of R≥1 of the increment δt. So the increments Wj(li)−Wj(li−1) becomes
Wj(lj)−Wi(lj−1)=iR∑k=iR−R+1dWk. |
Thus, the implementation of Milstein's higher order method for the considered model can be concluded with the following steps.
(1) Interval descretization of [0,T] into subintervals δt=TK>0: 0=l0<l1<⋯<lK=T in K with equal width with ln=nδt.
(2) Defining the proper initial data for the compartmental population groups of the model (S0,L0,B0,R0,C0).
(3) Defining the recursive formulas (Sj,Lj,Bj,Rj,Cj) for 1≤i≤K as described through the procedure of Milstein's scheme.
(4) Discretizing the Brownian paths to be used in the calculation of Wj(li)−Wj(li−1) while using the constant R≥1 and increment Δt.
(5) Finding of Wj(li)−Wj(li−1)=Wj(iRδt)−Wj((i−1)Rδt)=∑iRk=iR−R+1dWk.
We shall present the numerical experiment of the model with the aid of Milstein's Higher Order method as concluded in the above algorithm. According to Milstein's higher order method, we can assume two different sets of parameter values based on the sufficient analysis and calculations of the conditions that satisfy the derived analytical results in Theorems 4.1 and 4.2. We also assume some feasible biological initial sizes of the groups of the compartmental populations while implementing the algorithm. More precisely, the set of values of the epidemiological parameters for extinction analysis is assumed to be A1={Π,γ,β,ξ,α0,v,γ1,α1,γ2,p,q,η1,η2,α,ϑ1,ϑ2,ϑ3,ϑ4,ϑ5}, where the values are assigned as Π=0.9, γ=0.68, β=0.03, ξ=0.1, v=0.40, α=0.34, ϑ1=0.32, ϑ2=0.31, ϑ3=0.451, ϑ4=0.731, ϑ5=0.64, α=0.6, q=0.51, γ1=0.35, p=0.022, γ2=0.06, α1=0.65, η1=0.43 and η2=0.87, while the densities for the initial data were adjusted to be S(0)=0.50, L(0)=0.40, B(0)=0.30, R(0)=0.20, C(0)=0.10. Using these parametric values and initial data, we chose to implement Milstein's higher order scheme for the proposed epidemiological model (2.1) in Matlab and produce the output as depicted in Figures 2–6, which respectively demonstrate the temporal dynamics of susceptible, latent or asymptomatic and symptomatic, recovered individuals, as well as the amount of reservoir in the long term. This investigates the model extinction, which reveals that in the case of extinction analysis, the calculation of the stochastic reproductive parameters implies that RSa=0.68, RSb=0.96, RSc=0.57, and RSd=0.041, while the infected population, as well as the amount of reservoirs, goes to zero (vanishes) in the long term, as shown in Figures 3, 4 and 6. However, there will be always a non-infected population in the case of model extinction as shown in Figures 2 and 10.
On the other hand, the set A2 is assumed to be the set of parameters whose values for the model persistence have been adjusted to be γ=0.68, β=0.3, ξ=0.1, v=0.40, Π=0.9, α=0.14, ϑ1=0.31, ϑ2=0.41, ϑ3=0.251, ϑ4=0.731, ϑ5=0.64, α=0.6, q=0.51, γ1=0.35, p=0.022, γ2=0.06, α1=0.65, η1=0.43 and η2=0.87, while the initial data are same as that taken in the case of model extinction. We accordingly calculated the values of the reproductive parameters to be RSa=1.59, RSb=3.80, RSc=1.00 and RSd=1.05. With these parametric values and initial data, the implementation of Milstein's higher order algorithm generates the outputs shown in Figures 7–11, which depict the temporal dynamics of the model population whenever RSa>1, RSb>1, RSc>1 and RSd>1. We observed that in the case that the values of the reproductive numbers are greater than unity, there will always be infected individuals and an amount of reservoir in the community as shown in Figures 8, 9, and 11 respectively. Accordingly, the disease will persist and reach an endemic level. In addition, we also observed a strong impact of noise in coronavirus 2 transmission, where increasing noise intensity leads to the extinction; however, there is an inverse relation between the persistence and noise intensity.
In this work, we investigated the dynamics of SARS-CoV-2 virus transmission under the effects of asymptomatic and symptomatic infection phases, and various sources of randomness by assuming variation in all groups of the model compartmental populations. Keeping in view the asymptotic behavior of the disease, we categorized the infected group of individuals into two subcategories namely, asymptomatic and symptomatic populations, as both are responsible for the transmission of the disease. More precisely, we classified the total populations, into four epidemiological groups of susceptible, asymptomatically infected, symptomatically infected, and recovered individuals, and we assumed that both the asymptomatic and symptomatic individuals are causes of diseases transmission and the production of environmental reservoirs. We also used the probability-based disease transmission coefficient because susceptible individuals either become asymptomatic or symptomatic after getting infected. The model has been developed and studied in terms of its biological and mathematical feasibility to prove that the considered problem is well-possed. We then calculated the threshold quantities and derived the conditions for the model extinction and persistence. Finally, on the basis of theoretical results, we developed the algorithm with the aid of Milstein's higher order procedure to verify the model and show the validity of the results that have been obtained in the theoretical investigation of the model. We observed that the asymptotic behavior of the infected individuals and the white noise intensity are very influential in the analysis of the dynamics of SARS-CoV-2 virus transmission. The extinction of the disease is directly related to the intensity of noise while there is an inverse relationship between persistence and the intensity of white noise.
In the future, we will incorporate suitable time dependent control measures into the proposed model by using optimal control theory to design a mechanism that leads to elimination of the infection from the community. We will also study the fractional version of the model with the aid of fractional calculus. Further, we will use the Euler-Maruyama scheme and exponential Euler-Maruyama method for simulation purposes.
The authors declare that they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors would like to extend their appreciation to the Deanship of Scientific Research at King Khalid University for supporting this work through a research group program under grant no. R.G.P2/202/44. Further, the author Muhammad Bilal Riaz is thankful to the Ministry of Education, Youth and Sports of the Czech Republic for their support through the e-INFRA CZ (ID: 90254).
The authors declare no conflict of interest.
[1] |
Z. Wu, J. M. McGoogan, Characteristics of and important lessons from the coronavirus disease 2019 (Covid-19) outbreak in China: Summary of a report of 72 314 cases from the chinese center for disease control and prevention, JAMA, 323 (2020), 1239–1242. https://doi.org/10.1001/jama.2020.2648 doi: 10.1001/jama.2020.2648
![]() |
[2] |
L. Ferretti, C. Wymant, M. Kendall, L. Zhao, A. Nurtay, L. Abeler-Dörner, et al., Quantifying SARS-CoV-2 transmission suggests epidemic control with digital contact tracing, Science, 368 (2020), eabb6936. https://doi.org/10.1126/science.abb6936 doi: 10.1126/science.abb6936
![]() |
[3] |
A. V. Kamyad, R. Akbari, A. A. Heydari, A. Heydari, Mathematical modeling of transmission dynamics and optimal control of vaccination and treatment for hepatitis B virus, Comput. Math. Method. M., 2014 (2014), 475451. https://doi.org/10.1155/2014/475451 doi: 10.1155/2014/475451
![]() |
[4] |
S. Cai, Y. Cai, X. Mao, A stochastic differential equation sis epidemic model with two correlated brownian motions, Nonlinear Dynam., 97 (2019), 2175–2187. https://doi.org/10.1007/s11071-019-05114-2 doi: 10.1007/s11071-019-05114-2
![]() |
[5] |
A. Nwankwo, D. Okuonghae, A mathematical model for the population dynamics of malaria with a temperature dependent control, Differ. Equat. Dyn. Sys., 30 (2022), 719–748. https://doi.org/10.1007/s12591-019-00466-y doi: 10.1007/s12591-019-00466-y
![]() |
[6] |
A. Din, A. Khan, Y. Sabbar, Long-term bifurcation and stochastic optimal control of a triple-delayed Ebola virus model with vaccination and quarantine strategies, Fractal Fract., 6 (2022), 578. https://doi.org/10.3390/fractalfract6100578 doi: 10.3390/fractalfract6100578
![]() |
[7] |
P. R. S. Rao, M. N. Kumar, A dynamic model for infectious diseases: the role of vaccination and treatment, Chaos Soliton. Fract., 75 (2015), 34–49. https://doi.org/10.1016/j.chaos.2015.02.004 doi: 10.1016/j.chaos.2015.02.004
![]() |
[8] |
A. Akgül, S. H. Khoshnaw, A. S. Abdalrahman, Mathematical modeling for enzyme inhibitors with slow and fast subsystems, Arab J. Basic Appl. Sci., 27 (2020), 442–449. https://doi.org/10.1080/25765299.2020.1844369 doi: 10.1080/25765299.2020.1844369
![]() |
[9] |
A. Omame, D. Okuonghae, R. Umana, S. Inyama, Analysis of a co-infection model for HPV-TB, Appl. Math. Model., 77 (2020), 881–901. https://doi.org/10.1016/j.apm.2019.08.012 doi: 10.1016/j.apm.2019.08.012
![]() |
[10] |
D. Li, F. Wei, X. Mao, Stationary distribution and density function of a stochastic SVIR epidemic model, J. Franklin I., 359 (2022), 9422–9449. https://doi.org/10.1016/j.jfranklin.2022.09.026 doi: 10.1016/j.jfranklin.2022.09.026
![]() |
[11] |
S. Majee, S. Jana, S. Barman, T. K. Kar, Transmission dynamics of monkeypox virus with treatment and vaccination controls: A fractional order mathematical approach, Phys. Scripta, 98 (2023), 024002. https://doi.org/10.1088/1402-4896/acae64 doi: 10.1088/1402-4896/acae64
![]() |
[12] |
Y. Sabbar, A. Din, D. Kiouach, Influence of fractal-fractional differentiation and independent quadratic lévy jumps on the dynamics of a general epidemic model with vaccination strategy, Chaos Soliton. Fract., 171 (2023), 113434. https://doi.org/10.1016/j.chaos.2023.113434 doi: 10.1016/j.chaos.2023.113434
![]() |
[13] |
X. Zhai, W. Li, F. Wei, X. Mao, Dynamics of an HIV/AIDS transmission model with protection awareness and fluctuations, Chaos Soliton. Fract., 169 (2023), 113224. https://doi.org/10.1016/j.chaos.2023.113224 doi: 10.1016/j.chaos.2023.113224
![]() |
[14] |
A. J. Kucharski, T. W. Russell, C. Diamond, Y. Liu, J. Edmunds, S. Funk, et al., Early dynamics of transmission and control of COVID-19: A mathematical modelling study, Lancet Infect. Dis., 20 (2020), 553–558. https://doi.org/10.1016/S1473-3099(20)30144-4 doi: 10.1016/S1473-3099(20)30144-4
![]() |
[15] |
P. Samui, J. Mondal, S. Khajanchi, A mathematical model for COVID-19 transmission dynamics with a case study of india, Chaos Soliton. Fract., 140 (2020), 110173. https://doi.org/10.1016/j.chaos.2020.110173 doi: 10.1016/j.chaos.2020.110173
![]() |
[16] |
F. Bozkurt, A. Yousef, D. Baleanu, J. Alzabut, A mathematical model of the evolution and spread of pathogenic coronaviruses from natural host to human host, Chaos Soliton. Fract., 138 (2020), 109931. https://doi.org/10.1016/j.chaos.2020.109931 doi: 10.1016/j.chaos.2020.109931
![]() |
[17] |
A. G. M. Selvam, J. Alzabut, D. A. Vianny, M. Jacintha, F. B. Yousef, Modeling and stability analysis of the spread of novel coronavirus disease COVID-19, Int. J. Biomath., 14 (2021), 2150035. https://doi.org/10.1142/S1793524521500352 doi: 10.1142/S1793524521500352
![]() |
[18] |
T. Khan, G. Zaman, Y. El Khatib, Modeling the dynamics of novel coronavirus (COVID-19) via stochastic epidemic model, Results Phys., 24 (2021), 104004. https://doi.org/10.1016/j.rinp.2021.104004 doi: 10.1016/j.rinp.2021.104004
![]() |
[19] |
C. Bender, An Itô formula for generalized functionals of a fractional brownian motion with arbitrary hurst parameter, Stoch. Proc. Appl., 104 (2003), 81–106. https://doi.org/10.1016/S0304-4149(02)00212-0 doi: 10.1016/S0304-4149(02)00212-0
![]() |
[20] |
Q. Lei, Z. Yang, Dynamical behaviors of a stochastic SIRI epidemic model, Appl. Anal., 96 (2017), 2758–2770. https://doi.org/10.1080/00036811.2016.1240365 doi: 10.1080/00036811.2016.1240365
![]() |
[21] | Y. El-Khatib, A. Hatemi-J, Option valuation and hedging in markets with a crunch, J. Econ. Stud., 44 (2017). https://doi.org/10.1108/JES-04-2016-0083 |
[22] | K. Chung, The strong law of large numbers, The Regents of the University of California, 2008,145–156. |
[23] |
D. J. Higham, An algorithmic introduction to numerical simulation of stochastic differential equations, SIAM Rev., 43 (2001), 525–546. https://doi.org/10.1137/S0036144500378302 doi: 10.1137/S0036144500378302
![]() |
1. | Naveed Shahid, Ali Raza, Sana Iqbal, Nauman Ahmed, Emad Fadhal, Baboucarr Ceesay, Stochastic delayed analysis of coronavirus model through efficient computational method, 2024, 14, 2045-2322, 10.1038/s41598-024-70089-z |