Numerical solution of stochastic fractional differential equations is proposed.
The spatial approximation method is based on Galerkin and the temporal approximation is based on implicit Euler scheme.
The order on convergence of implicit method for SFDEs is 1.