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

NLEVSD - решение нелинейной задачи на собственные значения для системы дифференциальных уравнений первого порядка


Автор: Т.Бояджиев
rus
Вы
counter
посетитель.

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


Программа NLEVSD служит для решения нелинейной задачи на собственные значения для системы из NDE <= 32 обыкновенных дифференциальных уравнений первого порядка вида:

Fm ( x,y,y',p ) = 0,   m = 1,2, ... NDE
на интервале (a, b), где a < b, с граничными условиями:
L1 ( a, y(a) ) = 0,   l = 1,2,...,NBCL,
Rr( b, y(b) ) = 0,    r = NBCL+1,...,NDE.
Здесь p - искомый параметр.
Предполагается, что задано также "условие нормировки" в точке XP ( a <= XP <= b ):
H (y,p) = 0.

Для поиска решения { y(x), p } применяется непрерывный аналог метода Ньютона вместе с методом сплайн - коллокации в гауссовых узлах сетки.

Обращение:

CALL NLEVSD (X,Y,PRM,XP,NDE,NBCL,NPAR,NNOD,HX,W,U,V,YO,RHSU,RHSV,INTEGS,
TIMES,NITEM,ITER,EPS,FUN,DFUN,BCOND,DBCOND,PCOND,DPCOND,DETAILS,IFAIL)

X - сетка из NNOD узлов (задается пользователем);
Y - массив из 2NDE*NNOD элементов; на входе содержит начальное приближение, на выходе содержит решение y(x) в узлах сетки;
PRM - на выходе содержит собственное значение p;
NPAR - = 1
HX - массив из NNOD-1 шагов сетки (задается пользователем);
W - рабочий массив из 8*NNOD*NDE**2 элементов;
U,V,YO,RHSU,RHSV - рабочие массивы из 2*NDE*NNOD элементов;
INTEGS - INTEGER массив размера 2*NNOD-2;
TIMES - число обращений к подпрограмме;
NITEM - максимальное число итераций (задается пользователем);
ITER - число итераций;
EPS - абсолютная точность (задается пользователем);
BCOND - подпрограмма вычисления левых частей граничных условий;
DBCOND - подпрограмма вычисления производных Фреше левых частей граничных условий;
PCOND - подпрограмма вычисления левой части условия;
DPCOND - подпрограмма вычисления производной Фреше левой части условия;
DETAILS - (integer) уровень печати
IFAIL - флаг итераций (IFAIL=0 - итерационный процесс сходится за ITER итераций).

Архив программы с исходными текстами.




home up e-mail