Здравствуйте, помогите пожалуйста разобраться, что я делаю не так.
Используется разностная схема с центральной разностью для уравнения:

получается следующее:

где

-шаг по времени, а

- шаг по сетке.
Далее положим

и получим вот что

Ну и формируется трех диагональная матрица

Замечательно, это мы можем решить с помощью метода прогонки. Условия следующие: на правой половине расчетной области в начальный момент времени

а на левой

.

тогда:

Все верно? Правильно ли я мыслю? Конечно нет! Потому, что начинаю считать и получается полный бред. Где моя ошибка?