In this work, we present a collocation-based numerical approach for handling stochastic
delay differential equations. We approximate the solution function, and after that,
we carry out integrations between the predefined collocation points to achieve a mapping
from the delayed state to the present state. We build the first and second moment
mapping matrices based on the mapping, and we utilize the matrices to approximate
the stationary first and second moments and their stability. Numerical studies of
a first and second-order stochastic delay differential equation show the convergence
and time complexity of the stochastic collocation method. The last section covers
the issues and possible further improvements of the method.
