to estimate the solution of two-dimensional (2D) stochastic integral equations on hypercube areas. The main advantage of the suggested approach is that this algorithm can be easily implemented to estimate the solution of multidimensional stochastic integral equations defined on irregular domains. Also, it is established that the convergence order is proportional to hl X,D, where hX,D denotes fill distance parameter. Finally, to reveal accuracy, efficiency and applicability of our scheme two test problems are included.
