Mathematical Modelling of Rotavirus Disease Through Efficient Methods

Ali Raza

Department of Mathematics,Govt.Maulana Zafar Ali Khan Graduate College Wazirabad,Punjab Higher Education Department(PHED),Lahore,54000,Pakistan

Abstract: The design of evolutionary approaches has a vital role in the recent development of scientific literature.To tackle highly nonlinear complex problems,nonlinear ordinary differential equations,partial differential equations,stochastic differential equations, and many more may called computational algorithms.The rotavirus causes may include severe diarrhea,vomiting,and fever leading to rapid dehydration.By the report of the World Health Organization (WHO), approximately 600,000 children die worldwide each year,80 percent of whom live in developing countries.Two million children are hospitalized each year.In Asia,up to 45 percent of the children hospitalized for diarrhea may be infected with rotavirus.The rotavirus model is categorized into five-compartment like susceptible(S),breastfeeding(M),vaccinated(V),infected (I), and recovered children (R).Positivity, boundedness, equilibria,reproduction number,and stability results are part of the qualitative analysis of the model.After that, the design of the evolutionary approaches on the model predicts the efficiency,visualization,long-term behavior of the disease,and best results of the rotavirus disease.In the end,evolutionary computations are an appropriate tool for double-checking the qualitative examination of the model.

Keywords:Rotavirus disease;epidemic model;methods;comparison analysis

1 Literature Survey

Rotavirus is a double-stranded RNA virus.Rotavirus is the most common cause of diarrhea.Its spread is contagious and affects every child before the age of five.Rotavirus is round in shape when observed under a microscope.Rotavirus produces inflammation in the small bowel,especially in the stomach and intestine.Diarrhea is most common in South Asia and Africa.In 2004,2.5 billion cases of diarrhea were reported, resulting in 1.5 million deaths.More than half of the patients were from South Asia and Africa.Diarrhea is still a problem in the developing world.In America,the death rate from diarrhea is 10% merely while in South Asia and Africa, the death rate is 31.3%.In Pakistan,many children suffering from diarrhea are 27.7%in 2018.All the children were under five years of age.Burnett et al.presented a cost-effective vaccine to prevent the rotavirus disease[1].In 2020,Ahmed et al.developed a mathematical analysis of epidemic modeling of diarrhea like rotavirus[2].Lin et al.investigated the causes of diarrheal infectious disease rotavirus(Rotarix and Rotateq)in 2014[3].In 2019,Payne et al.presented a rotavirus vaccine to understand their effectiveness and recommend US infant immunization [4].Omondi et al.presented a mathematical model to explore the co-infection of malaria and rotavirus[5].In 2009,Effeltrerre et al.studied the dynamics of rotavirus to check the indirect effect of vaccination [6].Chan et al.presented epidemiology in 1998 of rotavirus infection to calculate and estimate the burden of disease under five years of age in Hong Kong [7].Linhares et al.investigated a rotavirus infection that affects children before five in Brazil by a longitudinal study in 2009[8].In 2010,Atchison et al.developed a determined age-structure model to investigate the transmission and effect of rotavirus vaccination in England and Wales[9].In 2013,Arnold et al.presented a theory of rotavirus-like infections [10].In 2020, Shuaib et al.developed a mathematical analysis to study the effect of rotavirus on childhood mortality and diarrheal disease[11].Guzel et al.presented a met-analysis in 2020 regarding awareness programs and preventive policies of rotavirus infection in younger children in Turkey [12].In 2020, Folorunso et al.investigated oral rotavirus vaccine in developed countries and its strategic preventive[13].Bibera et al.studied the evolution of the rotavirus vaccine[14].In 2020,Ilmi et al.investigated the dynamics of the rotavirus epidemic model with the effect of crowding infective individuals [15].In 2018, Shumetie et al.suggested improving child care to protect their children of a mother from diarrheal morbidity and rotavirus vaccination[16].Bennett studied the rotavirus epidemic model with the well-known assumption of mathematics[17].In 2016,Namawejje presented a mathematical model to check the affected regions of rotavirus and illustrate the effects on children and treatment[18].Omondi et al.studied a mathematical analysis and simulation of rotavirus to check the quality of vaccination in 2015, whether it’s working or not on the infection [19].The paper’s strategy is as follows: Section 2 describes the model’s design and qualitative analysis.Section 3 and 4 presented the model’s reproduction number and stability results,respectively.Section 5 investigated the structure of evolutionary approaches and their efficiency with comparative analysis.In the end,remarks and the conclusion of the paper are presented.

