[BACK]
Intelligent Automation & Soft Computing
DOI:10.32604/iasc.2022.024993
images
Article

Computational Approximations for Real-World Application of Epidemic Model

Shami A. M. Alsallami1, Ali Raza2,*, Mona Elmahi3, Muhammad Rafiq4, Shamas Bilal5, Nauman Ahmed6 and Emad E. Mahmoud7

1Department of Mathematical Sciences, College of Applied Science, Umm Al-Qura University, Makkah, 21955, Saudi Arabia
2Department of Mathematics, Govt. Maulana Zafar Ali Khan Graduate College Wazirabad, Punjab Higher Education Department (PHED), Lahore, 54000, Pakistan
3Mathematics Department, College of Science & Arts in Riyadh Alkkbra, Qassim University, Saudi Arabia
4Department of Mathematics, Faculty of Sciences, University of Central Punjab, Lahore, 54600, Pakistan
5Department of Mathematics, University of Sialkot, 51310, Sialkot, Pakistan
6Department of Mathematics and Statistics, The University of Lahore, 54590, Lahore, Pakistan
7Department of Mathematics, College of Science, Taif University, P. O. Box 11099, Taif, 21944, Saudi Arabia
*Corresponding Authors: Ali Raza. Email: Alimustasamcheema@gmail.com
Received: 07 November 2021; Accepted: 24 December 2021

Abstract: The real-world applications and analysis have a significant role in the scientific literature. For instance, mathematical modeling, computer graphics, camera, operating system, Java, disk encryption, web, streaming, and many more are the applications of real-world problems. In this case, we consider disease modeling and its computational treatment. Computational approximations have a significant role in different sciences such as behavioral, social, physical, and biological sciences. But the well-known techniques that are widely used in the literature have many problems. These methods are not consistent with the physical nature and even violate the actual behavior of the continuous model. The structural properties needed for such disciplines, as dynamical consistency, positivity, and boundedness, are the critical requirements of the models in these fields. We studied the transmission of Lassa fever dynamically and numerically. The model’s positivity, boundedness, reproduction number, equilibria, and local stability are investigated in dynamical analysis. In numerical analysis, we developed some explicit and implicit methods. Unfortunately, explicit methods like Euler and Runge Kutta are time-dependent and violate the physical properties of the disease. Then, the proposed implicit method for the said model, the non-standard finite difference, is independent of the time step, dynamically consistent, positive, and bounded. In the end, a comparison of methods is presented.

Keywords: Lassa fever disease; epidemic model; computational approximations; convergence analysis

1  Literature Survey

Lassa fever is an intense hemorrhagic illness caused by the Lassa virus (member of the arenavirus). Lassa virus carries a rat, which is very common in West Africa. It is also known as a zoonosis, which means that disease spreads from animal to human. People usually become infected with Lassa fever because of food and household items’ exposure to urine of infected Mastomys rats. Its symptoms are varied and include cardiac, neurological, and pulmonary problems. In some West Africa, this disease is endemic to the rodent population. It is most common in Liberia, Sierra Leone, Guinea, and Nigeria. Lassa fever is transmitted by playing, touching, and cutting up a rat’s dead body. In 2018, Usman et al. analyzed the Lassa fever virus infection for the transmission dynamics of qualitative and analytic activities of a mathematical model [1]. In 2020, Peter et al. studied Lassa fever’s dynamics, and the solution model stayed and verified boundedness and positivity of basic properties [2]. In 2017, Olayiwola et al. analyzed in the lowliest endemic countries people with an ultimate risk of infection and need for constant investigation develops commanding in the endemic region like Africa with Nigeria at the essential attention will help in no small processes in scheming the scourge [3]. Woyessa et al. investigated the private and public health facilities, control interventions, and prevent infection when feverish patients avoid nosocomial infections [4]. In 2013, Ajayi et al. have studied that Epidemic contained rejoinder schemes and testing to control exertions comprised fright between health staffs, insufficient/poor quality of defensive things, insufficient extra preparation, and poor local laboratory capability [5]. In 2019, Ilori et al. analyzed the activity supporting improving planned for patient care and LLC, emerging infectious diseases, and Medscape through applied epidemiological characteristics and clarifying factors associated with mortality [6]. Adewuyi et al. studied the disease have the probable of actuality organize an infectious menace that essential be controlled by way of biological weapon and currently, vaccine of Lassa fever no available and some natural problems occurs for development of vaccines so prevention the way by control the rodent [7]. Iroezindu et al. analyzed Lassa fever spread; challenges, letdown to use proper defensive tackle, stigmatization of associates, and absence of a purpose-built isolation facility [8]. In 2019, Amodu et al. have studied Lassa fever as an acute disease of scarceness, high endemicity by way of cooperated environmentally-friendly sanitation, and relics susceptible populations to community health problems in Nigeria. Public meeting defense for attractive prevention strategy remains particular sanitation [9]. In 2019, Kangbai et al. highlighted that seasonal epidemics use an effective treatment to make a stratagem, which can control procedures of Lassa fever prevention and control the connection between humans to rodents [10]. Makinde et al. investigated Lassa fever to identify when a nonconformity toward the prestige quo has happened [11]. In 2020, Zhao et al. analyzed quantify of this impact in Nigeria’s presence of Lassa fever significance measure the connotation among local precipitation and infection reproduction number which facts has probable elect applied as a bad sign for Lassa fever epidemics [12]. In 2017, Obabiyi et al. developed a mathematical model for transmission of Lassa fever dynamics with the behavior of susceptible humans, recovered humans divided the population into two parts such as human populations, and rodent population by using the positivity, bounded theorem, and suggested the stability hygiene of environment [13]. In 2018, Akpede et al. studied the necessities for achievement and enduring capacity of the control exertions in Nigeria and the sub-region. In wholly these, the Nigerian administration with NCDC necessity carries a huge responsibility for the organization, supply deployment, and support. If necessary, even persuade sub-regional administrations addicted to action. In addition, there should be expected through determined action [14]. In 2019, Mazzola et al. explored the Diagnostics necessary for acknowledging and controlling epidemics of LASV, unique prevailing and genetically various mediators of VHF, that use scenarios with different performance requirements for text complexity, sensitivity, specificity, and development time [15]. In 2019, Nwafor et al. examined the Lassa fever outbreak in Nigeria; the Health maintenance workers necessity, take a high index of doubt of the infection and follow IPC measures even though provided that maintenance for all patients. Explaining health maintenance workers by the new strategies mentioned above is also significant to reduce the menace of nosocomial transmission of Lassa fever [16]. In 2018, Shehu et al. studied that the Occurrence of rural to urban change of clinical and epidemiological reduced the Lassa fever cases during 2016 for morbidity and mortality [17]. In 2007, Ogbu et al. discussed the situation of Lassa fever in the sub-region of West Africa and suggested strategies for socioeconomic behavior that control the shortage of health care system [18]. In 2020, Tewogbola et al. analyzed the overview and discussed the main reasons it damaged the human population and recommended the control measure of Lassa fever [19]. In 2014, Ajayi et al. reported a case of 59 years that recovers without taking a vaccine such as ribavirin. The symptoms of this disease increase day by day because few people do not use the main precautions to control Lassa fever [20]. Some well-known numerical models related to diseases are studied [2126].

