In this article, the Laplace transformation method in connection with the variational iteration method will be used to solve approximately fuzzy random ordinary differential equations. After that, the sequence of approximated closed form iterated solutions is derived based on the general Lagrange multiplier evaluated using the well-known convolution theorem of the Laplace transformation method. In addition, two examples are given and solved to illustrate the reliability, efficiency and applicability of the proposed method, they are simulated using computer programs with two different generations of stochastic processes, namely the Wiener process or Brownian motion, which are 1000 and 10000, respectively.