|
CONTIN-NLIN Library "JINRLIB"
Author: E.V.Zemlyanaya You are
Language: Fortran
visitor here.
PROGRAM COMPLEX FOR NUMERICAL CONTINUATION
CONTIN-NLIN is the program for numerical continuation of solutions of the
following system of two nonlinear ODEs:
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 ,
where a x b, - vector of parameters, Fij and Gi (i,j=1,2) - given
functions, solution is continued in one of elements of vector ,
boundary conditions have the following form:
Di*ui'(a) + Qi*ui(a) = 0, DDi*ui'(b) + QQi*ui(b) = 0.
Here i=1,2, Di2 + Qi2 > 0, DDi2 + QQi2 > 0.
The continuation scheme through turning points is described in [1].
At each step of numerical continuation, the Newtonian iteration scheme [2]
with the 4th accuracy finite-difference approximation [3] is used.
References:
-----------
1. E.V.Zemlyanaya, I.V.Barashenkov. Mathematical Modeling. V.16, No.10,
pp.3-14, 2004.
2. I.V.Puzynin, I.V.Amirkhanov, E.V.Zemlyanaya et al. Particles & Nuclei.
v.30, No.1, pp.210-265, 1999.
3. E.V.Zemlyanaya, I.V.Puzynin, T.P.Puzynina. JINR Comm. P11-97-414, Dubna,
1997.
Sources, tests and the detailed description (in Russian, pdf) are submitted.
|