In this paper, a numerical method based on modified hat functions (MHFs) is investigated to find an approximate solution for a fractional-order system of stochastic integro-differential equations. The fractional and stochastic operational matrices of integration of these functions are employed to present the numerical approach. By using these operational matrices and properties of MHFs, the considered problem is transformed into a system of algebraic equations which can be easily solved by an iterative method. Also, error analysis of the proposed method is discussed. At the end, the accuracy and effectiveness of this approach are studied by some numerical examples.
