|
CONTIN-NLIN Библиотека "JINRLIB"
Автор: Е.В.Земляная Вы
Язык: Фортран
посетитель.
ПРОГРАММНЫЙ КОМПЛЕКС ЧИСЛЕННОГО ПРОДОЛЖЕНИЯ ПО ПАРАМЕТРУ
Комплекс программ CONTIN-NLIN реализует численное продолжение по параметру
применительно к следующей системе двух нелинейных уравнений:
u1"(x)+F11(u1,u2,x, )u1'(x)+F12(u1,u2,x, )u2'(x)+G1(u1,u2,x, ) = 0 ,
u2"(x)+F21(u1,u2,x, )u1'(x)+F22(u1,u2,x, )u2'(x)+G2(u1,u2,x, ) = 0 ,
где a x b, - вектор параметров, Fij и Gi (i,j=1,2) - функции,
обеспечивающие существование решения (они задаются пользователем).
Продолжение ведется по одному из элементов вектора .
Краевые условия на функции 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.
|