Numerical solutions of stochastic Lotka-Volterra equations via operational matrices