Numerical integration by the matrix method and evaluation of the approximation order of difference boundary value problems for non-homogeneous linear ordinary differential equations of the fourth order with variable coefficients

V. N. Maklakov ; M. A. Ilicheva

Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences, Tome 24 (2020) no. 1, pp. 137-162.

Résumé

The use of the second degree Taylor polynomial in approximation of derivatives by finite difference method leads to the second order approximation of the traditional grid method for numerical integration of boundary value problems for non-homogeneous linear ordinary differential equations of the second order with variable coefficients. The study considers a previously proposed method of numerical integration using matrix calculus which didn’t include the approximation of derivatives by finite difference method for boundary value problems of non-homogeneous fourth-order linear ordinary differential equations with variable coefficients. According to this method, when creating a system of difference equations, an arbitrary degree of the Taylor polynomial can be chosen in the expansion of the sought-for solution of the problem into a Taylor series. In this paper, the possible boundary conditions of a differential boundary value problem are written both in the form of derived degrees from zero to three, and in the form of linear combinations of these degrees. The boundary problem is called symmetric if the numbers of the boundary conditions in the left and right boundaries coincide and are equal to two, otherwise it is asymmetric. For a differential boundary value problem, an approximate difference boundary value problem in the form of two subsystems has been built. The first subsystem includes equations for which the boundary conditions of the boundary value problem were not used; the second one includes four equations in the construction of which the boundary conditions of the problem were used. Theoretically, the patterns between the order of approximation and the degree of the Taylor polynomial were identified. The results are as follows: a) the approximation order of the first and second subsystems is proportional to the degree of the Taylor polynomial used; b) the approximation order of the first subsystem is two units less than the degree Taylor polynomial with its even value and three units less with its odd value; c) the approximation order of the second subsystem is three units less than the degree Taylor polynomial regardless of both even-parity or odd-parity of this degree, and the degree of the highest derivative in the boundary conditions of the boundary value problem. The approximation order of the difference boundary value problem with all possible combinations of boundary conditions is calculated. The theoretical conclusions are confirmed by numerical experiments.

Détail

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

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 (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”, Intern. J. Comp. Math., 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”, Appl. Math. Sci., 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., “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 2. Boundary value problems with boundary conditions of the second and third kind”, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 21:1 (2017), 55–79 (In Russian) | DOI

[14] Maklakov V. N., Stelmakh Ya. G., “Numerical integration by the matrix method of boundary value problems for linear inhomogeneous ordinary differential equations of the third order with variable coefficients”, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 22:1 (2018), 153–183 (In Russian) | DOI | Zbl

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

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

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

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