Здравствуйте. Проблема такая. Я написала программу для решения задачи Римана методом SPH. Но она не работает. И сложность в том что определить в чем ошибка (в формуле или в коде) я не могу. Прошу помощи. Если у кого то есть рабочая программа, поделитесь пжл, чтобы имела возможность сравнить промежуточные результаты и найти ошибку.
Система уравнений:
Кратко описываю программу.
Имеется два класса. первый "частица", с параметрами - координата(одномерный случай), плотность, давление, энергия и скорость. Второй класс "метод sph", в нем все вычисления происходят.
Начальные значения параметров задаются в конструкторе. Остальные считываются с файла const.
Интегрирование проводится по координате, энергии и скорости. Методом эйлера.
Результат заносится в файл.
Самое интересное в том, что на первый шагах результат есть(не знаю правильный или нет), а потом значения неопределенные становятся.
Очень надеюсь на помощь. Я в полном отчаянии.
Вот ссылка на программу
http://rghost.ru/38307097