images
Computer Modeling in Engineering & Sciences
images

DOI: 10.32604/cmes.2024.052383

ARTICLE

Numerical Analysis of Bacterial Meningitis Stochastic Delayed Epidemic Model through Computational Methods

Umar Shafique1,*, Mohamed Mahyoub Al-Shamiri2, Ali Raza3, Emad Fadhal4,*, Muhammad Rafiq5,6 and Nauman Ahmed5,7

1Department of Mathematics, National College of Business Administration and Economics, Lahore, 54660, Pakistan
2Department of Mathematics, Applied College, Mahayl Assir, King Khalid University, Abha, 62529, Saudi Arabia
3Department of Physical Sciences, The University of Chenab, Gujrat, 50700, Pakistan
4Department of Mathematics & Statistics, College of Science, King Faisal University, P.O. Box 400, Al-Ahsa, 31982, Saudi Arabia
5Department of Computer Science and Mathematics, Lebanese American University, Beirut, 1102-2801, Lebanon
6Department of Mathematics, Faculty of Science and Technology, University of Central Punjab, Lahore, 54000, Pakistan
7Department of Mathematics and Statistics, University of Lahore, Lahore, 54000, Pakistan
*Corresponding Authors: Umar Shafique. Email: umarshafique5919@gmail.com; Emad Fadhal. Email: efadhal@kfu.edu.sa
Received: 31 March 2024; Accepted: 28 June 2024; Published: 20 August 2024

Abstract: Based on the World Health Organization (WHO), Meningitis is a severe infection of the meninges, the membranes covering the brain and spinal cord. It is a devastating disease and remains a significant public health challenge. This study investigates a bacterial meningitis model through deterministic and stochastic versions. Four-compartment population dynamics explain the concept, particularly the susceptible population, carrier, infected, and recovered. The model predicts the nonnegative equilibrium points and reproduction number, i.e., the Meningitis-Free Equilibrium (MFE), and Meningitis-Existing Equilibrium (MEE). For the stochastic version of the existing deterministic model, the two methodologies studied are transition probabilities and non-parametric perturbations. Also, positivity, boundedness, extinction, and disease persistence are studied rigorously with the help of well-known theorems. Standard and nonstandard techniques such as Euler Maruyama, stochastic Euler, stochastic Runge Kutta, and stochastic nonstandard finite difference in the sense of delay have been presented for computational analysis of the stochastic model. Unfortunately, standard methods fail to restore the biological properties of the model, so the stochastic nonstandard finite difference approximation is offered as an efficient, low-cost, and independent of time step size. In addition, the convergence, local, and global stability around the equilibria of the nonstandard computational method is studied by assuming the perturbation effect is zero. The simulations and comparison of the methods are presented to support the theoretical results and for the best visualization of results.

Keywords: Bacterial Meningitis disease; stochastic delayed model; stability analysis; extinction and persistence; computational methods

1  Introduction

Meningitis caused by bacteria results in inflammation of the membranes surrounding the brain and spinal cord. The most common bacteria that cause meningitis are Streptococcus pneumoniae, Neisseria meningitidis, and Haemophilus influenzae type b. The signs and symptoms of bacterial meningitis include a sharply rising temperature, headache, stiff neck, and light sensitivity. Serious complications such as brain damage or death are caused by bacterial meningitis, so swift medical attention is crucial. Specific bacterial strains that cause meningitis can be controlled by effective vaccination [1]. Medical officers usually prescribe medicines such as ceftriaxone or cefotaxime to target invading bacteria and cure bacterial meningitis. Corticosteroids such as dexamethasone may be dispensed to decrease inflammation. Shin et al. [2] studied a mathematical model incorporating hospitalization and treatment factors. Madaki et al. [3] designed a mathematical model of meningococcal meningitis in Nigeria. In [4,5], the study was based on acute bacterial meningitis in adults. In [6,7], the study was designed on reviewed reports from the WHO and a model that describes a bacterial disease. Mohanty et al. [8] explained the different types of vaccines and preventive measures to control bacterial meningitis in Sweden. Van et al. [9] investigated in great detail all the stages of bacterial meningitis, from its initiation to its outbreak, and formulated various preventive measures to control it. Rauti et al. [10] developed a mathematical model for preventing bacterial meningitis in animals that can be employed to prevent infections of this deadly disease. Hou et al. [11] presented a dynamic model for preventing bacterial meningitis in young children using metagenomic sequencing, in which this approach can largely control the disease. Signing et al. [12] created a mathematical model that can visualize the changes caused by bacterial meningitis disease and their effect and can be treated accordingly to control the disease. Afridi et al. [13] identified different preventive measures and vaccines while reviewing bacterial meningitis cases brought to a civil hospital in Pakistan. Buonomo et al. [14] reported that the meningitis outbreak in Nigeria is analyzed using a behavior-focused vaccination model. This strategy considers behavioral elements while evaluating immunization strategies. It contributes to developing more effective preventative methods tailored to Nigeria’s socioeconomic environment and sheds light on the dynamics of meningitis transmission. Babatunde et al. [15] employed an online tool using logistic regression to predict cerebrospinal meningitis incidence is currently being developed. The goal of this initiative is to enhance illness management and early diagnosis. The initiative uses statistical modeling to evaluate relevant data, which can give useful information to policymakers and healthcare practitioners. Bradley et al. [16] reviewed the English medical literature and introduced a method that can be utilized to control bacterial meningitis. Choi et al. [17] created a mathematical model using AI to control bacterial meningitis in the future. Chen et al. [18] stated that it is essential to identify clinical and laboratory parameters to control bacterial meningitis disease. Olu et al. [19] created a mathematical model that can be very useful to control the growing epidemic of bacterial meningitis in the population. In [20], the authors studied a new analyzing technique for nonlinear time fractional Cauchy reaction-diffusion model equations. The stochastic epidemic models with delays are constructed given the need to incorporate the proper transmission process of the epidemic in any populace. These models make use of such lags and lags to mimic the period within which a person gets exposed and becomes infectious and lags in the uptake and efficacy of the implemented interventions. These models help to explain why so many epidemics exhibit cyclic behavior, this is in terms of the waves and seasonality that characterize many epidemics, and delay, noise, and heterogeneity are shown to play crucial roles. Also, the stochastic with delay effect results close to the actual statistics of the problem. So, in the present study, we extend the deterministic model in the sense of stochastic with delay. Then our new system was based on stochastic delay differential equations (SDDEs), due to complexity, we used standard and nonstandard methods to understand the dynamics of bacterial meningitis.