2  Formulation of Lassa Fever Model

The variables and parameters are described of the lassa fever model as follows: SH(t) : denoted as the susceptible class at any time t, IH(t) : characterized as the infectious class at any time t, RH(t) : characterized as the recovered class at any time t, SR(t) : characterized as the susceptible rodent vectors at any time t, IR(t) : characterized as the infectious rodent vectors at any time t, NH(t) : characterized as whole humans’ population at any time t, m = NRNH : characterized as the number of infectious rodent vectors by the human host, α1 : described as the rate at which contagious rodent vectors and a susceptible class of humans interact with each other, α2 : defined as the force of infection, α3 : defined as the rate at which sensitive rodent vectors and an infectious class of humans interact with each other, τc : denoted the speed at which infectious human hosts comply with the drug, τnc : indicated the rate at which infectious human hosts do not comply with the drug, rc : denoted the rate at which infectious human hosts are educated to adhere to the medication, δ: indicated the rate of mortality of an infectious class, γ : indicated the rate at which humans may lose their immunity. The leading equations of the model are as follows:

dSH(t)dt=ΛHα1α2SH(t)IR(t)NH+γRH(t)+τncIH(t)μHSH(t),t0. (1)

dIH(t)dt=α1α2SH(t)IR(t)NHτcIH(t)rcIH(t)τncIH(t)δIH(t)μHIH(t),t0. (2)

dRH(t)dt=τcIH(t)+rcIH(t)γRH(t)μHRH(t),t0. (3)

dSR(t)dt=ΛRα1α3SR(t)IH(t)NHμRSR(t),t0. (4)

dIR(t)dt=α1α3SR(t)IH(t)NHμRIR(t),t0. (5)

where SH(0)0,IH(0)0,RH(0)0,SR(0)0,IR(0)0 .

2.1 Fundamental Properties of Model

We studied the feasible region, positivity, and boundedness of the model at any time, t ≥ 0, as follows:

Π={(SH,IH,RH,SR,IR)εR+5:SH+IH+RHΛHμH,SR+IRΛRμR,SH0,IH0,RH0,SR0,IR0} .

Lemma 1: The solutions (SH,IH,RH,SR,IR)εR+5 of Eqs. (1)(5) is positive at any time t ≥ 0, with given non-negative initial conditions.

Proof: It is clear from Eqs. (1)(5),

dSHdt|SH=0=ΛH+γRH+τncIH0 , dIHdt|IH=0=α1α2SHIRNH0 , dRHdt|RH=0=τcIH+rcIH0 , dSRdt|SR=0=ΛR0 , dIRdt|IR=0=α1α3SRIHNH0, as desired.

Lemma 2: Forgiven any non-negative initial conditions for the solution of the system (1)(5) is bounded if limtSupNH(t)ΛHμH,limtSupNR(t)ΛRμR.

Proof: let us consider the population function as follows:

NH(t)=SH+IH+RH,dNHdt=dSHdt+dIHdt+dRHdt,dNHdt=ΛHμHNH,

NH(t)=A+ΛHμH

By the Gronwall’s inequality, we get