2 Model Formulation

For any time,the whole populationE(t):is divided into the five compartments named:susceptible(S), breastfeeding (M), vaccinated (V), infected (I), and recovered (R).The flow map of rotavirus disease is presented in Fig.1.

Figure 1:Flow transmission of rotavirus disease[11]

The fixed ratios of the model are presented as follows:Λ(1-σ-η):represent the inclusion rate into the susceptible compartment,σΛ: represent the inclusion rate into breastfeeding compartment,ηΛ: represent the inclusion rate into the vaccinated box, Ψ: represent the breastfeeding rate of the vaccinated box, γ: represent the vaccination rate of the susceptible compartment, Φ: represent the vaccination rate of breastfeeding compartment, β: represent the nominal contact rate, ?: illustrate the waning rate of maternal antibodies from breast milk, ω: represent the waning rate of vaccine, ?:represent the reduction in the risk of infection due to maternal antibodies,ξ:represent the reduction in the risk of infection due to vaccination,τ:represent the disease mortality rate,μ:represent the natural death rate and k:represent the recovery rate.The governing equations of the model are as follows:

Here,we can observe from the system Eqs.(1)-(5),there is no termR(t),in first four equations.So,by using the assumption of epidemiological modeling,we can consider the analysis of the model without recovered compartment.Therefore,

with nonnegative conditionsS(0)≥0,M(0)≥0,V(0)≥0,I(0)≥0.

2.1 Properties

Theorem 1:The system results(6)-(9)with given initial conditions are positive for allt≥0.

Proof:By letting the Eq.(6),

Similarly,for Eqs.(7)-(9),we have

