This has motivated the consideration of random integral numerical methods to approximate infinite mean square convergent integrals. In fact, random Gauss-Hermite quadrature formulae are proposed to approximate the solution stochastic process in a more computable way.