This paper introduces analytical and numerical solutions of the nonlinear Langevin’s equation under square nonlinearity with stochastic non-homogeneity. The solution is obtained by using the Wiener-Hermite expansion with perturbation (WHEP) technique, and the results are compared with those of Picard iterations and the homotopy perturbation method (HPM). The WHEP technique is used to obtain up to fourth order approximation for different number of corrections. The mean and variance of the solution are obtained and compared among the different methods, and some parametric studies are done by using Matlab.