The structure of the paper is as follows: a brief overview and thorough analysis of bacterial meningitis-like disorders found in the literature are provided in Section 1. Establishing the delayed model and the ensuing mathematical analysis are the focus of Section 2. Section 3 presents the stochastic formulation of the model and the extinction and persistence of the model. The nonstandard computational method with delay and stochastic is introduced in Section 4, along with its global and local stability assessment. Numerical simulations and the presentation of the results are the explicit focus of Section 5. The final opinions provide a comprehensive conclusion of the work in Section 6.

2  Formulation of Model

The dynamics of population theory is a foundation for the model’s development. The classes susceptible S(t), carrier C(t), infected I(t), and recovered R(t) add up to the population N(t), and they use the law of mass action to represent the condition of bacterial meningitis (Fig. 1).

Figure 1: Flow of bacterial meningitis [3]

The parameter values of the constants are shown as follows in Table 1.

The law of mass action specifies the continuous model based on the presumptions. The transmission flow of bacterial meningitis-type disease is determined using the nonlinear delayed differential equations (DDEs) as follows:

S′(t)=Λ−βSIe−µτN−(µ+ρ)S(t)+ωR(t).(1)

C′(t)=βSIe−µτN−(γ+µ)C(t).(2)

I′(t)=γC(t)−(τ1+δ+μ)I(t).(3)

R′(t)=τ1I(t)+ρS(t)−(ω+µ)R(t).(4)

with the initial assumptions S(0)≥0,C(0)≥0,I(0)≥0,R(0)≥0, as well as the constraint that t≥0 and τ≤t.

2.1 Model Properties

For the critical analysis of the model to continue, all of the variables S(t),C(t),I(t), and R(t) must be non-negative. As a result, every time t that is greater than or equal to 0 and less than or equal to τ inside a feasible zone is covered by the conclusions drawn from the model’s analysis.

ℋ={(S,C,I,R)ϵR+4:N(t)≤Λµ,S≥0,C≥0,I≥0,R≥0}.

2.2 Analysis of Model

This section presents a concise analysis of the equilibria in the delayed model for bacterial meningitis, i.e., Meningitis-Free Equilibrium (MFE−D0), and Meningitis Existing Equilibrium (MEE−D∗). Hence,

D0=(S0,C0,I0,R0)=(Λ(ω+µ)[(ω+µ)(µ+ρ)−ωρ],0,0,Λρ[(ω+µ)(µ+ρ)−ωρ]), and D∗=(S∗,C∗,I∗,R∗).

S∗=N[(Λ+τ1I∗)βγe−µτ+ρN(µ+γ)(τ1+δ+μ)]βγe−μτ[βγe−μτ+N(µ+ρ)], C∗=(τ1+δ+μ)I∗γ,I∗=Λβγe−μτ−Nμ(µ+γ)(τ1+δ+μ)(µ+γ)(τ1+δ+μ)−τ1βγe−μτ,

R∗=τ1βγe−μτ+ρN(µ+γ)(τ1+δ+μ)βγe−μτ.

2.3 Reproduction Number

The systems (1)–(4) incorporate the advanced matrix approach to compute the threshold number R0. Neglected S′, the following can be obtained:

[C′I′R′]=[0βSe−μτN0000000][CIR]−[(µ+γ)00−γ(τ1+δ+μ)00−τ1(ω+µ)][CIR]

F=[0βSe−μτN0000000],G=[(µ+γ)00−γ(τ1+δ+μ)00−τ1(ω+µ)].

FG−1|D0=[βγS0e−μτN(µ+γ)(τ1+δ+μ)βS0e−μτN(τ1+δ+μ)0000000].

The spectral radius of FG−1|D0, which is also referred to as the threshold number, is defined as follows:

ℛ0=Λ(ω+µ)βγe−μτ(µ+γ)(τ1+δ+μ)[(ω+µ)(µ+ρ)−ωρ].(5)

3  Transition Probabilities of the Delayed Model

Let us examine the vector U(t)=[S(t),C(t),I(t),R(t)]T, along with the corresponding event probabilities as outlined in Table 2.

Here,

E∗[ΔU]=∑i=110PiMi=[Λ−βS(t)I(t)e−µτ−(µ+ρ)S(t)+ωR(t)βS(t)I(t)e−µτ−(γ+µ)C(t)γC(t)−(τ1+δ+μ)I(t)τ1I(t)+ρS(t)−(ω+µ)R(t)]Δt.

Var=E∗[ΔUΔUT]=∑i=110Pi[Mi][Mi]T.

=[P1+P2+P3+P4+P5−P20−P3−P5−P2P2+P6+P7−P600−P6P6+P8+P9−P8−P3−P50−P8P3+P5+P8+P10]Δt.

drift =H(U(t),t)=E∗[ΔU]Δt, diffusion =K(U(t),t)=E∗[ΔUΔUT]Δt, so

dU(t)=H(U(t),t)dt+K(U(t),t)dB(t).(6)

Eq. (6) is the stochastic differential equation, where B(t) is the Brownian motion.