NH(t)NH(0)+ΛHμH,t0,limtSupNR(t)ΛHμH

NR(t)=SR+IR,dNRdt=dSRdt+dIRdt,dNRdt=ΛRμRNR

NR(t)=B+ΛRμR

By the Gronwall’s inequality, we get

NR(t)NR(0)+ΛRμR,t0

limtSupNR(t)ΛRμR , as desired.

2.2 Steady States of Lassa fever Model

There are two steady states of Eqs. (1) to (5), as follows: disease-free equilibrium (DFE)=(SH,IH,RH,SR.IR)=(ΛHμH,0,0,ΛRμR,0) and endemic equilibrium (EE)=(SH,IH,RH,SR,IR) ,

RH=(τc+rc)IHγ+μH=A1IH,A1=(τc+rc)γ+μH,SH=ΛH+γA1IHA2IHμH,A2=τc+rc+δ+μH,IR=α1α3SRIHμR,

SR=ΛRα1α3IH+μR,IH=ΛHA4μRA4α1α3γA1+A2,A3=τc+rc+τnc+δ+μH,A4=A3μHμRα12α2α3ΛR.

3  Reproduction Number of Lassa Fever Model

In this section, we shall find the two types of matrices like transmission and transition by assuming the disease-free equilibria in the system (1)(5) by using the next-generation matrix method, furthermore, considering the infected classes as follows:

F=[α1α2ΛHμHNH0000000α1α2ΛRμRNH], V1=[μR(γ+μH)00μR(τc+rc)μR(τc+rc+τnc+δ+μH)000(τc+rc+τnc+δ+μH)(γ+μH)](τc+rc+τnc+δ+μH)(γ+μH)μR .

FV1=1(τc+rc+τnc+δ+μH)(γ+μH)μR[α1α2ΛHμHNH(μR)(γ+μH)0000000α1α2ΛRμRNH(τc+rc+τnc+δ+μH)(γ+μH)].

Thus, the dominant eigenvalue of the matrix is called reproduction number and denoted as follows:

R0=α1α2ΛHμHNH(τc+rc+τnc+δ+μH) .

4  Local Stability

In this section, we present two well-known theorems for stability in the sense of local. Again, consider the system (1)(5) as function of A,B,C,DandE as follows:

A=ΛHα1α2SHIRNH+γRH+τncIHμHSH. (6)

B=α1α2SHIRNHτcIHrcIHτncIHδIHμHIH. (7)

C=τcIH+rcIHγRHμHRH. (8)

D=ΛRα1α3SRIHNHμRSR. (9)

E=α1α3SRIHNHμRIR. (10)

Theorem 1:The disease-free equilibrium is locally asymptotically stable for the system (6)(10). If R0<1 and otherwise unstable if R0>1 .

Proof:First, we take the partial derivates of the system (6)(10) concerning state variables as follows:

ASH=α1α2IRNHμH,AIH=τnc,ARH=γ,ASR=0,AIR=α1α2SHNH .

BSH=α1α2IRNH,BIH=τcrcτncδμH,BRH=0,BSR=0.

BIR=α1α2SHNH,CSH=0CIH=τc+rc,CRH=γμH,CSR=0,

CIR=0,DSH=0,DIH=α1α3SRNH,DRH=0,DSR=α1α3IHNHμR,DIR=0,FSH=0,

FIH=α1α3SRNH,FRH=0,FSR=α1α3IHNH,FIR=μR.

Here, the Jacobian matrix as follows:

J=[α1α2IRNHμHτncγ0α1α2SHNHα1α2IRNHτcrcτncδμH00α1α2SHNH0τc+rcγμH000α1α3SRNH0α1α3IHNHμR00α1α3SRNH0α1α3IHNHμR]

The Jacobian matrix at the disease-free equilibria of the system (6)(10) as follows:

J(E0)=J(ΛHμH,0,0,ΛRμR,0)=[μHτncγ0α1α2ΛHμHNH0τcrcτncδμH00α1α2ΛHμHNH0τc+rcγμH000α1α3ΛRμRNH00μR00α1α3ΛRμRNH00μR] .

|JλI|=|μHλτncγ0α1α2ΛHμHNH0(τcrcτncδμH)λ00α1α2ΛHμHNH0τc+rcγμHλ000α1α3ΛRμRNH00μRλ00α1α3ΛRμRNH00μRλ|=0.

λ1=μH<0,λ2=μR<0,λ3=(γ+μH)<0.

|JλI|=|(τcrcτncδμH)λα1α2ΛHμHNHα1α3ΛRμRNHμRλ|=0 .

λ2+λ(a1+μR)+a1μRa2=0 .

where, a1=τc+rc+τnc+δ+μH,a2=(α1α2ΛHμHNH)(α1α3ΛRμRNH)

Since all the coefficients of the polynomial are positive, therefore, by using Routh Hurwitz Criteria for 2nd order, the disease-free equilibria are locally asymptotically stable.

Theorem 2:The endemic equilibrium is locally asymptotically stable for the system (6)(10) if R0>1 .

Proof: The Jacobian matrix at the endemic equilibria of the system (6)(10) is as follows:

J(E) = J(SH,IH,RH,SR,IR)=[α1α2IRNHμHτncγ0α1α2SHNHα1α2IRNHτcrcτncδμH00α1α2SHNH0τc+rcγμH000α1α3SRNH0α1α3IHNHμR00α1α3SRNH0α1α3IHNHμR] .

|JλI|=|B1μHλτncγ0B2B1B3λ00B20B4B5λ000B60B7μRλ00B60B7μRλ|=0 .

where, B1=α1α2IRNH,B2=α1α2SHNH,B3=τc+rc+τnc+δ+μH,B4=τc+rc,B5=γ+μH,B6=α1α3SRNH,B7=α1α3IHNH.

|JλI|=(B1μHλ)|B3λ00B2B4B5λ00B60B7μRλ0B60B7μRλ|B1|τncγ0B2B4B5λ00B60B7μRλ0B60B7μRλ|=0 .

|JλI|=(B1μHλ)(B5+λ)[(B3λ)|B7μRλ0B7μRλ|+B2|B6B7μRλB6B7|]+(B1)(γ)[(B4)|B6B7μRλB6B7|]+(B1)(B5+λ)[τnc|B7μRλ0B7μRλ|B2|B6B7μRλB6B7|]=0

λ5+λ4[2μR+B3+μH+B7+B5+B1]+λ3[B7μR+μR2+2B3μR+2B5μR+2B1μR+B3B7+B5B7+B1B7+B3B5+B1B3B2B6+B5μH+B3μH+B7μH+2μRμH+B1B5B1τnc]+λ2[B3B7μR+B5B7μR+B1B7μR+B3μR2+B5μR2+B1μR2+2B3B5μR+2B1B3μR+2B1B5μR+B3B5B7+B1B3B7+B1B5B7+B1B3B5+B1B4γ2B1μRτncB1B7τncB1B5τncB1B2B6+B3B7μH+B7B5μH2B2B6μH+2B3μRμH+2B5μRμH+B7μRμHB2B5B6]+λ[B3B5B7μR+B1B3B7μR+B1B5B7μR+B3B5μR2+B5μR2+B1B5μR2+2B1B3B5μRB2B6B5μR2B1B2B5B6+2B1B4γμRB1B2B6γ2B1B5τncμRB1B5B7τnc+B7B3B5μH+2B3B5μRμH2B1B2B7B6B1B2B6γ2B1B5τncμRB1B5B7τnc+B7B3B5μH+2B3B5μRμH2B1B2B7B62B1B2B6μR+B5B3B7μH+B7B5μRμH+B3B7μRμHB2B6μRμHB2B6B5μHB1B2B6]+B1B3B5B7μR+B1B3B5μR2+B1B4B7γμR+B1B4γμR2B1B5B7τncμR+B1B2B6γμRB1B5τncμR2+B1B3B7B52B1B2B5B6μR2B1B2B5B7B6B2B5B6μRμH+B5B3B7μRμH+B3B5μR2μH=0.

λ5a1+λ4a2+λ3a3+λ2a4+λa5+a6=0 .

where, a1=[1],a2=[2μR+B3+μH+B7+B5+B1].

a3=[B7μR+μR2+2B3μR+2B5μR+2B1μR+B3B7+B5B7+B1B7+B3B5+B1B3B2B6+B5μH+B3μH+B7μH+2μRμH+B1B5B1τnc]

a4=[B3B7μR+B5B7μR+B1B7μR+B3μR2+B5μR2+B1μR2+2B3B5μR+2B1B3μR+2B1B5μR+B3B5B7+B1B3B7+B1B5B7+B1B3B5+B1B4γ2B1μRτncB1B7τncB1B5τncB1B2B6+B3B7μH+B7B5μH2B2B6μH+2B3μRμH+2B5μRμH+B7μRμHB2B5B6].

a5=[B3B5B7μR+B1B3B7μR+B1B5B7μR+B3B5μR2+B5μR2+B1B5μR2+2B1B3B5μRB2B6B5μR+2B1B2B5B6μR+2B1B4γμRB1B2B6γ2B1B5τncμRB1B5B7τnc+B7B3B5μH+2B3B5μRμH2B1B2B7B62B1B2B6μR+B5B3B7μH+B7B5μRμH+B3B7μRμHB2B6μRμHB2B6B5μHB1B2B6]

a6=[B1B3B5B7μR+B1B3B5μR2+B1B4B7γμR+B1B4γμR2B1B5B7τncμR+B1B2B6γμRB1B5τncμR2+B1B3B7B52B1B2B5B6μR2B1B2B5B7B6B2B5B6μRμh+B5B3B7μRμH+B3B5μR2μH]

The Routh Hurwitz criteria of the 5th order are satisfied. Hence, the endemic equilibria are locally asymptomatically stable.

5  Computational Approximations

In this section, we present the well-known approximations like Euler, Runge Kutta, and non-standard finite difference for the system (1)(5) as follows:

5.1 Euler Approximation

The discretization of the system (1)(5) under the rules of the Euler approximation is as follows:

SHn+1=SHn+h[ΛHα1α2sHnIRnNH+γRHn+τncIHnμHSHn]. (11)

