|Computers, Materials & Continua |
Mathematical Modelling of Rotavirus Disease Through Efficient Methods
Department of Mathematics, Govt. Maulana Zafar Ali Khan Graduate College Wazirabad, Punjab Higher Education Department (PHED), Lahore, 54000, Pakistan
*Corresponding Author: Ali Raza. Email: firstname.lastname@example.org
Received: 09 January 2022; Accepted: 17 February 2022
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 . In 2020, Ahmed et al. developed a mathematical analysis of epidemic modeling of diarrhea like rotavirus . Lin et al. investigated the causes of diarrheal infectious disease rotavirus (Rotarix and Rotateq) in 2014 . In 2019, Payne et al. presented a rotavirus vaccine to understand their effectiveness and recommend US infant immunization . Omondi et al. presented a mathematical model to explore the co-infection of malaria and rotavirus . In 2009, Effeltrerre et al. studied the dynamics of rotavirus to check the indirect effect of vaccination . 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 . Linhares et al. investigated a rotavirus infection that affects children before five in Brazil by a longitudinal study in 2009 . In 2010, Atchison et al. developed a determined age-structure model to investigate the transmission and effect of rotavirus vaccination in England and Wales . In 2013, Arnold et al. presented a theory of rotavirus-like infections . In 2020, Shuaib et al. developed a mathematical analysis to study the effect of rotavirus on childhood mortality and diarrheal disease . Guzel et al. presented a met-analysis in 2020 regarding awareness programs and preventive policies of rotavirus infection in younger children in Turkey . In 2020, Folorunso et al. investigated oral rotavirus vaccine in developed countries and its strategic preventive . Bibera et al. studied the evolution of the rotavirus vaccine . In 2020, Ilmi et al. investigated the dynamics of the rotavirus epidemic model with the effect of crowding infective individuals . In 2018, Shumetie et al. suggested improving child care to protect their children of a mother from diarrheal morbidity and rotavirus vaccination . Bennett studied the rotavirus epidemic model with the well-known assumption of mathematics . In 2016, Namawejje presented a mathematical model to check the affected regions of rotavirus and illustrate the effects on children and treatment . 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 . 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 population : 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.
The fixed ratios of the model are presented as follows: : 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 : 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 term , 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 conditions .
And , be the feasible region.
Theorem 1: The system results (6)–(9) with given initial conditions are positive for all .
Proof: By letting the Eq. (6),
Similarly, for Eqs. (7)–(9), we have
, as desired.
Theorem 2: The solutions of the system (6)–(9) are bounded at any time and .
Proof: By considering the population function is as follows:
For large t→∞
, as desired.
The system (6)–(9) admits two types of equilibria as follows: rotavirus free equilibrium (RFE-,
rotavirus endemic equilibrium (REE-
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
are the transition and transmission matrices respectively.
The spectral radius of the is called the threshold number is as follows:
4 Stability Results
Theorem 3: The rotavirus-free equilibrium = is locally asymptotically stable (LAS) when < 1.
Proof: The Jacobian matrix at the rotavirus-free equilibrium is as follows:
, if .
By applying Routh-Hurwitz Criterion for 3rd order, and,
, if So, the rotavirus free equilibria is locally asymptotically stable.
Theorem 4: The rotavirus endemic equilibrium = (,,,) is locally asymptotically stable (LAS) when > 1.
Proof: The Jacobian matrix at the rotavirus endemic equilibrium is as follows:
By applying Routh-Hurwitz Criterion for 4th order, if 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
5.2 Runge-Kutta Approach
The system (6)–(9) is defined under Runge Kutta computation, as follows:
where, the difference of two consective values of time by and
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 .
5.4 Computational Outcomes
By using the parameters’ values presented in Tab. 1, we approximate the continuous model through the evolutionary approaches as follows:
5.5 Comparison Unit
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 long-term 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.
|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.|