erical integration for Hamiltonian systems,
including a derivation of the error bound for the leapfrog
integrator.</p>
</div>
<div id="metropolis-accept-step" class="section level3">
<h3>Metropolis Accept Step</h3>
<p>If the leapfrog integrator were perfect numerically, there would no
need to do any more randomization per transition than generating a
random momentum vector. Instead, what is done in practice to account
for numerical errors during integration is to apply a Metropolis
acceptance step, where th