Theorem 2:The solutions(S,M,V,of the system (6)-(9)are bounded at any time and limit→∞SupE(t)≤.

Proof:By considering the population function is as follows:

2.2 Equilibria

The system (6)-(9)admits two types of equilibria as follows: rotavirus free equilibrium (RFER1)=(S1,M1,V1,I1)=(S1,M1,V1,0),K1=μ[?(ω+γ+μ)+(φ+μ)[ω+Ψ+μ+γ)+Ψω].

rotavirus endemic equilibrium(REE-R2)=(S*,M*,V*,I*)

F1I*3+F2I*2+F3I*+F4=0.

F1=εξ(δ+κ+μ),F2=(δ+κ+μ)(μεξ+εξψ+εξγ+με+μξ+φξ+εω+ξ?)-Λβεξ,F3=(δ+κ+μ)(φω + ?ω + μφ + ξγ? + μξγ + φξψ + μξψ + εψω + μ2ξ + μεγ +μεω + μξ? + μ? + μω + μφξ + μ2+ φξγ + μ2ε + μεψ -Λβ(μσεξ +μηεξ-μσξ-μηε+εξψ+εξγ+με+μξ+φξ+εω+ξ?)),F4=(δ+κ+μ)(μφψ+μφω+μ2(γ+?+φ+?ω+μ+ψ+ω)+μψω+μγ?+μφγ)-Λβ(μ2(1+σε+ηξ-(η+σ))+μω+εψω+φω+?ω+σεω-μηεψ-μφη-μη?+μξγ+φξγ+ξγ?+μηξψ+μηξ?+μφηξ+μσεγ+μφσξ-μφσ-μσξγ-μσω+μεψ+μφ+μ?+φξψ).

3 Threshold Number

By using the next-generation matrix method,we calculate two types of matrices one is transition matrix,and the second is transmission matrix as follows:

After substituting the rotavirus free equilibrium,so we have

The spectral radius of theFG-1is called the threshold number is as follows:

4 Stability Results

Theorem 3:The rotavirus-free equilibrium=(S1,M1,V1,0)is locally asymptotically stable(LAS)whenR0＜1.

Proof:The Jacobian matrix at the rotavirus-free equilibrium is as follows:

λ3+λ2(3μ+2Ψ+γ+?+φ)+λ((μ+γ+Ψ)(μ+Ψ)+(?+φ+μ)(μ+Ψ)+(μ+γ+Ψ)(?+φ+μ)-?Ψ-ωγ)+((μ+γ+Ψ)(?+φ+μ)(μ+Ψ)-?Ψ(μ+Ψ)-ωΨφ-ωγ(?+φ+μ))=0.

By applying Routh-Hurwitz Criterion for 3rdorder,(3μ+2Ψ+γ+?+φ)＞0,((μ+γ+Ψ)(?+φ+μ)(μ+Ψ)-?Ψ(μ+Ψ)-ωΨφ-ωγ(?+φ+μ))＞0,and,

(3μ+2Ψ+γ+?+φ)((μ+γ+Ψ)(μ+Ψ)+(?+φ+μ)(μ+Ψ)+(μ+γ+Ψ)(?+φ+μ)-?Ψ-ωγ)＞((μ+γ+Ψ)(?+φ+μ)(μ+Ψ)-?Ψ(μ+Ψ)-ωΨφ-ωγ(?+φ+μ)),ifR0＜1.So,the rotavirus free equilibria is locally asymptotically stable.

Theorem 4:The rotavirus endemic equilibrium = (S*,M*,V*,I*)is locally asymptotically stable(LAS)whenR0＞1.

Proof:The Jacobian matrix at the rotavirus endemic equilibrium is as follows:

λ4+(μ+ γ + Ψ +(?βI*+(?+φ+μ)) +(β(S*+M*+ξV*) -(δ+k+μ)) +μ)λ3+((μ+γ+Ψ)(?βI*+(?+φ+μ))+(μ+γ+Ψ)(μ+Ψ)+(?βI*+(?+φ+μ))(μ+Ψ)-(μ+γ+Ψ)(β(S*+M*+ξV*)-(Δ+k+μ))-(μ+Ψ)(β(S*+M*+ξV*)-(Δ+k+μ))-(?βI*+(?+φ+μ))(β(S*+M*+ξV*)-(Δ+k+μ))-(-?βM*)?βI*-(-?βM*)βI*-ωγ-DβI*)λ2+((μ+γ+Ψ)(?βI*+(?+φ+μ))(μ+Ψ)-(μ+γ+Ψ)(?βI*+(?+φ+μ))(β(S*+M*+ξV*)-(Δ+k+μ))-(μ+γ+Ψ)(μ+Ψ)(β(S*+M*+ξV*)-(Δ+k+μ))-(μ+γ+Ψ)(-?βM*)?βI*-(?βI*+(?+φ+μ))(μ+Ψ)(β(S*+M*+ξV*)-(Δ+k+μ))+(-ξβI*)ξβI*-(-?βM*)φξβI*-(-?βM*)?βI*(μ+ Ψ) +?(μ+ Ψ)Ψ(β(S*+M*+ξV*) -(δ+k+μ)) -?(μ+ Ψ)Ψ +?(β(S*+M*+ξV*) -(δ+k+μ))Ψ +(-ξβI*)ξβI*Ψ -(-?βM*)γξβI*-(-?βM*)L(μ+Ψ) -ωΨT+ω(?βI*+(?+φ+μ))γ -ωγ(β(S*+M*+ξV*) -(δ+k+μ)) +ω(-ξβI*)βI*-βS*Ψ?βI*-βS*(?βI*+(?+φ+μ))βI*-βS*γξβI*-βS*βI*J)λ+((μ+γ+Ψ)KξβI*+(μ+γ+Ψ)(-?βM*)φξβI*-(μ+γ+Ψ)(?βI*+(?+φ+μ))(μ+Ψ)(β(S*+M*+ξV*)-(Δ+k+μ))-(μ+γ+Ψ)(-?βM*)?βI*(μ+Ψ)+?JΨ(β(S*+M*+ξV*)-(Δ+k+μ))+ωΨI(β(S*+M*+ξV*)-(Δ+k+μ))+ω(?βI*+(?+φ+μ))γ(β(S*+M*+ξV*)-(Δ+k+μ))+ω(?βI*+(?+φ+μ))(-ξβI*)βI*-ωGγ?βI*-βS*ΨIξβI*-βS*Ψ?βI*(μ+Ψ)-βS*(?βI*+(?+φ+μ))γξβI*-βS*(?βI*+(?+φ+μ))βI*(μ+Ψ))=0.

By applying Routh-Hurwitz Criterion for 4thorder,ifR0＞1.So,the rotavirus existing equilibria are locally asymptotically stable.

5 Evolutionary Approaches

To analyze the behavior of the continuous model(6)-(9),we use the evolutionary approaches such as Euler,Runge Kutta,and the nonstandard finite difference is as follows:

5.1 Euler Approach

The system(6)-(9)is selected under Euler computation,as follows:

where,the difference of two consective values of time by h.

5.2 Runge-Kutta Approach

The system(6)-(9)is defined under Runge Kutta computation,as follows:

Stage 1

Stage 2

Stage 3

Stage 4

Final stage

where,the difference of two consective values of time by h andn≥0.

5.3 Nonstandard Finite Difference Approach

The system(6)-(9)is defined under NSFD computation,as follows:

where,the difference of two consective values of time by h.

5.4 Computational Outcomes

By using the parameters’values presented in Tab.1,we approximate the continuous model through the evolutionary approaches as follows:

Table 1: Values of parameters

Figure 2:The graphical behavior of the continuous model at both equilibria of the rotavirus disease(a)Subpopulations for rota VFE at any time t(b)Subpopulations for rota VEE at any time t

Figure 3:Euler computation for the behavior of the continuous model at both equilibria of the rotavirus disease (a)convergent behavior for VFE at h = 0.01 (b)divergent behavior for VFE at h = 3 (c)concurrent behavior for VEE at h=0.01(d)divergent behavior for VEE at h=3

Figure 4:Runge Kutta computation for the behavior of the continuous model at both equilibria of the rotavirus disease(a)At h=0.01 for VFE(b)At h=4 for VFE(c)At h=0.01 for VEE(d)At h=4 for VEE

Figure 5: NSFD computation for the behavior of the continuous model at both equilibria of the rotavirus disease(a)Sub-populations at VEE when h=100(b)Sub-populations at VEE when h=100

5.5 Comparison Unit

Figure 6: (Continued)

Figure 6: Combined graphical behaviors of NSFD with Euler and Runge Kutta computations at different time-step sizes (a)Infected children for VFE at h = 0.01 (b)Infected children for VFE at h=3(c)Infected children for VEE at h=0.01(d)Infected children for VEE at h=4

6 Concluding Remarks

We analyze the dynamics of rotavirus disease via evolutionary approaches.Fig.2 predicts the behavior of the continuous model by using the command build algorithm ODE-45 to solve the system of differential equations.Figs.3a-3d indicates the solution of the model through the Euler approach at both equilibria,such as virus-free and virus existing.Unfortunately,when we try to predict the longterm behavior of the disease in the population, they show negative and unboundedness (exceeding the result from the total population).In the field of epidemiology, these results have no physical relevance.Figs.4a-4d,the well-known approach like Runge Kutta of order four implemented on the continuous model.No doubt,Runge Kutta has the best convergence compared to the Euler approach.But, Runge Kutta has the same issues when predicting long-term behavior, such as negativity,inconsistency,and many more.We construct the nonstandard finite difference approach on the model to overcome these issues.Our proposed construction always gives a positive solution,bounded,and dynamically consistent at any step size.This method predicts the long-term disease behavior and restores the fundamental properties of epidemiological modeling,as shown in Figs.5a-5b.To analyze the efficiency of the approaches,we draw a comparison of existing methods with the proposed system,as shown in Figs.6a-6d.In the end,our proposed policy,nonstandard finite difference,could be more effective in all other types of modeling.Furthermore,flexible tool on behalf of dynamical properties like stability,positivity,and boundedness displays the particular behavior of the continuous model.

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

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

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