This paper introduces a novel numerical technique for solving fractional stochastic differential equations with neutral delays. The method employs a stepwise collocation scheme with Jacobi poly-fractonomials to consider unknown stochastic processes. For this purpose, the delay differential equations are transformed into augmented ones without delays. This transformation makes it possible to use a collocation scheme improved with Jacobi poly-fractonomials to solve the changed equations repeatedly. At each iteration, a system of nonlinear equations is generated. Next, the convergence properties of the proposed method are rigorously analyzed. Afterward, the practical utility of the proposed numerical technique is validated through a series of test examples. These examples illustrate the method’s capability to produce accurate and efficient solutions.
