Ищите ошибку.
Бяка пропадает при уменьшении шага по времени. Видимо просто округление.
В связи с чем вопрос: как вообще в этом методе "по научному" шаг выбирают?
Кстати, вот это
Ну вот, к примеру симплектический метод Эйлера:

смотрится как-то странно. Я, конечно, не коопенгаген, но вроде бы основная идея этих интеграторов (насколько я их понял) в замене экспоненциала суммы двух некоммутирующих операторов чересполосицей экспоненциалов, действующих сдвигами. И в таковом случае для симплектического Эйлера должно быть
![$$\[
\begin{gathered}
q' = q + \dot q\left( p \right)\tau \hfill \\
p' = p + \dot p\left( {q'} \right)\tau \hfill \\
\end{gathered}
\]
$$ $$\[
\begin{gathered}
q' = q + \dot q\left( p \right)\tau \hfill \\
p' = p + \dot p\left( {q'} \right)\tau \hfill \\
\end{gathered}
\]
$$](https://dxdy-02.korotkov.co.uk/f/9/b/9/9b937a75f0a82df7f808bc4f884055fa82.png)
то есть, каждый сдвиг проводится последовательно с использованием последних обновленных значений переменных.
И еще вроде бы расщепление
![$\[
H\left( {p,q} \right) = A\left( p \right) + B\left( q \right)
\]
$ $\[
H\left( {p,q} \right) = A\left( p \right) + B\left( q \right)
\]
$](https://dxdy-04.korotkov.co.uk/f/f/4/8/f48c6626d535d8e437a2a2c37255490682.png)
принципиально. Или нэ?
(Оффтоп)
Тьфу ты, дьявольщина! Не заметил, что тему в пургу снесли...