IHn+1=IHn+h[α1α2SHnIRnNHτcIHnrcIHnτncIHnδIHnμHIHn]. (12)

RHn+1=RHn+h[τcIHn+rcIHnγRHnμHRHn]. (13)

SRn+1=SRn+h[ΛRα1α3SRnIHnNHμRSRn]. (14)

IRn+1=IRn+h[α1α3SRnIHnNHμRIRn]. (15)

where h is any discretization parameter and n0.

5.2 Runge-Kutta Approximation

The discretization of the system (1)(5) under the rules of the Runge Kutta approximation is as follows:

Stage#1

K1=h[ΛHα1α2sHnIRnNH+γRHn+τncIHnμHSHn] .

L1=h[α1α2SHnIRnNHτcIHnrcIHnτncIHnδIHnμHIHn] .

M1=h[τcIHn+rcIHnγRHnμHRHn] .

O1=h[ΛRα1α3SRnIHnNHμRSRn] .

P1=h[α1α3SRnIHnNHμRIRn] .

Stage#2

K2=h[ΛHα1α2(SHn+K12)(IRn+P12)NH+γ(RHn+M12)+τnc(IHn+L12)μH(SHn+K12)] .

L2=h[α1α2(SHn+K12)(IRn+P12)NHτc(IHn+L12)rc(IHn+L12)τnc(IHn+L12)δ(IHn+L12)μH(IHn+L12)] .

M2=h[τc(IHn+L12)+rc(IHn+L12)γ(RHn+M12)μH(RHn+M12)] .

O2=h[ΛRα1α3(SRn+O12)(IHn+L12)NHμR(SRn+O12)] .

P2=h[α1α3(SRn+O12)(IHn+L12)NHμR(IRn+P12)] .

Stage#3

K3=h[ΛHα1α2(SHn+K22)(IRn+P22)NH+γ(RHn+M22)+τnc(IHn+L22)μH(SHn+K22)] .

L3=h[α1α2(SHn+K22)(IRn+P22)NHτc(IHn+L22)rc(IHn+L22)τnc(IHn+L22)δ(IHn+L22)μH(IHn+L22)] .

M3=h[τc(IHn+L22)+rc(IHn+L22)γ(RHn+M22)μH(RHn+M22)] .

O3=h[ΛRα1α3(SRn+O22)(IHn+L22)NHμR(SRn+O22)] .

P3=h[α1α3(SRn+O22)(IHn+L22)NHμR(IRn+P22)] .

Stage#4

K4=h[ΛHα1α2(SHn+K3)(IRn+P3)NH+γ(RHn+M3)+τnc(IHn+L3)μH(SHn+K3)] .

L4=h[α1α2(SHn+K3)(IRn+P3)NHτc(IHn+L3)rc(IHn+L3)τnc(IHn+L3)δ(IHn+L3)μH(IHn+L3)] .

M4=h[τc(IHn+L3)+rc(IHn+L3)γ(RHn+M3)μH(RHn+M3)] .

O4=h[ΛRα1α3(SRn+O3)(IHn+L3)NHμR(SRn+O3)] .

P4=h[α1α3(SRn+O3)(IHn+L3)NHμR(IRn+P3)] .

Final stage

SHn+1=SHn+16[K1+2K2+2K3+K4]. (16)

IHn+1=IHn+16[L1+2L2+2L3+L4]. (17)

RHn+1=RHn+16[M1+2M2+2M3+M4]. (18)

SRn+1=SRn+16[O1+2O2+2O3+O4]. (19)

IRn+1=IRn+16[P1+2P2+2P3+P4]. (20)

where h is any discretization parameter and n0.

5.3 Non-standard Finite Difference Approximation

The discretization of the system (1)(5) under the rules of the non-standard finite difference scheme is as follows [27]:

SHn+1=SHn+hΛH+γhRHn+hτncIHn1+hα1α2IRnNH+μHh. (21)

IHn+1=IHn+hα1α2SHnIRnNH1+hτc+rch+hτnc+δh+μHh. (22)

RHn+1=RHn+hτcIHn+rcIHn1+γh+μHh. (23)

SRn+1=SRn+hΛR1+hα1α3IHnNH+μRh. (24)

IRn+1=IRn+hα1α3SRnIHnNH1+μRh. (25)

where h is any discretization parameter and n0.

5.4 Convergence Analysis

Considering the functions, A, B, C, D, and E at the system (21)(25) as follows:

A=SH+hΛH+γhRH+hτncIH1+hα1α2IRNH+μHh,

 B=IH+hα1α2SHIRNH1+hτc+rch+hτnc+δh+μHh,

C=RH+hτcIH+rcIH1+γh+μHh, D=SR+hΛR1+hα1α3IHNH+μRh, E=IR+hα1α3SRIHNH1+μRh.

The elements of Jacobian matrix as follows:

ASH=11+μHh,AIH=hτnc1+μHh,ARH=γh1+μHh,ASR=0,AIR=(SH+hΛH+γhRH+hτncIH)(hα1α2NH)(1+hα1α2IRNH+μHh)2.

BSH=hα1α2IRNH1+hτc+rch+hτnc+δh+μHh,BIH=11+hτc+rch+hτnc+δh+μHh,BRH=0,BSR=0,

