this paper introduces a numerical method for solving the vasicek model by using a stochastic operational matrix based on the triangular functions (tfs) in combination with the collocation method. the method is stated by using conversion the vasicek model to a stochastic nonlinear system of $2m+2$ equations and $2m+2$ unknowns. finally, the error analysis and some numerical examples are provided...