This paper deals with the problem of correct asymptotic dissipation
of mean energy functional related to numerical integration of systems of uncoupled
linear oscillators under random perturbations. It is shown that the drift-implicit
trapezoidal method provides numerical approximations which possess the correct asymptotic behavior of their mean energy functional compared to that of the underlying
exact solution as integration time t advances to infinity.