BIR=hα1α2SHNH1+hτc+rch+hτnc+δh+μHh,CSH=0,CIH=hτc+rc1+γh+μHh,CRH=11+γh+μHh,CSR=0,

CIR=0,DSH=0,DIH=(SR+hΛR)(hα1α3NH)(1+hα1α3IHNH+μRh)2,DRH=0,DSR=11+μRh,DIR=0,ESH=0,

EIH=hα1α3SRNH1+μRh,ERH=0,ESR=hα1α3IHNH1+μRh,EIR=11+μRh .

Theorem 3: For n0, the eigenvalues of the Jacobian matrix at disease-free equilibrium for the system (21)(25) lie in the unit circle if R0<1 .

Proof: The Jacobian matrix at disease-free equilibrium (DFE-E0 ) = (ΛHμH,0,0,ΛRμR,0) is as follows:

|J(E0)λ|=|11+μHhλhτnc1+μHhγh1+μHh0(ΛHμH+hΛH)(hα1α2NH)(1+μHh)2011+hτc+rch+hτnc+δh+μHhλ00hα1α2ΛHμHNH1+hτc+rch+hτnc+δh+μHh0hτc+rc1+γh+μHh11+γh+μHhλ000(ΛRμR+hΛR)(hα1α3NH)(1+μRh)2011+μRhλ00hα1α3ΛRμRNH1+μRh0011+μRhλ|=0 .

λ1=|11+μHh|<1,λ2=|11+γh+μHh|<1,λ3=|11+hμR|<1.

|J(E0)λ|=|(11+hτc+rch+hτnc+δh+μHh)λhα1α2ΛHμHNH1+hτc+rch+hτnc+δh+μHhhα1α3ΛRμRNH1+μRh11+μRhλ|=0 .

P1=TraceofJ=(11+hτc+rch+hτnc+δh+μHh)+11+μRh

P2=DeterminantofJ=((11+hτc+rch+hτnc+δh+μHh)(11+μRh))(hα1α3ΛRμRNH1+μRh)(hα1α2ΛHμHNH1+hτc+rch+hτnc+δh+μHh) .

Lemma 1:For the quadratic equation λ2P1λ+P2==0 , |λi|<1,i=1,2,3, if and only if the following conditions are satisfied:

(i). 1+P1+P2>0.

(ii). 1P1+P2>0.

(iii). P2<1.

Proof: The proof is straight forward.

Theorem 6: For n0, the eigenvalues of the Jacobian matrix at endemic equilibrium for the system (21)(25) lie in the unit circle if R0>1 .

Proof:The Jacobian matrix at the endemic equilibria E1=(SH,IH,RH,SR,IR) as follows:

JJ(E*)=J(SH*,IH*,RH*,SR*,IR*)=[11+μHhhτnc1+μHhγh1+μHh0(SH*+hΛH+γhRH*+hτncIH*)(hα1α2NH)(1+hα1α2IR*NHμHh)2hα1α2IR*NH1+hτc+rch+δh+μHh11+hτc+rch+δh+μHh00hα1α2SH*NH1+hτc+rch+δh+μHh0hτc+rc1+γh+μHh11+γh+μHh000(SR*+hΛR)(hα1α3NH)(1+hα1α3IH*NH+μRH)2011+μRh00hα1α3SR*NH1+μRh0hα1α3IH*NH1+μRh11+μRh].

Hence, the largest eigenvalue of the Jacobian is less than one, ultimately remaining will also lie in the unit circle when R0>1 . Thus, endemic equilibrium is stable.

5.5 Computational Approximations

By using the values of the parameters as presented in Tab. 1. The diagrams for the system (1)(5) for disease-free equilibrium (DFE) and endemic equilibrium (EE) plotted with MATLAB software as follows:

images

5.6 Comparison Section

imagesimages

Figure 5: Combine graphical behaviors of NSFD with Euler and Runge Kutta methods at different time-step sizes (a) Comparison of Euler and NSFD at h = 0.1 (b) Comparison of Euler and NSFD at h = 1 (c) Comparison of Runge Kutta and NSFD at h = 0.1 (d) Comparison of Runge Kutta and NSFD at h = 1

6  Results and Discussion

We investigated the transmission dynamics of Lassa fever disease in humans and rats through the study. The critical point is modeling, terminology related to epidemiology, and Lassa fever disease. Dynamical analysis of the model is investigated. Computational analysis, including well-known methods, is presented. Mostly, methods are valid for only tiny time step sizes. But inappropriately flop for huge time step sizes like Euler and Runge Kutta. Our proposed scheme (NSFD) remains convergent for step sizes like h = 100. Furthermore, Tab. 2 shows the efficiency of the numerical methods.

images

7  Conclusion

