{"title":"An unconditional boundary and dynamics preserving scheme for the stochastic epidemic model","authors":"Ruishu Liu, Xiaojie Wang, Lei Dai","doi":"10.1007/s10092-024-00606-z","DOIUrl":null,"url":null,"abstract":"<p>In the present article, we construct a logarithm transformation based Milstein-type method for the stochastic susceptible-infected-susceptible (SIS) epidemic model evolving in the domain (0, <i>N</i>). The new scheme is explicit and unconditionally boundary and dynamics preserving, when used to solve the stochastic SIS epidemic model. Also, it is proved that the scheme has a strong convergence rate of order one. Different from existing time discretization schemes, the newly proposed scheme for any time step size <span>\\(h>0\\)</span>, not only produces numerical approximations living in the entire domain (0, <i>N</i>), but also unconditionally reproduces the extinction and persistence behavior of the original model, with no additional requirements imposed on the model parameters. Numerical experiments are presented to verify our theoretical findings.</p>","PeriodicalId":9522,"journal":{"name":"Calcolo","volume":"13 1","pages":""},"PeriodicalIF":1.4000,"publicationDate":"2024-08-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Calcolo","FirstCategoryId":"100","ListUrlMain":"https://doi.org/10.1007/s10092-024-00606-z","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS","Score":null,"Total":0}
引用次数: 0
Abstract
In the present article, we construct a logarithm transformation based Milstein-type method for the stochastic susceptible-infected-susceptible (SIS) epidemic model evolving in the domain (0, N). The new scheme is explicit and unconditionally boundary and dynamics preserving, when used to solve the stochastic SIS epidemic model. Also, it is proved that the scheme has a strong convergence rate of order one. Different from existing time discretization schemes, the newly proposed scheme for any time step size \(h>0\), not only produces numerical approximations living in the entire domain (0, N), but also unconditionally reproduces the extinction and persistence behavior of the original model, with no additional requirements imposed on the model parameters. Numerical experiments are presented to verify our theoretical findings.
期刊介绍:
Calcolo is a quarterly of the Italian National Research Council, under the direction of the Institute for Informatics and Telematics in Pisa. Calcolo publishes original contributions in English on Numerical Analysis and its Applications, and on the Theory of Computation.
The main focus of the journal is on Numerical Linear Algebra, Approximation Theory and its Applications, Numerical Solution of Differential and Integral Equations, Computational Complexity, Algorithmics, Mathematical Aspects of Computer Science, Optimization Theory.
Expository papers will also appear from time to time as an introduction to emerging topics in one of the above mentioned fields. There will be a "Report" section, with abstracts of PhD Theses, news and reports from conferences and book reviews. All submissions will be carefully refereed.