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 , где axb, - вектор параметров, 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. |