The non-standard finite difference scheme was designed for the communication dynamics of Lassa fever disease. Unfortunately, the earlier methods, like Euler and Runge Kutta of order 4th, are unsuitable because they depend on time step size. So, Euler and Runge Kutta are tentatively convergent. When we increase the time step size, the graph of Euler and Runge Kutta gives variation in result from time to time they display divergent. The new well-known numerical scheme, like the non-standard finite difference scheme independent of time step size. The NSFD scheme is a comfortable tool on behalf of dynamical properties like stability, positivity, boundedness and shows the exact behavior of the continuous model. The graphical behavior of ODE-45, Euler, Runge Kutta, NSFD schemes and comparison of schemes are given in Figs. 1a, 1b, Figs. 2a, 2b, Figs. 3a, 3b, Figs. 4a, 4b and Figs. 5a5d respectively. In the end, we could extend this idea to all types of nonlinear and complex models. In the future, we could develop our analysis into fuzzy epidemic models and many other types of modeling as given in [2731].

images

Figure 1: Combine graphical behaviors of the Lassa fever disease (a) Sub-populations at disease-free equilibrium (DFE) (b) Subpoulations at endemic equilibrium (EE)

images

Figure 2: Euler method for the behavior of infected rats at different time-step sizes (a) Infected rats at h = 0.01 (b) Infected rats at h = 1

images

Figure 3: Runge Kutta method for the behavior of infected rats at different time-step sizes (a) The behavior of infected rats at time step size h = 0.1 (b) The behavior of infected rats at time step size h = 1

images

Figure 4: NSFD method for the behavior of infected rats at different time-step sizes (a) The behavior of Infected rats for EE at h = 0.1 (b) The behavior of infected rats for EE at h = 1000

Acknowledgement: Thanks, our families and colleagues who supported us morally.

Funding Statement: The authors received no specific funding for this study.

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

References

 1.  S. Usman and I. I. Adamu, “Modelling the transmission dynamics of the Lassa fever infection,” Mathematical Theory and Modeling, vol. 8, no. 5, pp. 43–68, 2018. [Google Scholar]

 2.  O. J. Peter, A. I. Abioye, F. A. Oguntolu, T. A. Owolabi, M. O. Ajisope et al., “Modelling and optimal control analysis of Lassa fever disease,” Informatics in Medicine Unlocked, vol. 20, no. 1, pp. 01–11, 2020. [Google Scholar]

 3.  J. O. Olayiwola and A. S. Bakarey, “Epidemiological trends of Lassa fever outbreaks and insights for future control in Nigeria,” International Journal of Tropical Disease and Health, vol. 24, no. 4, pp. 01–14, 2017. [Google Scholar]

 4.  A. B. Woyessa, L. Maximore, D. Keller, J. Dogba, M. Pajibo et al., “Lesson learned from the investigation and response of Lassa fever outbreak,” BMC Infectious Diseases, vol. 19, no. 610, pp. 01–08, 2019. [Google Scholar]

 5.  N. A. Ajayi, C. G. Nwigwe, B. N. Azuogu, B. N. Onyire, E. U. Nwonwu et al., “Containing a Lassa fever epidemic in a resource-limited setting: Outbreak description and lessons learned from Abakaliki,” International Journal of Infectious Diseases, vol. 17, no. 1, pp. 1011–1017, 2013. [Google Scholar]

 6.  E. A. Ilori, Y. Furuse, O. B. Ipadeola, C. C. Dan-Nwafor, A. Abubakaret et al., “Epidemiologic and clinical features of Lassa fever outbreak in Nigeria,” Emerging Infectious Diseases, vol. 25, no. 6, pp. 1066–1076, 2019. [Google Scholar]

 7.  G. M. Adewuyi, A. Fowotade and B. T. Adewuyi, “Lassa fever: Another infectious menace,” Journal of Clinical Microbiology, vol. 10, no. 3, pp. 144–155, 2009. [Google Scholar]

 8.  M. O. Iroezindu, U. S. Unigwe, C. C. Okwara, G. A. Ozoh, A. C. Ndu et al., “Lessons learnt from the management of a case of Lassa fever and follow-up of nosocomial primary contacts in Nigeria during ebola virus disease outbreak in West Africa,” Tropical Medicine and International Health, vol. 20, no. 10, pp. 1424–1430, 2015. [Google Scholar]

 9.  S. E. Amodu and S. O. Fapohunda, “Lassa fever and the Nigerian experience: A review,” European Journal of Biological Research, vol. 9, no. 3, pp. 2449–8955, 2019. [Google Scholar]

10. J. B. Kangbai, F. K. Kamara, R. C. Lahai and F. Gebeh, “Lassa fever in post-ebola sierra leone socio-demographics and case fatality rates of in-hospital patients admitted at the Kenema Government hospital Lassa fever ward between 2016-2018,” Research Square, vol. 1, no. 1, pp. 01–17, 2019. [Google Scholar]

11. O. A. Makinde, “As ebola winds down, Lassa fever reemerges yet again in West Africa,” Journal of Infection in Developing Countries, vol. 10, no. 2, pp. 199–200, 2016. [Google Scholar]

12. S. Zhao, S. S. Musa, H. Fu, D. He and J. Qin, “Large-scale Lassa fever outbreaks in Nigeria: Quantifying the association between disease reproduction number and local rainfall,” Epidemiology and Infection, vol. 148, no. 1, pp. 01–14, 2019. [Google Scholar]

13. O. S. Obabiyi and A. A. Onifade, “Mathematical model for Lassa fever transmission dynamics with variable human and reservoir population,” International Journal of Differential Equations and Applications, vol. 16, no. 1, pp. 01–26, 2017. [Google Scholar]

