|
PROGRAM LIBRARY JINRLIBPROGON4 - solution of the boundary problem for the ordinary differential equationAuthor: E.V.Zemlyanaya |
|
Language: Fortran Subroutine PROGON4 solves the ODE in the following form
Z"(x)+F(x)*Z'(x)+G(x)*Z(x)=K(x)
at the interval [a,b] with the boundary conditions:
D1*Z'(x=a)+F1*Z(x=a) = G1
D2*Z'(x=b)+F2*Z(x=b) = G2 The generalized Numerov's method is used for the 4-th order finite difference approximation. This approach was suggested and implemented in SLIPH4 [1] in the framework of Newtonian scheme for numerical solution of the Sturm-Liouville problem. Later, this method was generalized for the system of two ODEs in the PROGS2H4 [2]. References:
Usage: CALL PROGON4(N,H,F,G,UK,D1,D2,F1,F2,G1,G2,Z,AV,BV), where:
Sources and description (in Russian) are submitted. |