  Two specialized algorithms for the numerical integration of the equations of
motion of a Brownian walker obeying detailed balance are introduced. The
algorithms become symplectic in the appropriate limits, and reproduce the
equilibrium distributions to some higher order in the integration time step.
Comparisons with other existing integration schemes are carried out both for
static and dynamical quantities.
