Estimation of the Order of the Matrix Method Approximation of Numerical Integration of Boundary-Value Problems for Inhomogeneous Linear Ordinary Differential Equations of the Second Order

V. N. Maklakov

Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences, no. 3 (2014), pp. 143-160.

Résumé

Using the first three terms of Taylor expansion of the required function in the approximate derivative by finite differences leads to the second order approximation of the traditional numerical quadrature method of boundary value problems for linear ordinary second order differential equations with variable coefficients. The paper shows previously proposed numerical quadrature method using tools of matrix calculus where the approximate derivative by finite differences was not used. Agreeing to above method the arbitrary number of terms of Taylor expansion for the required solution may be used when compiling the difference equation system. When using the three first terms of expansion the difference equation system coincided with the traditional system. The estimation of residuals and the order of approximation depending on the number of the used terms of Taylor expansion is given. It is theoretically shown that for the boundary value problem with boundary conditions of the first kind the approximation method order increases in direct proportion with the increasing in the number of members used in Taylor series expansion only for odd values of this number. For even values of this number the order of approximation coincides with the order of approximation for the number less by unit of the odd values. For boundary value problems with boundary conditions of the second and third kinds the order of approximation was directly proportional to the number of used terms in the Taylor series expansion of the required solution of the problem regardless of evenness. In these cases the order of approximation of the boundary points and therefore the whole problem turned out to be one unit less than the order for the inner points of the grid for the interval of integration. The method of approximation order increase at the boundary points up to the approximation order in the inner points of the grid is presented. The theoretical conclusions are confirmed by a numerical experiment for a boundary value problem with boundary conditions of the first and third kinds.

Détail

Keywords: ordinary differential equations of second order, boundary value problems, boundary conditions of the first, second and third types, order of approximation, numerical methods, Taylor polynomials.

Bibliographie

[1] Keller H. B., “Accurate Difference Methods for Nonlinear Two-point Boundary Value Problems”, SIAM J. Numer. Anal., 11:2 (1974), 305–320 | DOI | MR | Zbl

[2] Lentini M., Pereyra V., “A Variable Order Finite Difference Method for Nonlinear Multipoint Boundary Value Problems”, Mathematics of Computation, 28:128 (1974), 981–1003 | DOI | MR | Zbl

[3] 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, Academic Press, New York, 1975, 27–88 | DOI | MR

[4] Godunov S. K., Ryaben'kii V. S., Difference schemes. An introduction to the underlying theory, Studies in Mathematics and its Applications, 19, North-Holland, Amsterdam, New York, Oxford, Tokyo, 1987, xvii+489 pp. | DOI | MR | MR | Zbl | Zbl

[5] Samarskii A. A., The theory of difference schemes, Pure and Applied Mathematics, 240, Marcel Dekker, New York, NY, 2001, 786 pp. | DOI | MR | MR | Zbl | Zbl

[6] Formaleev V. F., Reviznikov D. L., Chislennye metody [Numerical Methods], Fizmatlit, Moscow, 2004, 400 pp. (In Russian)

[7] Boutayeb A., Chetouani A., “Global Extrapolations of Numerical Methods for a Parabolic Problem with Nonlocal Boundary Conditions”, International Journal of Computer Mathematics, 80:6 (2003), 789–797 | DOI | MR | Zbl

[8] Boutayeb A., Chetouani A., “A Numerical Comparison of Different Methods Applied to the Solution of Problems with Non Local Boundary Conditions”, Applied Mathematical Sciences, 1:44 (2007), 2173–2185 | MR

[9] V. P. Radchenko, A. A. Usov, “Modified grid method for solving linear differential equation equipped with variable coefficients based on Taylor series”, Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2008, no. 2(17), 60–65 (In Russian) | DOI | MR

[10] Maklakov V. N., “Numerical integration of the mixed boundary value problems for the second order inhomogeneous linear ordinary differential equations by a matrix method”, Sovremennyi nauchnyi vestnik, 2013, no. 16 (155), 72–78 (In Russian)

[11] Maklakov V. N., Usov A. A., “Numerical integration of the boundary value problems for the second order nonlinear ordinary differential equations by a matrix method with the use of iterative procedures”, Proceedings of the Ninth All-Russian Scientific Conference with international participation (21–23 May 2013). Part 3, Matem. Mod. Kraev. Zadachi, Samara State Technical Univ., Samara, 2013, 35–42 (In Russian)

[12] Turchak L. I., Osnovy chislennykh metodov [Foundations of numerical methods], Nauka, Moscow, 1987, 320 pp. (In Russian) | MR | Zbl

[13] Zaks L., Statisticheskoe otsenivanie [Statistical Estimation], Statistika, Moscow, 1976, 598 pp. (In Russian) | MR