https://bjlkeng.github.io/posts/hamiltonian-monte-carlo/ "...we cannot exactly simulate this system because, in general, Hamilton's equations do not yield a closed form solution."
Ну, если вы вооружены
статьёй под названием Generating Power of Lazy Semantics, то вас это не должно останавливать. Из Гамильтониана вполне можно получить систему уравнений вида a=f(x,v), где a это вектор ускорения, v - вектор скорости и x это вектор координат, при этом x это интеграл от скорости по времени с константой интегрирования x(0).
Что позволяет получать сколь угодно близкие приближения к точному решению.