Abstract
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.