3.1 Euler Maruyama Method

This section discusses a conventional numerical technique for approximating a stochastic delayed model’s result of (6). In this regard, Iq={0,1,2,3,…,q} is admitted for each q ∈ N. Let N∈N, and with the effect of time ∆t divides the partition into equal intervals [0, T] with constant delay. As follows:

0=to<t1<t2<…<tN=T,

for each n ∈ IN. Needless to mention tn=τn, for each n ∈ IN. Moreover, Un=U(tn) is accepted whenever n ∈ IN and U = S, E, I, R. Also, it can be set

ΔBn =B(tn+1) −B(tn),∀n∈ IN−1.

The mean of each ∆ Bn follows a normal distribution with a variance of one and an average of zero.

3.2 Stochastic Delayed Model

The system of stochastic delayed differential equations (SDDEs) is a mathematical model that describes the evolution of a set of variables over time, where the equations involve both deterministic time delays and stochastic (random) components. Where the stochastic term σi:(i=1,2,3,4)d(B(t)) introduces randomness into the system of differential equations as follows:

dS(t)=[Λ−βS(t)I(t)e−µτN−(µ+ρ)S(t)+ωR(t)]dt+σ1S(t)d(B(t)).(7)

dC(t)=[βS(t)I(t)e−µτN−(γ+µ)C(t)]dt+σ2C(t)d(B(t)).(8)

dI(t)=[γC(t)−(τ1+δ+μ)I(t)]dt+σ3I(t)d(B(t)).(9)

dR(t)=[τ1I(t)+ρS(t)−(ω+µ)R(t)]dt+σ4R(t)d(B(t)).(10)

Subject to non-negative initial conditions S(t)=S0≥0, C(t)=C0≥0,I(t)=I0≥0,R(t)=R0≥0.

3.2.1 Positivity and Boundedness

Let us examine a probability space represented as (Ω,F,P) and filtered by {Ft}tϵR. All P-null sets that fulfill the criteria of being both right-continuous and growing are encompassed inside Fo. Represent

U(t)=(S(t),C(t),I(t),R(t)).(11)

The norm |U(t)|=S2(t)+C2(t)+I2(t)+R2(t), and denote O3,1(R4×(0,∞);R+).

This study analyzes the collection of all non-negative functions V (U, t) defined on R4×(0,∞), which possess continuous second-order differentiability in U and first-order differentiability in t. Within this specific framework, the differential operator Q is precisely defined, which is associated with a three-dimensional stochastic differential equation.

dU(t)=H(U,t)dt+K(U,t)dW(t).(12)

As, Q=∂∂t+∑i=14Hi(U,t)∂∂ui+12∑i,j=14(KT(U,t)K(U,t))i,j×∂2∂Ui∂Uj

If Q acts on a function V∈O3,1(R4×(0,∞);R+), then

QV(U,t)=Vt(U,t)+VU(U,t)H(U,t)+12Trace(KT(U,t)VUU(U,t)K(U,t)),

where T means transportation.

Theorem 1. For a given model (7)–(10) and initial values (S(0),C(0),I(0),R(0))∈R+4, there exists a single solution (S(t),C(t),I(t),R(t)) on t≥0,τ≤t, and this solution will always remain in R+4 with a probability of one.

Proof: Based on Ito’s formula, the model (7)–(10) has a unique local positive solution in the interval [0,τe], and the time at which it explodes is represented by τe. Since all the coefficients in the model above satisfy the local Lipschitz criterion.

Next, this study demonstrates that the provided model (7)–(10) possesses a solution in the global sense, specifically when τe approaches infinity with practically certain probability.

Assume that mo=0 is a suitably high value such that S(0),C(0),I(0), and R(0) fall within the range [1mo,mo]. Let’s define a sequence of interrupting times for each integer m≥mo.

τm=inf{tϵ[0,τe]:S(t)∉(1m,m)orC(t)∉(1m,m)orI(t)∉(1m,m)orR(t)∉(1m,m)}.(13)

where this study sets infϕ=∞(ϕ represents the empty set). Since τm is non-decreasing as m→∞,

τ∞=limm→∞τm.(14)

Thereafter, the limit τ∞ is less than or equal to τe nearly surely. Now, this study must demonstrate that τ∞ approaches infinity with almost certain probability.

If this condition is not met, then there are positive values T and ε, where T is more than zero, and ε is between 0 and 1, such that

P{τ∞≤T}>ε.(15)

there is an integer m1>m0 such that

P{τm≤T}≥ε,∀m≥m1.(16)

Define a O4− function f:R+4→R+ by

f(S,C,I,R)=(S−1−ln⁡S)+(C−1−ln⁡C)+(I−1−ln⁡I)+(R−1−ln⁡R).(17)

The following can be calculated using Ito’s formula:

df(S,C,I,R)=(1−1S)dS+(1−1C)dC+(1−1I)dI+(1−1R)dR+σ12+σ22+σ32+σ422dt.

df(S,C,I,R)=(1−1S)((Λ−βS(t)I(t)e−µτN−(µ+ρ)S(t)+ωR(t))dt+σ1S(t)d(B(t)))+(1−1C)((βS(t)I(t)e−µτN−(γ+µ)C(t))dt+σ2C(t)d(B(t)))+(1−1I)((γC(t)−(τ1+δ+μ)I(t))dt+σ3I(t)d(B(t)))+(1−1R)(τ1I(t)+ρS(t)−(ω+µ)R(t))dt+σ4R(t)d(B(t))

df(S,C,I,R)≤[Λ+4μ+ρ+τ1+δ+γ+ω+σ12+σ22+σ32+σ422]dt+[σ1S(t)+σ2C(t)+σ3I(t)σ4R(t)]d(B(t)).(18)

For simplification, this study assumes

