The 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.

Исследованы вопросы устойчивости и сходимости предложенного ранее матричного метода численного интегрирования краевых задач с граничными условиями первого, второго и третьего рода для линейных неоднородных обыкновенных дифференциальных уравнений второго порядка с переменными коэффициентами. Использование многочленов Тейлора произвольных степеней позволило повысить порядок аппроксимации метода до произвольного натурального числа и отказаться от аппроксимации производных конечными разностями. При выборе второй степени многочлена Тейлора уравнения метода совпали с известными уравнениями традиционного метода численного интегрирования краевых задач, в котором производные аппроксимированы конечными разностями. Показано, что достаточный критерий устойчивости при использовании в методе многочленов Тейлора степени три и выше совпадает с достаточным критерием устойчивости традиционного метода сеток для численного интегрирования краевых задач с граничными условиями первого, второго и третьего рода. Теоретически установлено, что скорость сходимости матричного метода интегрирования краевых задач с граничными условиями первого рода пропорциональна степени используемого многочлена Тейлора в случае, когда эта степень является чётной, и пропорциональна числу, меньшему степени на единицу, когда эта степень является нечётной; при интегрировании краевых задач с граничными условиями второго и третьего рода скорость сходимости метода пропорциональна степени используемого многочлена Тейлора независимо от её чётности и меньше неё на единицу. Полученные теоретические результаты подтверждены численным экспериментом.