In this paper, a numerical approach is proposed based on the variation-of-constants formula for the numerical discretization Langevin-type equations. Linear and non-linear cases are treated separately. The proofs of ...