N1=Λ+4μ+ρ+τ1+δ+γ+ω+σ12+σ22+σ32+σ422, Then Eq. (18) could be written as

df(S,C,I,R)≤N1dt+[σ1S(t)+σ2C(t)+σ3I(t)σ4R(t)]d(B(t)).(19)

Assume N1 be a positive constant. Integrating from 0 to τm∧τ yields

∫0τm∧τdf(S(t),C(t),I(t),R(t))≤∫0τm∧τN1ds+∫0τm∧τ[σ1S(t)+σ2C(t)+σ3I(t)σ4R(t)]d(B(t)).(20)

where τm∧τ=min(τm,T), the taking the expectations leads to

Ef(S(τm∧τ),C(τm∧τ),I(τm∧τ),R(τm∧τ))≤f(S(0),C(0),I(0),R(0))+N1T.(21)

Set Ωm={τm≤T} for m>m1 and from (16), then P(Ωm≥ε). For every l∈Ωm there are some i such that ui(τm,l) equals either m or 1m for i=1,2,3.

Hence, f(S(τm,l),C(τm,l),I(τm,l),R(τm,l)) is less than min{m−1−ln⁡m,1m−1−ln⁡1m}.

Thereafter, the following can be obtained:

f(S(0),C(0),I(0),R(0))+N1T≥E(IΩm(l)V(S(τm),C(τm),I(τm),R(τm)))≥{min(m−1−ln⁡m,1m−1−ln⁡1m)}.(22)

The indicator function is represented by IΩm(l) of Ωm. Letting m→∞ leads to the contradiction

∞=f(S(0),C(0),I(0),R(0))+N1T<∞.

as desired.

3.2.2 Extinction and Persistence

Definition: Consider B(t) as a Brownian motion and I(t) as an Ito drift-diffusion process that obeys the stochastic differential equation:

dI(t)=[γC(t)−(τ1+δ+μ)I(t)]dt+σ3I(t)d(B(t)).

If f(I,t)∈C2(ℛ2,ℛ) then f(I(t),t) is also an Ito drift-diffusion process, which satisfies as follows:

d(f(I(t),t))=∂f∂t(I(t),t)dt+f′((I(t),t))dB(t)+12f′′((I(t),t))dB(t)2.

Let us introduce ℛ0S=ℛ0d−σ322(τ1+δ+μ).

Theorem 2. If ℛ0S < 1 and σ32<Λ(ω+µ)βγSe−μτ(µ+γ)(τ1+δ+μ)[(ω+µ)(µ+ρ)−ωρ], then the number of infected individuals in the system will exponentially approach zero (limt→∞I(t)=0).

Proof: Given the initial data (S(0),C(0),I(0),R(0))∈R+4, this study determines if the system (7)–(10) fulfills the stochastic differential equation and allows for a solution in the form of (S(t),C(t),I(t),R(t)).

dI(t)=[γC(t)−(τ1+δ+μ)I(t)]dt+σ3cI(t)d(B(t)).

Let, f(I)=ln⁡(I), then

dln(I)=f′(I)dI+12f′′(I)I2σ2dt.

dln⁡(I)=1IdI+12(−1I2)I2σ2dt.

dln⁡(I)=(γC(t)I(t)−(τ1+δ+μ)−12σ32)dt+σ3cd(B(t)).

ln⁡(I)=ln⁡I(0)+∫0t(γC(t)I(t)−(τ1+δ+μ)−12σ32)dt+∫0tσ3cd(B(t)),

Notice that, W(t)=∫0tσ3cd(B(t)) with W(0)=0.

Ifσ32>Λ(ω+µ)βγSe−μτ(µ+γ)(τ1+δ+μ)[(ω+µ)(µ+ρ)−ωρ],

ln⁡(I)>(Λ(ω+µ)βγSe−μτ(µ+γ)(τ1+δ+μ)[(ω+µ)(µ+ρ)−ωρ]−(τ1+δ+μ)−12Λ(ω+µ)βγSe−μτ(µ+γ)(τ1+δ+μ)[(ω+µ)(µ+ρ)−ωρ])t+W(t)+ln⁡I(0),

ln⁡(I)t>(Λ(ω+µ)βγSe−μτ2(µ+γ)(τ1+δ+μ)[(ω+µ)(µ+ρ)−ωρ]−(τ1+δ+μ))+W(t)t+ln⁡I(0)t,

limt→∞ln⁡(I)t>(Λ(ω+µ)βγSe−μτ2(µ+γ)(τ1+δ+μ)[(ω+µ)(µ+ρ)−ωρ]−(τ1+δ+μ))>0,withlimt→∞W(t)t=0,

Ifσ32<Λ(ω+µ)βγSe−μτ(µ+γ)(τ1+δ+μ)[(ω+µ)(µ+ρ)−ωρ], then

ln⁡(I(t))<(Λ(ω+µ)βγSe−μτ(µ+γ)(τ1+δ+μ)[(ω+µ)(µ+ρ)−ωρ]−(τ1+δ+μ)−12σ32)t+W(t)+ln⁡I(0),

ln⁡(I)t<(τ1+δ+μ)(Λ(ω+µ)βγSe−μτ(µ+γ)(τ1+δ+μ)2[(ω+µ)(µ+ρ)−ωρ]−1)+W(t)t+ln⁡I(0)t,

limt→∞supln⁡(I)t<(τ1+δ+μ)(ℛ0S−1), when ℛ0S<1, then limt→∞supln⁡(I)t≤0 can be obtained,

limt→∞I(t)=0, as desired.

4  Nonstandard Computational Method

The proposed non-parametric perturbation model’s first Eq. (10) can be stated using an unconventional computing technique; specifically, Eqs. (10)–(13) could be solved using a stochastic nonstandard finite difference technique.