14. G. O. Akpede, D. A. Asogun, S. A. Okogbenin and P. O. Okokhere, “Lassa fever outbreaks in Nigeria,” Expert Review of Anti-Infective Therapy, vol. 16, no. 9, pp. 663–667, 2018. [Google Scholar]

15. L. T. Mazzola and C. Kelly-Cirino, “Diagnostics for Lassa fever virus: A genetically diverse pathogen found in low-resource settings,” BMJ Global Health, vol. 4, no. 1, pp. 01–12, 2019. [Google Scholar]

16. C. C. D. Nwafor, Y. Furuse, E. A. Ilori, O. Ipadeola, K. O. Akabike et al., “Measures to control protracted large Lassa fever outbreak in Nigeria,” Euro Surveill, vol. 24, no. 1, pp. 01–05, 2019. [Google Scholar]

17. N. Y. Shehu, S. S. Gomerep, S. E. Isa, K. O. Iraoyah, J. Mafuka et al., “Nigeria-the changing epidemiology and clinical presentation,” Frontier Public Health, vol. 6, no. 1, pp. 01–06, 2018. [Google Scholar]

18. O. Ogbu, E. Ajuluchukwu and C. J. Uneke, “Lassa fever in West African sub-region: An overview,” Journal of Vector-Borne Diseases, vol. 44, no. 1, pp. 01–11, 2007. [Google Scholar]

19. P. Tewogbola and N. Aung, “Lassa fever: History, causes, effects, and reduction strategies,” International Journal of One Health, vol. 6, no. 2, pp. 95–99, 2020. [Google Scholar]

20. N. A. Ajayi, K. N. Ukwaja, N. A. Ifebunandu, R. Nnabu, F. I. Onwe et al., “Lassa fever-full recovery without ribavirin treatment: A case report,” African Health Sciences, vol. 14, no. 1, pp. 1072–1077, 2014. [Google Scholar]

21. W. Shatanawi, A. Raza, M. S. Arif, M. Rafiq, M. Bibi et al., “Essential features preserving dynamics of stochastic dengue model,” Computer Modeling in Engineering and Sciences, vol. 126, no. 1, pp. 201–215, 2021. [Google Scholar]

22. A. Raza, M. S. Arif, M. Rafiq, M. Bibi, M. Naveed et al., “Numerical treatment for stochastic computer virus model,” Computer Modeling in Engineering and Sciences, vol. 120, no. 2, pp. 445–465, 2019. [Google Scholar]

23. M. S. Arif, A. Raza, K. Abodayeh, M. Rafiq and A. Nazeer, “A numerical efficient technique for the solution of susceptible infected recovered epidemic model,” Computer Modeling in Engineering and Sciences, vol. 124, no. 2, pp. 477–491, 2020. [Google Scholar]

24. W. Shatanawi, A. Raza, M. S. Arif, M. Rafiq, M. Bibi et al., “Essential features preserving dynamics of stochastic dengue model,” Computer Modeling in Engineering and Sciences, vol. 126, no. 1, pp. 201–215, 2021. [Google Scholar]

25. M. A. Noor, A. Raza, M. S. Arif, M. Rafiq, K. S. Nisar et al., “Non-standard computational analysis of the stochastic COVID-19 pandemic model: An application of computational biology,” Alexandria Engineering Journal, vol. 61, no. 1, pp. 619–630, 2021. [Google Scholar]

26. K. Abodayeh, A. Raza, M. S. Arif, M. Rafiq, M. Bibi et al., “Numerical analysis of stochastic vector-borne plant disease model,” Computers, Materials and Continua, vol. 63, no. 1, pp. 65–83, 2020. [Google Scholar]

27. R. E. Mickens, “Dynamic consistency: A fundamental principle for constructing non-standard finite difference schemes for differential equations,” Journal of Difference Equations and Applications, vol. 11, no. 7, pp. 645–653, 2005. [Google Scholar]

28. A. Raza, A. Ahmadian, M. Rafiq, S. Salahshour and M. Ferrara, “An analysis of a nonlinear susceptible-exposed-infected-quarantine-recovered pandemic model of a novel coronavirus with delay effect,” Results in Physics, vol. 21, no. 1, pp. 01–07, 2021. [Google Scholar]

29. G. Neha, A. Ghosh, S. P. Mondal, M. Y. Bajuri, A. Ahmadian et al., “Identification of dominant risk factor involved in spread of COVID-19 using hesitant fuzzy MCDM methodology,” Results in Physics, vol. 21, no. 1, pp. 01–18, 2021. [Google Scholar]

30. Z. Muhammad, K. Shah, F. Nadeem, M. Y. Bajuri, Ali Ahmadian et al., “Threshold conditions for global stability of disease-free state of COVID-19,” Results in Physics, vol. 21, no. 1, pp. 01–08, 2021. [Google Scholar]

31. A. Ahmadian, N. Senu, F. Larki, S. Salahshour, M. Suleiman et al., “A Legendre approximation for solving a fuzzy fractional drug transduction model into the bloodstream,” Recent Advances on Soft Computing and Data Mining, vol. 1, no. 1, pp. 25–34, 2014. [Google Scholar]

images 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.