Balajied Me Syrti, Anuradha Devi, Ankur Jyoti Kashyap
{"title":"Analysis of Stability, Sensitivity Index and Hopf Bifurcation of Eco-Epidemiological SIR Model under Pesticide Application","authors":"Balajied Me Syrti, Anuradha Devi, Ankur Jyoti Kashyap","doi":"10.5614/cbms.2023.6.2.4","DOIUrl":"https://doi.org/10.5614/cbms.2023.6.2.4","url":null,"abstract":"In this paper, a deterministic SIR plant mathematical model is proposed and analysed with the application of pesticides as a control measure. The primary purpose of this model is to study the role of pesticides in controlling disease prevalence in plant populations. The total plant population is subdivided into three categories: susceptible, infected, and recovered. Pesticides are considered to be applied to both susceptible and infected populations to prevent the spread of infection to unaffected plant populations. It is considered that plant populations can be recovered only through the use of pesticides. To ensure the biological validity and well-defined nature of the model, the positivity, boundedness, uniqueness and existence of solutions are analysed. The basic reproduction number (R0) of the infection is determined and observed that the disease-free equilibrium state is locally asymptotically stable whenever (R0) is less than unity and unstable otherwise. The sensitivity analysis of the basic reproduction number is carried out, and it is observed that the value of R0 decreases as the value of the death rate and the recovery rate of plants increases. Moreover, it is revealed that above a critical parameter value of the infective induce rate, the population starts oscillating periodically, and the endemic equilibrium state becomes unstable. Finally, numerical simulations are conducted in MATLAB software to compare the analytical findings. Overall, the results obtained from this analysis are both novel and significant, making them an intriguing and potentially valuable contribution to the field of theoretical ecology.","PeriodicalId":33129,"journal":{"name":"Communication in Biomathematical Sciences","volume":"108 18","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138959159","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Rolly Czar Joseph Castillo, Victoria May Mendoza, Jose Ernie Lope, Renier Mendoza
{"title":"Modeling Infectious Disease Trend using Sobolev Polynomials","authors":"Rolly Czar Joseph Castillo, Victoria May Mendoza, Jose Ernie Lope, Renier Mendoza","doi":"10.5614/cbms.2023.6.2.2","DOIUrl":"https://doi.org/10.5614/cbms.2023.6.2.2","url":null,"abstract":"Trend analysis plays an important role in infectious disease control. An analysis of the underlying trend in the number of cases or the mortality of a particular disease allows one to characterize its growth. Trend analysis may also be used to evaluate the effectiveness of an intervention to control the spread of an infectious disease. However, trends are often not readily observable because of noise in data that is commonly caused by random factors, short-term repeated patterns, or measurement error. In this paper, a smoothing technique that generalizes the Whittaker-Henderson method to infinite dimension and whose solution is represented by a polynomial is applied to extract the underlying trend in infectious disease data. The solution is obtained by projecting the problem to a finite-dimensional space using an orthonormal Sobolev polynomial basis obtained from Gram-Schmidt orthogonalization procedure and a smoothing parameter computed using the Philippine Eagle Optimization Algorithm, which is more efficient and consistent than a hybrid model used in earlier work. Because the trend is represented by the polynomial solution, extreme points, concavity, and periods when infectious disease cases are increasing or decreasing can be easily determined. Moreover, one can easily generate forecast of cases using the polynomial solution. This approach is applied in the analysis of trends, and in forecasting cases of different infectious diseases.","PeriodicalId":33129,"journal":{"name":"Communication in Biomathematical Sciences","volume":" 22","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138962353","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Blood Glucose Control on Diabetic Patient Type I using Sliding Mode Adaptive Control","authors":"Aminatus Sa'adah, Prihantini Prihantini","doi":"10.5614/cbms.2023.6.2.1","DOIUrl":"https://doi.org/10.5614/cbms.2023.6.2.1","url":null,"abstract":"Diabetes is a metabolic disorder due to insufficient insulin synthesis or inadequate insulin sensitivity. The Bergman’s minimal model describes the dynamics of blood glucose levels in type 1 diabetics. The model has control inputs in the form of insulin injections and covers external disturbance factors in the form of meal disturbances. This research developed a control design using an sliding mode adaptive control to reduce blood glucose levels in hyperglycemic patients and keep it within normal glucose levels. Sliding mode adaptive control is an adaptive controller updates the model based on measured performance while in operation. A numerical simulation of the proposed controller is carried out by giving eating disorders three times, namely at breakfast, lunch, and dinner. Based on the numerical simulation, to lower the high blood glucose in the hyperglicemic patient, the insulin injection should be given starting at 30 minutes before breakfast for the next four hour, with a maximal dose of injection is 13 mU/min. It can decrease the high blood pressure until 54.83%.","PeriodicalId":33129,"journal":{"name":"Communication in Biomathematical Sciences","volume":" 8","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138962511","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Dynamical Behaviour of Pro- and Anti-Inflammatory Cytokines during Pathogenesis of Atherosclerosis","authors":"Asish Adak, Arpita Devi, Praveen Kumar Gupta","doi":"10.5614/cbms.2023.6.2.3","DOIUrl":"https://doi.org/10.5614/cbms.2023.6.2.3","url":null,"abstract":"The role of anti-cytokines in atherosclerosis is to reduce inflammation in the intima. In some situations, certain anti-inflammatory cytokines like TGF-beta and IL-6 have shown the characteristics like a pro-inflammatory cytokines, which are showing different natures. In this study, a dynamical atherosclerosis model is proposed in the form of reaction-diffusion equation with consideration of immune cells, pro-inflammatory cytokines, and anti-inflammatory cytokines. The existence and uniqueness of the solutions are discussed for the proposed reaction dynamical system. The three equilibrium points, non-inflammatory, chronic, and coexistence, and their local stability are also determined for the model. Bellman and Cooke’s theorem is applied to illustrate the global stability at the coexistence equilibrium point. The effects of pro- and anti-inflammatory cytokines have also been discussed. The analytical and numerical studies evidently indicate that inflammation behaves differently when a certain number of anti-inflammatory cytokines behave like pro-inflammatory cytokines. The numerical simulations are demonstrated for different impacts of the reduction rate of macrophages due to the presence of anti-inflammatory cytokines, inhibition time, and the portion of anti-inflammatory cytokines behaving like pro-inflammatory cytokines through graphically. The results of this study suggest that chronic inflammation of the disease is likely to persist when a high concentration of ox-LDL and moderate concentration of cytokines are present in the intima. Coexistence inflammation is characterized by a high concentration of ox-LDL, moderate concentration of pro-inflammatory and high concentration of anti-cytokines; whereas a non-inflammatory condition would persevere if a low concentration of ox-LDL has been present in the intima.","PeriodicalId":33129,"journal":{"name":"Communication in Biomathematical Sciences","volume":" 493","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138960439","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A Numerical Study of Codimension-Two Bifurcations of an SIR-Type Model for COVID-19 and Their Epidemiological Implications","authors":"Livia Owen, Jonathan Hoseana, Benny Yong","doi":"10.5614/cbms.2023.6.2.6","DOIUrl":"https://doi.org/10.5614/cbms.2023.6.2.6","url":null,"abstract":"We study the codimension-two bifurcations exhibited by a recently-developed SIR-type mathematical model for the spread of COVID-19, as its two main parameters —the susceptible individuals’ cautiousness level and the hospitals’ bed-occupancy rate— vary over their domains. We use AUTO to generate the model’s bifurcation diagrams near the relevant bifurcation points: two Bogdanov-Takens points and two generalised Hopf points, as well as a number of phase portraits describing the model’s orbital behaviours for various pairs of parameter values near each bifurcation point. The analysis shows that, when a backward bifurcation occurs at the basic reproduction threshold, the transition of the model’s asymptotic behaviour from endemic to disease-free takes place via an unexpectedly complex sequence of topological changes, involving the births and disappearances of not only equilibria but also limit cycles and homoclinic orbits. Epidemiologically, the analysis confirms the importance of a proper control of the values of the aforementioned parameters for a successful eradication of COVID-19. We recommend a number of strategies by which such a control may be achieved.","PeriodicalId":33129,"journal":{"name":"Communication in Biomathematical Sciences","volume":" 48","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138963416","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Afifurrahman Afifurrahman, Mohd Hafiz Mohd, Farah Aini Abdullah
{"title":"The Interplay of Common Noise and Finite Pulses on Biological Neurons","authors":"Afifurrahman Afifurrahman, Mohd Hafiz Mohd, Farah Aini Abdullah","doi":"10.5614/cbms.2023.6.2.5","DOIUrl":"https://doi.org/10.5614/cbms.2023.6.2.5","url":null,"abstract":"The response of neurons is highly sensitive to the stimulus. The stimulus can be associated with a direct injection in vitro experimentation (e.g., time dependent and independent inputs); or post-synaptic potentials resulting from the interaction of many neurons. A typical incoming stimulus resembles a noise which in principle can be described as a random variable. In computational neuroscience, the noise has been extensively studied for different setups. In this study, we investigate the effect of noisy inputs in a minimal network of two identical leaky integrate-and-fire (LIF) neurons interacting with finite pulses. In particular, we consider a Gaussian white noise as a standard function for stochastic modelling of neurons, while taking into account the pulse width as an elementary component for the signal transmission. By exploring the role of noise and finite pulses, the two neurons show a synchronous spiking behaviour characterized by fluctuations in the interspike intervals. Above some critical values the synchronous regime collapses onto asynchronous dynamics. The abrupt change in such dynamics is accompanied by a hysteresis, i.e., the coexistence of synchronous and asynchronous firing behaviour.","PeriodicalId":33129,"journal":{"name":"Communication in Biomathematical Sciences","volume":" 12","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138961491","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Nadia Nadia, Afdol Zikri, Sila Rizqina, K. Sukandar, M. Fakhruddin, C. J. Tay, N. Nuraini
{"title":"Data-Driven Generating Operator in SEIRV Model for COVID-19 Transmission","authors":"Nadia Nadia, Afdol Zikri, Sila Rizqina, K. Sukandar, M. Fakhruddin, C. J. Tay, N. Nuraini","doi":"10.5614/cbms.2023.6.1.6","DOIUrl":"https://doi.org/10.5614/cbms.2023.6.1.6","url":null,"abstract":"The COVID-19 (SARS-CoV-2) vaccine has been extensively implemented through large-scale programs in numerous countries as a preventive measure against the resurgence of COVID-19 cases. In line with this vaccination effort, the Indonesian government has successfully inoculated over 74% of its population. Nevertheless, a significant decline in the duration of vaccine-induced immunity has raised concerns regarding the necessity of additional inoculations, such as booster shots. Prior to proceeding with further inoculation measures, it is imperative for the government to assess the existing level of herd immunity, specifically determining whether it has reached the desired threshold of 70%. To shed light on this matter, our objective is to ascertain the herd immunity level following the initial and subsequent vaccination programs, while also proposing an optimal timeframe for conducting additional inoculations. This study utilizes COVID-19 data from Jakarta and employs the SEIRV model, which integrates time-dependent parameters and incorporates an additional compartment to represent the vaccinated population. By formulating a dynamic generator based on the cumulative cases function, we are able to comprehensively evaluate the analytical and numerical aspects of all state dynamics. Simulation results reveal that the number of individuals protected by the vaccine increases following the vaccination program; however, this number subsequently declines due to the waning effect of the vaccine. Our estimates indicate that the vaccination program in Jakarta has achieved herd immunity levels exceeding 70% from October 2021 to February 2022, thus underscoring the necessity of rolling out further inoculations no later than February 2022.","PeriodicalId":33129,"journal":{"name":"Communication in Biomathematical Sciences","volume":" ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-07-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45275565","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
B. Handari, D. Aldila, Evllyn Tamalia, S. Khoshnaw, M. Shahzad
{"title":"Assessing The Impact of Medical Treatment and Fumigation on The Superinfection of Malaria: A Study of Sensitivity Analysis","authors":"B. Handari, D. Aldila, Evllyn Tamalia, S. Khoshnaw, M. Shahzad","doi":"10.5614/cbms.2023.6.1.5","DOIUrl":"https://doi.org/10.5614/cbms.2023.6.1.5","url":null,"abstract":"Malaria is a disease caused by the parasite Plasmodium, transmitted by the bite of an infected female Anopheles. In general, five species of Plasmodium that can cause malaria. Of the five species, Plasmodium falciparum and Plasmodium vivax are two species of Plasmodium that can allow malaria superinfection in the human body. Typically, the popular intervention for malaria eradication is the use of fumigation to control the vector population and provide good medical services for malaria patients. Here in this article, we formulate a mathematical model based on a host-vector interaction. Our model considering two types of plasmodium in the infection process and the use of medical treatment and fumigation for the eradication program. Our analytical result succeeds in proving the existence of all equilibrium points and how their existence and local stability criteria depend not only on the control reproduction number but also in the invasive reproduction number. This invasive reproduction number represent how one plasmodium can dominate other plasmodium. Our sensitivity analysis shows that fumigation is the most influential parameter in determining all control reproduction numbers. Furthermore, we find that the order in which numerous intervention measures are taken will be very crucial to determine the level of success of our malaria eradication program.","PeriodicalId":33129,"journal":{"name":"Communication in Biomathematical Sciences","volume":" ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-07-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48700282","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Mathematical Model and Study of Viral Hepatitis among Population in Afghanistan","authors":"A. Hasmani, Burhanuddin Safi, A. Das","doi":"10.5614/cbms.2023.6.1.4","DOIUrl":"https://doi.org/10.5614/cbms.2023.6.1.4","url":null,"abstract":"Despite availability of strategies against viral hepatitis, it is still a serious disease, which millions of people are already infected with, hence it yet needs to be focused on. As an attempt, we formulated a single mathematical model describing behaviour of all strains of viral hepatitis, presented in the literature. The basic reproduction number(R_0) at disease free equilibrium point is computed, feasible region has been determined. For local stability of the model, R_0 has been taken into account and for global stability of the model Lyapunov method is followed. The model is then applied to the data available for Afghanistan for the year 2020. Based on the data, values of the parameters are estimated, using Minimum Mean Absolute Error (MAE) method. Numerical simulation is performed to support the model and then the results are plotted and represented graphically. One-at-a-time sensitivity analysis (OAT) method is used for sensitivity analysis and involved parameters have been examined for the propose of sensitivity analysis, it indicated that infection rates of acute and chronic states of viral hepatitis are the most sensitive and critical parameters. It has been observed that large number of populations can become infected followed by small increment of infection rates. It has also been noticed that, entire population of Afghanistan could become infected, if no prevention measures were taken. The model presented in this paper is useful for forecasting outbreak by viral hepatitis and it can further be modified by including prevention measures.","PeriodicalId":33129,"journal":{"name":"Communication in Biomathematical Sciences","volume":" ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-07-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41714963","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Maya Rayungsari, A. Suryanto, W. M. Kusumawinahyu, I. Darti
{"title":"Nonstandard Numerical Scheme for a Predator-Prey Model Involving Predator Cannibalism and Refuge","authors":"Maya Rayungsari, A. Suryanto, W. M. Kusumawinahyu, I. Darti","doi":"10.5614/cbms.2023.6.1.2","DOIUrl":"https://doi.org/10.5614/cbms.2023.6.1.2","url":null,"abstract":"In this study, we implement a Nonstandard Finite Difference (NSFD) scheme for a predator-prey model involving cannibalism and refuge in predator. The scheme which is considered as a discrete dynamical system is analyzed. The performed analysis includes the determination of equilibrium point and its local stability. The system has four equilibrium points, namely the origin, the prey extinction point, the predator extinction point, and the coexistence point, which have exactly the same form and existence conditions as those in continuous system. The local stability of each first three equilibrium points is consistent with the one in continuous system. The stability of the coexistence point depends on the integration time step size. Nevertheless, the NSFD scheme allows us to choose the integration time step size for the solution to converge to a feasible point more flexible than the Euler and 4th order Runge-Kutta schemes. These are shown via numerical simulations.","PeriodicalId":33129,"journal":{"name":"Communication in Biomathematical Sciences","volume":" ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45012533","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}