dS(t)=[Λ−βS(t)I(t)e−µτN−(µ+ρ)S(t)+ωR(t)]dt+σ1S(t)d(B(t)).(23)

For the stochastic NSFD approach, the structure of the equation is as follows:

Sn+1−Snh=[Λ−βSnIne−µτN−(µ+ρ)Sn+ωRn+σ1SnΔBn](24)

The system (10)–(13) can be broken down using the stochastic NSFD process, and the complete system can then be written as follows:

Sn+1=Sn+h[Λ+ωRn+σ1SnΔBn]1+h(βIne−µτ+(µ+ρ))(25)

Cn+1=Cn+h[βSnIne−µτ+σ2CnΔBn]1+h(γ+μ)(26)

In+1=In+h[γCn+σ3InΔBn]1+h(τ1+δ+μ)(27)

Rn+1=Rn+h[τ1In+σ4RnΔBn]1+h(ω+μ)(28)

where n=0,1,2,…, and the discretization gap is denoted by “h.”

4.1 Convergence Analysis

Theorem 3. There is only one positive solution (S,C,I,R)∈R4,∀n>0 for any initial value (S(0),C(0),I(0),R(0))∈R4 for Eq. (25) through Eq.(28) with ΔBn=0.

Proof: The proof is easily verifiable due to the non-negative nature of the restriction of biological problems.

Theorem 4. For the region ℋ={(Sn,Cn,In,Rn)∈ℛ+4:Sn+Cn+In+Rn=N≤Λµ,Sn≥0,Cn≥0,In≥0,Rn≥0}. For every, n≥0 is an area of equations that is feasible and positive invariant (25) to (28) with ΔBn=0.

Proof: The system (25) to (28) can be deconstructed as follows:

Sn+1−Snh=[Λ−βSnIne−µτN−(µ+ρ)Sn+ωRn],

Cn+1−Cnh=[βSnIne−µτN−(γ+µ)Cn],

In+1−Inh=[γCn−(τ1+δ+μ)In],

Rn+1−Rnh=[τ1In+ρSn−(ω+µ)Rn],

After adding the above system of equations, the following can be obtained:

(Sn+1+Cn+1+In+1+Rn+1)−(Sn+Cn+In+Rn)h≤Λ−µ(Sn+Cn+In+Rn).

(Sn+1+Cn+1+In+1+Rn+1)−(Sn+Cn+In+Rn)≤hΛ−hµ(Sn+Cn+In+Rn).

(Sn+1+Cn+1+In+1+Rn+1)≤Λµ, as desired.

4.2 Local Stability Analysis of Nonstandard Computational Method

Consider the right-hand side of the Eqs. (25)–(28) functions F, G, H, and P and ΔBn=0.

F=S+h[Λ+ωR]1+h(βIe−µτ+(µ+ρ)),G=C+h[βSIe−µτ]1+h(γ+μ),H=I+h[γC]1+h(τ1+δ+μ),P=R+h[τ1I+ρS]1+h(ω+μ).

It is commonly understood that a system of the forms (25)–(28) converges to the model’s optimal state if and only if the Jacobian’s spectral radius, J,

JD=[∂F∂S∂F∂E∂F∂I∂F∂R∂G∂S∂G∂E∂G∂I∂G∂R∂H∂S∂P∂S∂H∂E∂P∂E∂H∂I∂P∂I∂H∂R∂P∂R](29)

∂F∂S=11+h(βIe−µτ+(µ+ρ)),∂F∂C=0,∂F∂I=−hβe−µτ(S+h[Λ+ωR])(1+h(βIe−µτ+(µ+ρ)))2,∂F∂R=hω1+h(βIe−µτ+(µ+ρ)),∂G∂S=0,∂G∂C=11+h(γ+μ),∂G∂I=hβSe−µτ1+h(γ+μ),∂G∂R=0,∂H∂S=0,∂H∂C=hγ1+h(τ1+δ+μ),∂H∂I=11+h(τ1+δ+μ),∂H∂R=0,∂P∂S=hρ1+h(ω+μ),∂P∂C=0,∂P∂I=hτ11+h(ω+μ),∂P∂R=11+h(ω+μ).

Theorem 5. The meningitis-free equilibrium (MFE−D0), D0=(S0,C0,I0,R0) of the nonstandard computational method is stable asymptotically in the sense of local if ℛ0<1 and ΔBn=0.

Proof: The Jacobian matrix (29) at D0 is

JD0=[11+h(µ+ρ)0−hβe−µτ((Λ(ω+µ)[(ω+µ)(µ+ρ)−ωρ])+h[Λ+ω(Λρ[(ω+µ)(µ+ρ)−ωρ])])(1+h(µ+ρ))2hω1+h(µ+ρ)011+h(γ+μ)hβ(Λ(ω+µ)[(ω+µ)(µ+ρ)−ωρ])e−µτ1+h(γ+μ)00hγ1+h(τ1+δ+μ)11+h(τ1+δ+μ)0hρ1+h(ω+μ)0hτ11+h(ω+μ)11+h(ω+μ)](30)

λ1=1−h2ρω1+h(µ+ρ)<1,λ2=1−h2βγ(Λ(ω+µ)[(ω+µ)(µ+ρ)−ωρ])e−µτ1+h(γ+μ)<1,λ3=11+h(ω+μ)<1,λ4=11+h(τ1+δ+μ)<1.

Using the definition of R0, it is possible to show that if R0<1, then λ2<1, and D0 is L.A.S. on the contrary, it is obviously to verify that λ2>1, if R0>1, which shows that D0 is unstable.

Theorem 6. The meningitis existing equilibrium (MEE−D∗), D∗=(S∗,C∗,I∗,R∗) of the nonstandard computational method is stable asymptotically in the sense of local if ℛ0>1 and ΔBn=0.

Proof: The Jacobian matrix (29) at D∗ is

