Second-order integrators for Langevin equations with holonomic constraints

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.