This study unveils a novel approach, integrating the variational iteration method and numerical integration, to address the n-th order fuzzy random ordinary differential equations' linear fuzzy initial value problems. The robustness of the variational iteration method, a proven and reliable technique, ensures the effectiveness of the proposed approach. The sequence of approximations generated by this method is scrutinized to confirm its convergence towards the exact solution, demonstrating the method's precision. Two distinct examples, each with a different number of Brownian motion generations, are simulated to elucidate the practical application of the proposed approach. The outcomes affirm the method's reliability and efficiency in tackling such complex mathematical problems.