JD∗=[11+h(βI∗e−µτ+(µ+ρ))0−hβe−µτ(S∗+h[Λ+ωR∗])(1+h(βI∗e−µτ+(µ+ρ)))2hω1+h(βI∗e−µτ+(µ+ρ))011+h(γ+μ)hβS∗e−µτ1+h(γ+μ)00hγ1+h(τ1+δ+μ)11+h(τ1+δ+μ)0hρ1+h(ω+μ)0hτ11+h(ω+μ)11+h(ω+μ)]

λ1=11+h(τ1+δ+μ)<1, λ2=hβS∗e−µτ1+h(γ+μ)<1.

|JD∗−λ|=|(1+h(βI∗e−µτ+(µ+ρ)))−h2βγI∗e−µτ(S∗+h[Λ+ωR∗])(1+h(βI∗e−µτ+(µ+ρ)))(1−h2βγS∗e−µτ)−λhω1+h(βI∗e−µτ+(µ+ρ))hρ+h3βγτ1I∗e−µτ(1+h(ω+μ))(1−h2βγS∗e−µτ)11+h(ω+μ)−λ|=0

A1=TrceofJD∗=(1+h(βI∗e−µτ+(µ+ρ)))−h2βγI∗e−µτ(S∗+h[Λ+ωR∗])(1+h(βI∗e−µτ+(µ+ρ)))(1−h2βγS∗e−µτ)+11+h(ω+μ).

A2=DeterminentofJD∗=(((1+h(βI∗e−µτ+(µ+ρ)))−h2βγI∗e−µτ(S∗+h[Λ+ωR∗])(1+h(βI∗e−µτ+(µ+ρ)))(1−h2βγS∗e−µτ))×(11+h(ω+μ)))+((hω1+h(βI∗e−µτ+(µ+ρ)))(hρ+h3βγτ1I∗e−µτ(1+h(ω+μ))(1−h2βγS∗e−µτ)))

Lemma 1. For the quadratic equation λ2−A1λ+A2=0, |λi|<1, i=1,2. If and only if the following conditions are satisfied:

(i) 1+A1+A2>0.

(ii) 1−A1+A2>0.

(iii) A2<1.

Proof: The proof is straightforward.

4.3 Global Stability Analysis of Nonstandard Computational Method

Theorem 7. The meningitis-free equilibrium (MFE−D0), D0=(S0,C0,I0,R0) of the nonstandard computational method is stable asymptotically in the sense of global if ℛ0<1 and ΔBn=0.

Proof: Let us consider the Lyapunov function as follows:

ϰm+1=1h[S0f(SmS0)+Cm+Im+Rm].

where f(x)=x−1−ln(x)≥f(1)=0. From the system, it can be obtained:

Δϰm=ϰm+1−ϰm=1h[S0f(Sm+1S0)−S0f(SmS0)+(Cm+1−Cm)+(Im+1−Im)+(Rm+1−Rm)].

ϰm+1=ϰm=1h[(Sm+1−Sm)−S0ln⁡(Sm+1Sm)+(Sm+1−Sm)].

Using Entasu et al. [21] criterion ln⁡(S2S1)≥S2−S1S2.

ϰm+1−ϰm≤1h[(Sm+1−Sm)−S0(Sm+1−SmSm+1)+(Cm+1−Cm)+(Im+1−Im)+(Rm+1−Rm)].

ϰm+1−ϰm≤−ΛSm+1S0(Sm+1−S0)2−(γ+µ)Cm+1−(τ1+δ+μ)Im+1−(ω+µ)Rm+1.

Definitely, ϰm+1−ϰm≤0 if R0<1 and ϰm+1−ϰm=0 iff Sm+1=S0,Cm+1=Im+1=Rm+1=0. Thus, the system is globally asymptotically stable.

Theorem 8. The meningitis existing equilibrium (MEE−D∗), D∗=(S∗,C∗,I∗,R∗) of the nonstandard computational method is stable asymptotically in the sense of global if ℛ0>1 and ΔBn=0.

Proof: In order to analyze the sufficient condition for global stability, the Lyapunov functions are defined as follows:

ℱm+1=[S∗f(SmS∗)+C∗f(CmC∗)+I∗f(ImI∗)+R∗f(RmR∗)].

where f(a)=a−1−ln(a)≥f(1)=0. From the system, the following can be obtained:

Δℱm=ℱm+1−ℱm=S∗f(Sm+1S∗)−S∗f(SmS∗)+C∗f(Cm+1C∗)−C∗f(CmC∗)+I∗f(Im+1I∗)−I∗f(ImI∗)+R∗f(Rm+1R∗)−R∗f(RmR∗).

ℱm+1−ℱm=(Sm+1−Sm)−S∗ln⁡(Sm+1S∗)+(Cm+1−Cm)−C∗ln⁡(Cm+1C∗)+(Im+1−Im)−I∗ln⁡(Im+1I∗)+(Rm+1−Rm)−R∗ln⁡(Rm+1R∗).

ℱm+1−ℱm≤hSm+1[(Sm+1−S∗)(Λ−βSm+1Ime−µτN−(µ+ρ)Sm+1+ωRm)]+hCm+1[(Cm+1−C∗)(βSmIme−µτN−(γ+µ)Cm+1)]+hIm+1[(Im+1−I∗)(γCm−(τ1+δ+μ)Im+1)]+hRm+1[(Rm+1−R∗)(τ1Im+ρSm−(ω+µ)Rm+1)].

ℱm+1−ℱm≤−hΛ(Sm+1−S∗)2Sm+1S∗−hω(Sm+1−S∗)2Sm+1S∗−βSmIme−µτN(Cm+1−C∗)2Cm+1C∗−hγCm(Cm+1−C∗)2Cm+1C∗−h(τ1Im+ρSm)(Rm+1−R∗)2Rm+1R∗.

