Construction of operational matrices based on linear cardinal B-spline functions for solving fractional stochastic integro-differential equation

The main purpose of this paper is to develop a new method based on operational matrices of the linear cardinal B-spline (LCB-S) functions to numerically solve of the fractional stochastic integro-differential (FSI-D) equations. To reach this aim, LCB-S functions are introduced and their properties are considered, briefly. Then, the operational matrices based on LCB-S functions are constructed, for the first time, including the fractional Riemann-Liouville integral operational matrix, the stochastic integral operational matrix, and the integer integral operational matrix. The main characteristic of the new scheme is to convert the FSI-D equation into a linear system of algebraic equations which can be easily solved by applying a suitable method. Also, the convergence analysis and error estimate of the proposed method are studied and an upper bound of error is obtained. Numerical experiments are provided to show the potential and efficiency of the new method. Finally, some numerical results, for various values of perturbation in the parameters of the main problem are presented which can indicate the stability of the suggested method.