Numerical integration by the matrix method of boundary value problems for linear inhomogeneous ordinary differential equations of the third order with variable coefficients

V. N. Maklakov ; Ya. G. Stelmakh

Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences, Tome 22 (2018) no. 1, pp. 153-183.

Résumé

The use of the Taylor polynomial of the second degree when approximating the derivatives by finite differences leads to the second order of approximation of the traditional method of nets in the numerical integration of second-order ordinary differential equations with variable coefficients. In the study of boundary-value problems for the third-order ordinary differential equations with variable coefficients, we offer the previously proposed method of numerical integration, using the means of the matrix calculus, in which approximation of the derivatives by finite differences was not used. According to this method, in the construction of a system of difference equations, an arbitrary power of the Taylor polynomial in the expansion of the desired solution of the problem in a Taylor series can be chosen. The disparity is calculated and an estimate of the order of approximation of the method is given depending on the chosen degree of the Taylor polynomial using the four-point pattern. The regularities between the order of approximation of the matrix method and the degree of the used Taylor polynomial are theoretically revealed. We found out that the order of approximation is proportional to the degree of the used Taylor polynomial and less by two than it. We propose a procedure for constructing a fictitious boundary condition that allows us to construct a closed system of difference equations for the matrix method of numerical integration. The system of difference equations is divided into two subsystems: the first subsystem consists of two equations, the first of which contains the given value of the derivative in the boundary conditions of the problem, the second one contains the value calculated from the fictitious boundary condition; the second subsystem consists of the remaining difference equations of the constructed closed system. The disparity is calculated and an estimate of the order of approximation of the method is given depending on the chosen degree of the Taylor polynomial using the five-point pattern. The regularities between the order of approximation of the matrix method and the degree of the used Taylor polynomial are theoretically revealed. The following is revealed: a) the order of approximation of the first subsystem, the second subsystem with an even value of the degree of the Taylor polynomial and the whole problem is proportional to this degree and less than it by two; b) the order of approximation of the second subsystem with an odd value of the degree of the Taylor polynomial is proportional to this degree and less than it by one.

Détail

Keywords: ordinary differential equations, boundary value problems, approximation order, numerical methods, Taylor polynomials.

Bibliographie

[1] Radchenko V. P., Usov A. A., “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 [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2008, no. 2(17), 60–65 (In Russian) | DOI

[2] 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

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

[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, Part I: Survey Lectures, ed. A. K. Aziz, Academic Press, New York, 1975, 27–88 | DOI | MR

[5] Godunov S. K., Ryabenki V. S., Theory of Difference Schemes: An Introduction, Wiley, New York, 1964, xii+289 pp. | MR

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

[7] Samarskii A. A., Teoriia raznostnykh skhem [The Theory of Difference Schemes], Nauka, Moscow, 1977, 656 pp. (In Russian) | MR

[8] Samarskii A. A., Gulin A. V., Chislennye metody [Numerical methods], Nauka, Moscow, 1973, 432 pp. (In Russian) | MR

[9] Samarskii A. A., Gulin A. V., Ustoichivost' raznostnykh skhem [The stability of difference schemes], Nauka, Moscow, 1973, 416 pp. (In Russian) | MR

[10] Boutayeb A., Chetouani A., “Global Extrapolations Of Numerical Methods For Solving A Parabolic Problem With Non Local Boundary Conditions”, International Journal of Computer Mathematics, 80:6 (2003), 789–797 | DOI | MR | Zbl

[11] 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 http://www.m-hikari.com/ams/ams-password-2007/ams-password41-44-2007/boutayebAMS41-44-2007.pdf | MR

[12] Maklakov V. N., “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”, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], no. 3(36), 143–160 (In Russian) | DOI

[13] Maklakov V. N., “Convergence of the matrix method of numerical integration of the boundary value problems for linear nonhomogeneous ordinary differential second order equations with variable coefficients”, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 19:3 (2015), 559–577 (In Russian) | DOI | Zbl

[14] Fichtenholz G. M., Differential- und Integralrechnung. I [Differential and integral calculus. I], Hochschulbücher für Mathematik [University Books for Mathematics], 61, VEB Deutscher Verlag der Wissenschaften, Berlin, 1986, xiv+572 pp. (In German) | MR | Zbl

[15] Kurosh A., Higher algebra, Mir Publ., Moscow, 1972, 428 pp. | MR | MR | Zbl

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

[17] Kamke E., Spravochnik po obyknovennym differentsial'nym uravneniiam [Manual of ordinary differential equations], Nauka, Moscow, 1976, 576 pp. (In Russian) | MR

[18] Maklakov V. N., “The evaluation of the order of approximation of the matrix method for numerical integration of the boundary value problems for systems of linear non-homogeneous ordinary differential equations of the second order with variable coefficients. Message 1. Boundary value problems with boundary conditions of the first kind”, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 20:3 (2016), 389–409 (In Russian) | DOI | Zbl