Hence, Δℱm is a monotonic decreasing sequence. Thus, D∗ is globally asymptotically stable. Applying the same technique indicates that limn→∞(ℱm+1−ℱm)=0.

5  Comparison

This section compares known numerical approaches to a stochastic nonstandard computational method with existing methods such as Euler Maruyama, stochastic Euler, and stochastic Runge Kutta in the literature for the particular real-world problem.

5.1 Results

Fig. 2a,b compares the infected class of the nonstandard computational method in the sense of delay with stochastic and the Euler Maruyama method. At ℎ = 0.01, both methods demonstrate convergence in Fig. 2a. When the step size increases to ℎ = 0.5, the Euler Maruyama method exhibits divergence, whereas the nonstandard computational method in the sense of delay with stochastic maintains convergence, as illustrated in Fig. 2b. Similarly, Fig. 2c,d compares the infected class of the nonstandard computational method in the sense of delay with stochastic and stochastic Euler methods. At ℎ = 0.01, both methods converge in Fig. 2c; however, with an increment in step size to ℎ = 0.7, the stochastic Euler method diverges, while the nonstandard computational method in the sense of delay with stochastic method-maintained convergence, as depicted in Fig. 2d. Figs. 2e,f compares the infected class of the nonstandard computational method in the sense of delay with stochastic and stochastic Runge Kutta methods. At ℎ = 0.1, convergence is a numerical analysis of the stochastic delayed bacterial meningitis epidemic model observed in both methods in Fig. 2e; however, with an increase in step size to ℎ = 1.0, the stochastic Runge Kutta method diverges. In contrast, the nonstandard computational method in the sense of delay with the stochastic method continues to converge, as indicated in Fig. 2f. Fig. 3a depicts the impact of delay on the susceptible class of the model at various values 𝜏 = 0.1, 0.2, 0.3, 0.4, 0.5. Fig. 3b shows the effect of delay on the infected class of the model at different values 𝜏 = 0.1, 0.2, 0.3, 0.4, 0.5, illustrating a gradual decrease in the disease from the infected class over time. Fig. 4 shows a comparison of the double delay effect with the reproduction number. The increase in delay value may decrease the value of the reproduction number and vice versa. The decrease in reproduction number shows that the disease has been controlled in the population.

Figure 2: (a) The comparison behavior of Euler Maruyama and stochastic NSFD for a particular infected class at h=0.01 (convergent); (b) The comparison behavior of Euler Maruyama and stochastic NSFD for a particular infected class at h=0.5 (divergent); (c) The comparison behavior of stochastic Euler and stochastic NSFD for a particular infected class at h=0.01 (convergent); (d) The comparison behavior of stochastic Euler and stochastic NSFD for a particular infected class at h=0.7 (divergent); (e) The comparison behavior of stochastic Runge Kutta and stochastic NSFD for a particular infected class at h=0.1 (convergent); (f) The comparison behavior of stochastic Runge Kutta and stochastic NSFD for a particular infected class at h=1 (divergent)

Figure 3: Delay effects plots: (a) The effect of delay of different values on the susceptible class at h=100; (b) The effect of delay of different values on the infected class at h=100

Figure 4: Time plot of the effect of time delay (τ) with reproduction number (R0)

6  Conclusion

The study thoroughly evaluates the mathematical analysis of delayed epidemic models for meningitis using reliable delay approaches. Subpopulations are divided by the model into four categories: susceptible, carrier, infected, and recovered. The model’s equilibria, feasible region, and reproduction number are investigated rigorously. After that, two techniques, transition probabilities, and non-parametric perturbation, are implemented successfully on the existing deterministic model of meningitis bacterial disease in the sense of delay effect. The critical points of the newly stochastic version of the model, such as positivity, boundedness, extinction (means on which condition disease will become extinct from the population), and persistence (means on which condition disease will persist in the population), are studied in support of theorems. The Euler Maruyama, stochastic Euler, and stochastic Runge Kutta methods are successfully implemented on both stochastic versions of the model. Unfortunately, as observed in a comparison section, these methods have many discrepancies, such as negative, unbounded, and inconsistent results, and even work for a short period. The nonstandard computation method is implemented on the non-parametric version of the deterministic model to overcome these discrepancies. In this study, the proposed method ensures the biological properties of real-world problems. The most effective, efficient, and precise technique is the nonstandard computational method with stochastic in the sense of delay. In addition, the local and global stabilities of nonstandard computational methods are studied rigorously around the equilibria of the model using the Lyapunov theory after assuming the perturbation is zero. Accordingly, the nonstandard computation method with delay and stochastic is a good agreement for studying such complex models, and their results are close to nature by fulfilling the biological properties of the model.

Acknowledgement: All authors are grateful for the valuable comments of anonymous reviewers. Furthermore, all authors reviewed the results and approved the final version of the manuscript. The authors extend their appreciation to the Deanship of Research and Graduate Studies at King Khalid University for funding this work through Large Research Project under Grant Number RGP2 /302/45. Also, this work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia (Grant Number A426).

Funding Statement: The authors extend their appreciation to the Deanship of Research and Graduate Studies at King Khalid University for funding this work through large Research Project under Grant Number RGP2 /302/45. Also, this work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia (Grant Number A426).

Author Contributions: Ali Raza: Conceptualization, Investigation, Data curation, Writing—original draft. Mohamed Mahyoub Al-Shamiri: Project administration, Conceptualization, Investigation, Funding acquisition. Emad Fadhal: Data curation, Resources, Investigation, Funding acquisition. Muhammad Rafiq and Nauman Ahmed: Resources, Supervision, Project administration. Umar Shafique: Software, Writing—original draft & review & editing. All authors reviewed the results and approved the final version of the manuscript.

Availability of Data and Materials: All the data used and analyzed is available in the manuscript.

