“…It has been determined for some types of orthogonal, biorthogonal, and interpolate polynomials, such as flatlet oblique multiwavelets [25][26][27][28], B-spline cardinal functions [10,29], Chebyshev cardinal functions [30][31][32], Chebyshev polynomials [33], and Legendre polynomials [34] and some the other operational matrix methods [35][36][37][38][39][40][41]. In addition, recently, the creation of an operational matrix of the stochastic integral has become a popular approach for solving stochastic problems, and many researchers have tried to construct a simple and reliable operational matrix (for more details, see [7,[10][11][12][13][14][15][16][17][18][19][20][21]42]).…”