  A natural method for the introduction of second-order derivatives of the log
likelihood into MCMC algorithms is introduced, based on Taylor expansion of the
Langevin equation followed by exact solution of the truncated system.