Ethics Approval: Not applicable.

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

References

1. Kim KS. Acute bacterial meningitis in infants and children. Lancet Infect Dis. 2010;10(1):32–42. doi:10.1016/S1473-3099(09)70306-8.

2. Shin SH, Kim KS. Treatment of bacterial meningitis: an update. Expert Opin Pharmacother. 2012;13(15):2189–206. doi:10.1517/14656566.2012.724399.

3. Madaki UY, Shuaibu A, Umar MI. A mathematical model for the dynamics of bacterial meningitis (Meningococcal meningitisa case study of Yobe State Specialist Hospital, Damaturu, Nigeria. Gadau J Pure Appl Sci. 2023;2(2):113–29. doi:10.54117/gjpas.v2i2.19.

4. Klein M, Abdel-Hadi C, Bühler R, Grabein B, Linn J, Nau R, et al. German guidelines on community-acquired acute bacterial meningitis in adults. Neurol Res Pract. 2023;5(1):01–11. doi:10.1186/s42466-023-00264-6.

5. Chazuka Z, Madubueze CE, Chukwu CW, Chikuni SM. Stability analysis of a bacterial meningitis model with saturated incidence and treatment default. Math Models Comput Simul. 2023;15(2):323–37. doi:10.1134/S2070048223020187.

6. Crankson MV, Olotu O, Afolabi AS, Abidemi A. Modeling the vaccination control of bacterial meningitis transmission dynamics: a case study. Math Model Control. 2023;3(4):416–34. doi:10.3934/mmc.2023033.

7. Turkun C, Golgeli M, Atay FM. A mathematical interpretation for outbreaks of bacterial meningitis under the effect of time-dependent transmission parameters. Nonlinear Dyn. 2023;12(2):01–18.

8. Mohanty S, Kostenniemi UJ, Silfverdal SA, Salomonsson S, Iovino F, Sarpong EM et al. Increased risk of long-term disabilities following childhood bacterial meningitis in Sweden. JAMA Netw Open. 2024;7(1):402–15. doi:10.1001/jamanetworkopen.2023.52402.

9. Van RN, Tubiana S, Broucker T, Cedric J, Roy C, Meyohas MC, et al. Persistent headaches one year after bacterial meningitis: prevalence, determinants, and impact on quality of life. Eur J Clin Microbiol Infect Dis. 2023;42(12):1459–67. doi:10.1007/s10096-023-04673-y.

10. Rauti R, Navok S, Biran D, Tadmor K, Leichtmann-Bardoogo Y, Ron EZ, et al. Insight on bacterial newborn meningitis using a neurovascular-unit-on-a-chip. Microbiol Spectr. 2023;11(3):e0123323. doi:10.1128/spectrum.01233-23.

11. Hou Y, Zhang M, Jiang Q, Yang Y, Liu J, Yuan K, et al. Microbial signatures of neonatal bacterial meningitis from multiple body sites. Front Cell Infect Microbiol. 2023;13:1169101. doi:10.3389/fcimb.2023.1169101.

12. Signing F, Tsanou B, Lubuma J, Bowong S. Effects of periodic aerosol emission on the transmission dynamics of neisseria meningitis A. Math Comput Simul. 2023;204(3):43–70. doi:10.1016/j.matcom.2022.06.005.

13. Afridi JK, Karim R, Idrees M, Yar S, Lala G, Gul H. Assessment of neonatal meningitis in a unique presentation with sepsis. J Postgrad Med Inst. 2023;37(2):130–34.

14. Buonomo B, Marca RD. A behavioural vaccination model with application to meningitis spread in Nigeria. Appl Math Model. 2024;125(2):334–50. doi:10.1016/j.apm.2023.09.031.

15. Babatunde RS, Adigun O, Folasade O. Towards prediction of cerebrospinal meningitis disease occurrence using logistics regression-a web-based application. FUDMA J Sci. 2023;7(4):337–43. doi:10.33003/fjs-2023-0704-1957.

16. Bradley JS, Bulitta JB, Cook R, Yu PA, Iwamoto C, Hesse EM, et al. Central nervous system antimicrobial exposure and proposed dosing for anthrax meningitis. Clin Infect Dis. 2024;78(6):01–18. doi:10.1093/cid/ciae093.

17. Choi BK, Choi YJ, Sung M, Ha W, Chu MK, Kim WJ, et al. Development and validation of an artificial intelligence model for the early classification of the aetiology of meningitis and encephalitis: a retrospective observational study. Clin Med. 2023;61:34–45. doi:10.1016/j.eclinm.2023.102051.

18. Chen Y, Luo C, Zhou G, Wang H, Dai K, Wu W, et al. The discrimination between autoimmune glial fibrillary acidic protein astrocytopathy and tuberculous meningitis. Mult Scler Relat Disord. 2024;6(1):105527–32. doi:10.1016/j.msard.2024.105527.

19. Olu OT, Michael OO, Nathaniel OS, Abraham OA, Emmanuel FS, James AK, et al. Propagation dynamics of meningitis disease based on complex network modeling. Univ J Public Health. 2023;11(3):324–31. doi:10.13189/ujph.2023.110306.

20. Ahmad H, Khan TA, Ahmad I, Stanimirovic PS, Chu YM. A new analyzing technique for nonlinear time fractional Cauchy reaction-diffusion model equations. Results Phys. 2020;19:103462. doi:10.1016/j.rinp.2020.103462.

21. Enatsu Y, Nakata Y, Muroya Y, Izzo G, Vecchio A. Global dynamics of difference equations for SIR epidemic models with a class of nonlinear incidence rates. J Differ Equ Appl. 2012;18(7):1163–81. doi:10.1080/10236198.2011.555405.

images
Copyright © 2024 The Authors. Published by Tech Science Press.
This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.