We propose a numerical scheme for the integration of the Langevin equation which is second-order accurate. More importantly, we indicate how to generalize this scheme to situations where holonomic constraints are added and show that the resulting scheme remains second-order accurate. (c) 2006 Elsevier B.V. All rights reserved.