ОБЪЕДИНЕННЫЙ   ИНСТИТУТ   ЯДЕРНЫХ   ИССЛЕДОВАНИЙ
lit БИБЛИОТЕКА   ПРОГРАММ   JINRLIB

CONTIN-NLIN - программный комплекс численного продолжения по параметру


Автор: Е.В.Земляная
eng
Вы
counter
посетитель.

Язык: Фортран


Комплекс программ CONTIN-NLIN реализует численное продолжение по параметру применительно к следующей системе двух нелинейных уравнений:

u1"(x)+F11(u1,u2,x,a)u1'(x)+F12(u1,u2,x,a)u2'(x)+G1(u1,u2,x,a) = 0 ,
u2"(x)+F21(u1,u2,x,a)u1'(x)+F22(u1,u2,x,a)u2'(x)+G2(u1,u2,x,a) = 0 ,

где alexleb, a - вектор параметров, Fij и Gi (i,j=1,2) - функции,
обеспечивающие существование решения (они задаются пользователем).
Продолжение ведется по одному из элементов вектора a.
Краевые условия на функции u1 и u2 имеют вид:

Di*ui'(a) + Qi*ui(a) = 0,      DDi*ui'(b) + QQi*ui(b) = 0.
Здесь i=1,2, Di2 + Qi2 > 0,    DDi2 + QQi2 > 0.

В программе реализована описанная в [1] схема численного продолжения с возможностью продолжения по параметру через точки поворота. На каждом шаге численного продолжения используется ньютоновская итерационная схема [2] c конечно-разностной аппроксимацией четвертого порядка [3].

Архив программы включает в себя программу contin-nlin.f и тесты
(главная программа + п/п пользователя, входные и выходные файлы).

Имеется подробное описание (в формате .pdf).

Литература:

  1. Е.В.Земляная, И.В.Барашенков. Математическое моделирование, т.16,
    вып.10, сс.3-14, 2004.
  2. И.В.Пузынин, И.В.Амирханов, Е.В.Земляная и др. ЭЧАЯ, Т.30, вып.1,
    сс.210-265, 1999.
  3. Е.В.Земляная, И.В.Пузынин, Т.П.Пузынина. Сообщение ОИЯИ P11-97-414,
    Дубна, 1997.



home up e-mail