The second order Langevin equation and numerical simulations