Numerical solution of stochastic Volterra integral equations via a stochastic operational matrix based on block pulse functions.