Numerical solution of second-order stochastic dierential equations with Gaussian random parameters

In this paper, we present the numerical solution of ordinary dierential equations (or SDEs), from each order especially second-order with time-varying and Gaussian random coecients. We indicate a complete analysis for second-order equations in special case of scalar linear second-order eq uations (damped harmonic oscillators with additive or multi- plicative noises). Making stochastic dierential equations system from this equation, it could be approximated or solved numerically by dierent numerical methods. In the case of linear stochastic dierential equations system by Computing fundamental matrix of this system, it could be calculated based on the exact solution of this system. Finally, this stochastic equa- tion is solved by numerically method like Euler-Maruyama and Milstein. Also its Asymptotic stability and statistical concepts like expectation and variance of solutions are discussed.