{"title":"Mitigating negative effects of EBHSV-infected eastern cottontail invasion in Italy using Z-type control on a four-population system","authors":"Francesca Acotto , Fabiana Camattari , Ezio Venturino","doi":"10.1016/j.matcom.2025.01.024","DOIUrl":"10.1016/j.matcom.2025.01.024","url":null,"abstract":"<div><div>The introduction in Italy of eastern cottontails (<em>Sylvilagus floridanus</em>) for hunting purposes has influenced the local predator–prey dynamics of the red foxes (<em>Vulpes vulpes</em>) and native European hares (<em>Lepus europaeus</em>). Although no direct competition seems to occur between the two lagomorphs, the cottontail invasion damages the indigenous hare population. Indeed, the invasive lagomorphs cause hyperpredation of red foxes on native hares and are also carriers of viruses and parasites. This paper focuses on the situation in which EBHSV-infected eastern cottontails are introduced in a region of virus-free European hares. To avoid the extinction of native lagomorphs and to contain the invasive ones, we look at two possible biological control actions using the Z-type control method on a four-population reference system. In particular, we consider an indirect control of the invasive prey acting on predators and a combination of this indirect control with direct control on native prey. The corresponding Z-controlled models are investigated analytically and numerically. In both cases, the Z-type control significantly reduces the number of equilibria and the convergence of the cottontail population to the desired state is ensured. The hare survival, instead, is guaranteed only in the second case. Overall, mathematically speaking, the second Z-type control action seems the best solution. However, in the cases in which the indirect control on cottontails allows the native prey survival, this control may be preferable since it involves only one control function and seems more practicable. In any case, in the choice of control action, ecosystem managers need to consider each specific situation, taking into account various elements from a biological and practical point of view.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"233 ","pages":"Pages 117-136"},"PeriodicalIF":4.4,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143226661","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Excitability and synchronization of vanadium dioxide memristor-inspired neurons","authors":"Yan Shao , Fuqiang Wu , Qingyun Wang","doi":"10.1016/j.matcom.2025.01.022","DOIUrl":"10.1016/j.matcom.2025.01.022","url":null,"abstract":"<div><div>Neuromorphic devices play a significant role in exploiting the dynamical analogy between various physical circuits and neuronal systems. Mott memristive systems made from vanadium dioxide (VO<sub>2</sub>) are candidates for neuromorphic computations because of their ability to better reproduce neuron-like functions/behaviors and operate at low transition energy. In this paper, we revisit a prior work on the VO<sub>2</sub> memristor-inspired neuron to not only identify the exact biophysical mechanisms of this system’s features, but also reproduce the more extensive neuron-like dynamical behaviors as two types of excitability and spiking by combining one-parameter bifurcations with two-parameter panels. Based on the Lyapunov stability theorem, a novel criterion for exponential synchronization is acquired in coupled VO<sub>2</sub> memristor-inspired neurons by introducing two controllers. Meanwhile, it is demonstrated via numerical simulations and hardware circuits. Results provide the intersection of electronic physics and theoretical neuroscience from the nonlinear dynamics point of view.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"233 ","pages":"Pages 99-116"},"PeriodicalIF":4.4,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143226658","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Split-step ϑ integrator for generalized stochastic Volterra integro-differential equations","authors":"Hassan Ranjbar","doi":"10.1016/j.matcom.2025.01.021","DOIUrl":"10.1016/j.matcom.2025.01.021","url":null,"abstract":"<div><div>This study investigates the theoretical analysis and numerical approximation of generalized stochastic Volterra integro-differential equations. First, we examine the existence, uniqueness, boundedness and Hölder continuity of the exact solutions for generalized SVIDEs. To numerically solve, the split-step <span><math><mi>ϑ</mi></math></span> integrator is proposed. We then demonstrate the boundedness of the numerical solution. Further, it has been shown that the scheme is strongly convergent with order half. Numerical experiments are carried out to support our findings.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"233 ","pages":"Pages 165-186"},"PeriodicalIF":4.4,"publicationDate":"2025-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143226659","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Multi-parameter bifurcations in a discrete Ricker-type predator–prey model with prey immigration","authors":"Karima Mokni, Hajar Mouhsine, Mohamed Ch-Chaoui","doi":"10.1016/j.matcom.2025.01.020","DOIUrl":"10.1016/j.matcom.2025.01.020","url":null,"abstract":"<div><div>This study examines a discrete-time prey–predator model featuring a Ricker-type growth function and immigration effects to uncover the dynamics shaping ecosystem stability. Through detailed bifurcation analysis, we identify codimension-one bifurcations, including transcritical, period-doubling, and Neimark–Sacker bifurcations, as well as codimension-two bifurcations involving 1:2, 1:3, and 1:4 resonances. Our results reveal that low immigration rates stabilize the system, ensuring predictable population dynamics, while exceeding critical thresholds induces complex behaviors, such as periodic oscillations and chaos. We numerically analyze the dynamics associated with 1:2, 1:3, and 1:4 resonances, utilizing two-parameter bifurcation diagrams and basins of attraction to illustrate the transitions and stability boundaries. These findings highlight the dual role of immigration in stabilizing and destabilizing ecosystems, offering valuable insights for ecological modeling, management, and conservation strategies.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"233 ","pages":"Pages 39-59"},"PeriodicalIF":4.4,"publicationDate":"2025-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143226047","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Stability analysis and optimal control of tumour-immune interaction problem using fractional order derivative","authors":"Tarekegn Dinku , Boka Kumsa , Jyotirmoy Rana , Aiyappan Srinivasan","doi":"10.1016/j.matcom.2024.12.028","DOIUrl":"10.1016/j.matcom.2024.12.028","url":null,"abstract":"<div><div>In this study, we propose a tumour-immune interaction model using Caputo–Fabrizio fractional order derivative. The conditions for the well-posedness of the solution are examined. The stability of the endemic equilibrium point is derived and its stability is proved using Routh–Hurwitz criteria. The solution is approximated using a shifted Legendre polynomial at Gauss–Legendre collocation points, which is compared with the numerical results of the Adams–Bashforth scheme in the interval <span><math><mrow><mo>[</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>]</mo></mrow></math></span>. We have also proposed a fractional optimal control problem and proved the necessary optimality conditions. The optimal system is solved using the forward–backward sweep method (FBSM) with the Adams–Bashforth predictor–corrector numerical method. We have demonstrated that the antigenicity of tumours plays a crucial role in activating immune cells, suggesting that enhancing tumour antigenicity could improve immunotherapeutic outcomes. The effects of fractional-order derivatives and the proliferation rate of the Michaelis–Menten term are observed. Moreover, the impact of other model parameters on the system is highlighted through numerical results. Finally, the reduction in tumour cells and the increase of active immune cells are demonstrated in the presence of optimal control.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"233 ","pages":"Pages 187-207"},"PeriodicalIF":4.4,"publicationDate":"2025-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143226660","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"News of IMACS","authors":"","doi":"10.1016/S0378-4754(25)00027-8","DOIUrl":"10.1016/S0378-4754(25)00027-8","url":null,"abstract":"","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"231 ","pages":"Page 352"},"PeriodicalIF":4.4,"publicationDate":"2025-01-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143133917","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Influence of price elasticity of demand on monopoly games under different returns to scale","authors":"Xiaoliang Li , Jing Yang , Ally Quan Zhang","doi":"10.1016/j.matcom.2025.01.017","DOIUrl":"10.1016/j.matcom.2025.01.017","url":null,"abstract":"<div><div>This paper examines a monopoly market featured by a general isoelastic demand function. With a quadratic cost function for the monopolist, we explore how the price elasticity of demand influences monopolistic behavior under different (decreasing, constant, and increasing) returns to scale. The combination of the general isoelastic demand and quadratic cost functions leads to a transcendental equilibrium equation, making closed-form solutions unattainable. To address this challenge, we develop an innovative approach that leverages the specific structure of marginal revenue and cost to conduct a comprehensive comparative static and stability analysis. Additionally, we introduce two dynamic models based on distinct adjustment mechanisms: gradient and local monopolistic approximation (LMA). Our findings reveal that the LMA model is more stable in both parameter and state spaces compared to the gradient model. Notably, we prove that the unique non-vanishing equilibrium of the LMA model is globally asymptotically stable.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"233 ","pages":"Pages 75-98"},"PeriodicalIF":4.4,"publicationDate":"2025-01-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143226663","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"IMACS Calendar of Events","authors":"","doi":"10.1016/S0378-4754(25)00028-X","DOIUrl":"10.1016/S0378-4754(25)00028-X","url":null,"abstract":"","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"231 ","pages":"Page 353"},"PeriodicalIF":4.4,"publicationDate":"2025-01-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143133918","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Alberto M. Bersani , Alessandro Borri , Giovanna Tomassetti , Pierluigi Vellucci
{"title":"Asymptotic expansions beyond the tQSSA for the double phosphorylation mechanism","authors":"Alberto M. Bersani , Alessandro Borri , Giovanna Tomassetti , Pierluigi Vellucci","doi":"10.1016/j.matcom.2025.01.019","DOIUrl":"10.1016/j.matcom.2025.01.019","url":null,"abstract":"<div><div>In this paper we study the asymptotic properties of the mathematical model of the double phosphorylation (dephosphorylation) enzymatic reaction, or futile cycle. Starting from the total quasi-steady state approximation (tQSSA), and applying singular perturbation techniques, we determine the inner and outer solutions and the corresponding matched expansions, up to the first order, in terms of an appropriate perturbation parameter <span><math><mi>ϵ</mi></math></span> (related to the kinetic constants and initial conditions of the model). Some numerical results are discussed.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"233 ","pages":"Pages 137-164"},"PeriodicalIF":4.4,"publicationDate":"2025-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143226662","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Néstor García-Chan , Lino J. Alvarez-Vázquez , Aurea Martínez , Miguel E. Vázquez-Méndez
{"title":"A nonconservative macroscopic traffic flow model in a two-dimensional urban-porous city","authors":"Néstor García-Chan , Lino J. Alvarez-Vázquez , Aurea Martínez , Miguel E. Vázquez-Méndez","doi":"10.1016/j.matcom.2025.01.016","DOIUrl":"10.1016/j.matcom.2025.01.016","url":null,"abstract":"<div><div>In this paper we propose a novel traffic flow model based on understanding the city as a porous media, this is, streets and building-blocks characterizing the urban landscape are seen now as the fluid-phase and the solid-phase of a porous media, respectively. Moreover, based in the interchange of mass in the porous media models, we can model the interchange of cars between streets and off-street parking-spaces. Therefore, our model is not a standard conservation law, being formulated as the coupling of a non-stationary convection–diffusion–reaction PDE with a Darcy–Brinkman–Forchheimer PDE system. To solve this model, the classical Galerkin <span><math><msub><mrow><mi>P</mi></mrow><mrow><mn>1</mn></mrow></msub></math></span> finite element method combined with an explicit time marching scheme of strong stability preserving type was enough to stabilize our numerical solutions. Numerical experiences on an urban-porous domain inspired by the city of Guadalajara (Mexico) allow us to simulate the influence of the porosity terms on the traffic speed, the traffic flow at rush-valley hours, and traffic congestion due to the lack of parking spaces.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"233 ","pages":"Pages 60-74"},"PeriodicalIF":4.4,"publicationDate":"2025-01-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143226664","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}