Помощь с симплектическими интеграторами

Я пытаюсь разработать физическое моделирование, и я хочу реализовать метод симплектического интегрирования четвертого порядка . Проблема в том, что я, должно быть, ошибаюсь в математике, поскольку моя симуляция вообще не работает при использовании симплектического интегратора (по сравнению с интегратором Рунге-Кутты четвертого порядка, который достаточно хорошо работает для симуляции). Я искал это в Google вечно, и все, что я могу найти, - это научные статьи на эту тему. Я попытался адаптировать метод, использованный в статьях, но безуспешно. Я хочу знать, есть ли у кого-нибудь исходный код для моделирования, в котором используются симплектические интеграторы, предпочтительно для моделирования гравитационного поля, но подойдет любой симплектический интегратор. На каком языке находится исходный текст, не имеет большого значения, но я был бы признателен за язык, использующий синтаксис в стиле C. Спасибо!