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 axb, - 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. |