- Journal Home
- Volume 22 - 2025
- Volume 21 - 2024
- Volume 20 - 2023
- Volume 19 - 2022
- Volume 18 - 2021
- Volume 17 - 2020
- Volume 16 - 2019
- Volume 15 - 2018
- Volume 14 - 2017
- Volume 13 - 2016
- Volume 12 - 2015
- Volume 11 - 2014
- Volume 10 - 2013
- Volume 9 - 2012
- Volume 8 - 2011
- Volume 7 - 2010
- Volume 6 - 2009
- Volume 5 - 2008
- Volume 4 - 2007
- Volume 3 - 2006
- Volume 2 - 2005
- Volume 1 - 2004
Int. J. Numer. Anal. Mod., 22 (2025), pp. 755-776.
Published online: 2025-08
Cited by
- BibTex
- RIS
- TXT
The main goal of this paper is to develop an improved stochastic $\theta$-scheme as a numerical method for stochastic Volterra integral equations (SVIEs) with double weakly singular kernels and demonstrate that the stability of the proposed scheme is affected by the kernel parameters. To overcome the low computational efficiency of the stochastic $\theta$-scheme, we employed the sum-of-exponentials (SOE) approximation. Then, the mean square stability of the proposed scheme with respect to a convolution test equation is studied. Additionally, based on the stability conditions and the explicit structure of the stability matrices, analytical and numerical stability regions are plotted and compared with the split-step $\theta$-method and the $\theta$-Milstein method. The results confirm that our approach aligns significantly with the expected physical interpretations.
}, issn = {2617-8710}, doi = {https://doi.org/10.4208/ijnam2025-1033}, url = {http://global-sci.org/intro/article_detail/ijnam/24292.html} }The main goal of this paper is to develop an improved stochastic $\theta$-scheme as a numerical method for stochastic Volterra integral equations (SVIEs) with double weakly singular kernels and demonstrate that the stability of the proposed scheme is affected by the kernel parameters. To overcome the low computational efficiency of the stochastic $\theta$-scheme, we employed the sum-of-exponentials (SOE) approximation. Then, the mean square stability of the proposed scheme with respect to a convolution test equation is studied. Additionally, based on the stability conditions and the explicit structure of the stability matrices, analytical and numerical stability regions are plotted and compared with the split-step $\theta$-method and the $\theta$-Milstein method. The results confirm that our approach aligns significantly with the expected physical interpretations.