Language: Fortran

In this work the description of the computational scheme and procedure
realizing the sweep method for solving the system of two ordinary
differential equations with 2th type boundary conditions at the ends points
of given interval is represented. The boundary problem is approximated
with the help of the difference scheme of O(h^4) precision on the mesh
with constant step h. This problem occurs in schemes of linearity of
nonlinear boundary problems (for example, in newtonian iterative schemes),
and represented here approach gives the possibility to receive the precision
of approximation with power O(h^4).

Sources and the detailed description (in Russian) are submitted.