Ищите ошибку.
Бяка пропадает при уменьшении шага по времени. Видимо просто округление.
В связи с чем вопрос: как вообще в этом методе "по научному" шаг выбирают?
Кстати, вот это
Ну вот, к примеру симплектический метод Эйлера:
смотрится как-то странно. Я, конечно, не коопенгаген, но вроде бы основная идея этих интеграторов (насколько я их понял) в замене экспоненциала суммы двух некоммутирующих операторов чересполосицей экспоненциалов, действующих сдвигами. И в таковом случае для симплектического Эйлера должно быть
то есть, каждый сдвиг проводится последовательно с использованием последних обновленных значений переменных.
И еще вроде бы расщепление
принципиально. Или нэ?
(Оффтоп)
Тьфу ты, дьявольщина! Не заметил, что тему в пургу снесли...