Такие системы получаются жесткие, поэтому в лоб обычно метода не проходит.
В принципе, судя по статьям солвер в ссылке достаточно хороший просто бестолково документация написана
| Версия для печати темы
Нажмите сюда для просмотра этой темы в оригинальном формате |
| Форум программистов > Fortran > Программа для решения систем уравнений численно |
| Автор: MathFortran 22.4.2014, 13:14 |
| Добрый день! Есть система алгебраических уравнений, которую надо решать численно (методом Ньютона, например). Я нашел программу здесь: http://www.zib.de/en/numerik/software/ant/nleq1.html Но никак не могу разобраться как внутри ввести нужные мне уравнения. Буду рад любой помощи. В частности, если кто-нибудь подскажет, где найти более прозрачную (user-friendly) программку. |
| Автор: Фантом 22.4.2014, 16:09 |
| Вообще говоря, метод Ньютона проще написать самому, это исходник на пару килобайт. Сколько уравнений решать надо? |
| Автор: MathFortran 22.4.2014, 16:15 |
| Там система уравнений химической кинетики для нахождения равновесия. Такие системы получаются жесткие, поэтому в лоб обычно метода не проходит. В принципе, судя по статьям солвер в ссылке достаточно хороший просто бестолково документация написана |
| Автор: Фантом 22.4.2014, 16:46 | ||
Уравнения химической кинетики - дифференциальные, а не алгебраические. Жесткие системы, естественно, тоже. Поэтому применение к ним метода Ньютона - это что-то, мягко говоря, странное. |
| Автор: MathFortran 22.4.2014, 17:18 |
| Чтобы найти равновесные концентрации, производные от концентраций по времени приравнивают к нулю. Отсюда получается система алгебраических уравнений, которая в целом тоже достаточно неприятная в вычислительном плане. |
| Автор: Фантом 22.4.2014, 22:26 | ||
Тогда это не совсем кинетика.  Ну ладно, с этим понятно. А все-таки - сколько уравнений в системе? |
| Автор: MathFortran 23.4.2014, 08:53 |
| немного: 5-6 уравнений |
| Автор: Фантом 23.4.2014, 15:53 |
| Тогда это стрельба из пушки по воробьям, проще написать свою. Если есть возможность, дайте ссылку на уравнения и величины параметров. |
| Автор: tzirechnoy 24.4.2014, 23:08 |
| Честно говоря, по-моему проще как раз наоборот не писать свою, а взять какой-нибудь octave. Или scilab. В котором это делается в число уравнений плюс 3 строчек, и описано в примерах для начинающих. |
| Автор: Фантом 25.4.2014, 00:50 | ||
Своя будет не намного длиннее. |
| Автор: kovalkos 30.4.2014, 07:53 |
| Решай методом Гаусса:)) Если метод Ньютона не подходит:)) |
| Автор: Santik 1.3.2016, 07:53 |
| А библиотеку IMSL почему не используем??? |
