|
|
|
MathFortran |
|
|||
Новичок Профиль Группа: Участник Сообщений: 4 Регистрация: 22.4.2014 Репутация: нет Всего: нет |
Добрый день!
Есть система алгебраических уравнений, которую надо решать численно (методом Ньютона, например). Я нашел программу здесь: http://www.zib.de/en/numerik/software/ant/nleq1.html Но никак не могу разобраться как внутри ввести нужные мне уравнения. Буду рад любой помощи. В частности, если кто-нибудь подскажет, где найти более прозрачную (user-friendly) программку. |
|||
|
||||
Фантом |
|
|||
Вы это прекратите! Профиль Группа: Участник Клуба Сообщений: 1516 Регистрация: 23.3.2008 Репутация: 5 Всего: 49 |
Вообще говоря, метод Ньютона проще написать самому, это исходник на пару килобайт.
Сколько уравнений решать надо? |
|||
|
||||
MathFortran |
|
|||
Новичок Профиль Группа: Участник Сообщений: 4 Регистрация: 22.4.2014 Репутация: нет Всего: нет |
Там система уравнений химической кинетики для нахождения равновесия.
Такие системы получаются жесткие, поэтому в лоб обычно метода не проходит. В принципе, судя по статьям солвер в ссылке достаточно хороший просто бестолково документация написана |
|||
|
||||
Фантом |
|
|||
Вы это прекратите! Профиль Группа: Участник Клуба Сообщений: 1516 Регистрация: 23.3.2008 Репутация: 5 Всего: 49 |
Уравнения химической кинетики - дифференциальные, а не алгебраические. Жесткие системы, естественно, тоже. Поэтому применение к ним метода Ньютона - это что-то, мягко говоря, странное. |
|||
|
||||
MathFortran |
|
|||
Новичок Профиль Группа: Участник Сообщений: 4 Регистрация: 22.4.2014 Репутация: нет Всего: нет |
Чтобы найти равновесные концентрации, производные от концентраций по времени приравнивают к нулю.
Отсюда получается система алгебраических уравнений, которая в целом тоже достаточно неприятная в вычислительном плане. |
|||
|
||||
Фантом |
|
|||
Вы это прекратите! Профиль Группа: Участник Клуба Сообщений: 1516 Регистрация: 23.3.2008 Репутация: 5 Всего: 49 |
||||
|
||||
MathFortran |
|
|||
Новичок Профиль Группа: Участник Сообщений: 4 Регистрация: 22.4.2014 Репутация: нет Всего: нет |
немного: 5-6 уравнений
|
|||
|
||||
Фантом |
|
|||
Вы это прекратите! Профиль Группа: Участник Клуба Сообщений: 1516 Регистрация: 23.3.2008 Репутация: 5 Всего: 49 |
Тогда это стрельба из пушки по воробьям, проще написать свою. Если есть возможность, дайте ссылку на уравнения и величины параметров.
|
|||
|
||||
tzirechnoy |
|
|||
Эксперт Профиль Группа: Завсегдатай Сообщений: 1173 Регистрация: 30.1.2009 Репутация: нет Всего: 16 |
Честно говоря, по-моему проще как раз наоборот не писать свою, а взять какой-нибудь octave. Или scilab. В котором это делается в число уравнений плюс 3 строчек, и описано в примерах для начинающих.
|
|||
|
||||
Фантом |
|
|||
Вы это прекратите! Профиль Группа: Участник Клуба Сообщений: 1516 Регистрация: 23.3.2008 Репутация: 5 Всего: 49 |
||||
|
||||
kovalkos |
|
|||
Новичок Профиль Группа: Участник Сообщений: 2 Регистрация: 29.4.2014 Где: Луганск Репутация: нет Всего: нет |
Решай методом Гаусса:)) Если метод Ньютона не подходит:))
|
|||
|
||||
Santik |
|
|||
Новичок Профиль Группа: Участник Сообщений: 45 Регистрация: 13.3.2012 Где: Мирный (Якутия) Репутация: нет Всего: нет |
А библиотеку IMSL почему не используем???
|
|||
|
||||
0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей) | |
0 Пользователей: | |
« Предыдущая тема | Fortran | Следующая тема » |
|
По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности Powered by Invision Power Board(R) 1.3 © 2003 IPS, Inc. |