Computational scheme for solving nonlinear fractional stochastic differential equations with delay

This paper studies the numerical solution of fractional stochastic delay differential equations driven by Brownian motion. The proposed algorithm is based on linear B-spline interpolation. The convergence and the numerical performance of the method are analyzed. The technique is adopted for determining the statistical indicators of stochastic responses of fractional Langevin and Mackey-Glass models with stochastic excitations.