Journal of Samara State Technical University, Ser. Physical and Mathematical SciencesJournal of Samara State Technical University, Ser. Physical and Mathematical Sciences1991-86152310-7081Samara State Technical University2047210.14498/vsgtu1426Research ArticleConvergence of the matrix method of numerical integration of the boundary value problems for linear nonhomogeneous ordinary differential second order equations with variable coefficientsMaklakovVladimir N(Cand. Phys. & Math. Sci.), Associate Professor, Dept. of Higher Mathematics and Applied Informaticsmakvo63@yandex.ruSamara State Technical University1509201519355957714022020Copyright © 2015, Samara State Technical University2015The problems of stability and convergence of previously proposed matrix method of numerical integration of boundary value problems with boundary conditions of the first, second and third kinds of nonhomogeneous linear ordinary differential second order equations with variable coeffcients are considered. Using of the Taylor polynomials of arbitrary degrees allowed to increase the approximation order of the method to an arbitrary natural number and to refuse from the approximation of derivatives by finite differences. When choosing the second degree Taylor polynomials the equation of the method coincided with the known equations of the traditional method of numerical integration of the boundary value problems where the derivatives are approximated by finite differences. It was shown that a suffcient criterion of stability when used in the method of Taylor polynomials of the third degree and more coincides with the suffcient criterion of stability of the traditional grid method for the numerical integration of boundary value problems with boundary conditions of the first, second and third kind. Theoretically, it is established that the degree of convergence of the matrix method for integration of boundary value problems with boundary conditions of the first kind is proportional to the degree of the used Taylor polynomials in the case, when the degree is even, and is proportional to the number that is one less than the degree if it is odd; when integrating the boundary value problems with boundary conditions of the second and third kind the degree of convergence of the method is proportional to the degree of the used Taylor polynomials regardless of its parity and one less than it. The obtained theoretical results are confirmed by numerical experiments.ordinary differential equationsboundary value problemsapproximation orderstabilityconvergencenumerical methodsTaylor polynomialsобыкновенные дифференциальные уравнениякраевые задачипорядок аппроксимацииустойчивостьсходимостьчисленные методымногочлены Тейлора[Keller H. B. Accurate difference methods for nonlinear two-point boundary value problems // SIAM J. Numer. Anal., 1974. vol. 11, no. 2. pp. 305-320. doi: 10.1137/0711028.][Lentini M., Pereyra V. A variable order finite difference method for nonlinear multipoint boundary value problems // Math. Comp., 1974. vol. 28, no. 128. pp. 981-1003. doi: 10.1090/s0025-5718-1974-0386281-4.][Keller H. B. Numerical solution of boundary value problems for ordinary differential equations: Survey and some resent results on difference methods / Numerical solutions of boundary value problems for ordinary differential equations; ed. A. K. Aziz. New York: Academic Press, 1975. pp. 27-88. doi: 10.1016/b978-0-12-068660-5.50007-7.][Годунов С. К., Рябенький В. С. Разностные схемы. М.: Наука, 1977. 439 с.][Формалеев В. Ф., Ревизников Д. Л. Численные методы. М.: Физматлит, 2004. 400 с.][Самарский А. А. Теория разностных схем. М.: Наука, 1977. 656 с.][Самарский А. А., Гулин А. В. Численные методы. М.: Наука, 1989. 432 с.][Самарский А. А., Гулин А. В. Устойчивость разностных схем. М.: Наука, 1973. 416 с.][Boutayeb A., Chetouani A. Global extrapolations of numerical methods for a parabolic problem with nonlocal boundary conditions // International Journal of Computer Mathematics, 2003. vol. 80, no. 6. pp. 789-797. doi: 10.1080/0020716021000039209.][Boutayeb A., Chetouani A. A numerical comparison of different methods applied to the solution of problems with non local boundary conditions // Applied Mathematical Sciences, 2007. vol. 1, no. 44. pp. 2173-2185, http://www.m-hikari.com/ams/ams-password-2007/ams-password41-44-2007/boutayebAMS41-44-2007.pdf.][Радченко В. П., Усов А. А. Модификация сеточных методов решения линейных дифференциальных уравнений с переменными коэффициентами на основе тейлоровских разложений // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2008. № 2(17). С. 60-65. doi: 10.14498/vsgtu646.][Маклаков В. Н. Оценка порядка аппроксимации матричного метода численного интегрирования краевых задач для линейных неоднородных обыкновенных дифференциальных уравнений второго порядка // Вестн. Сам. гос. техн. ун-та. Сер. Физ.мат. науки, 2014. № 3(36). С. 143-160. doi: 10.14498/vsgtu1364.][Закс Л. Статистическое оценивание. М.: Статистика, 1976. 598 с.][Турчак Л. И. Основы численных методов. М.: Наука, 1987. 320 с.]