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

---

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

u₁"(x)+F₁₁(u₁,u₂,x,)u₁'(x)+F₁₂(u₁,u₂,x,)u₂'(x)+G₁(u₁,u₂,x,) = 0 ,
u₂"(x)+F₂₁(u₁,u₂,x,)u₁'(x)+F₂₂(u₁,u₂,x,)u₂'(x)+G₂(u₁,u₂,x,) = 0 ,

где axb, - вектор параметров, F_ij и G_i (i,j=1,2) - функции,
обеспечивающие существование решения (они задаются пользователем).
Продолжение ведется по одному из элементов вектора .
Краевые условия на функции u₁ и u₂ имеют вид:

D_i*u_i'(a) + Q_i*u_i(a) = 0,      DD_i*u_i'(b) + QQ_i*u_i(b) = 0.
Здесь i=1,2, D_i² + Q_i² > 0,    DD_i² + QQ_i² > 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.

---