Replacing the finite difference methods for nonlinear two-point boundary value problems by successive application of the linear shooting method

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:358

Article

Filipov, Stefan M.¹  Gospodinov, Ivan D.¹  Farago, Istvan^2,3 

[1] Univ Chem Technol & Met, Dept Comp Sci, Sofia, Bulgaria

[2] Budapest Univ Technol & Econ, Inst Math, Dept Differential Equat, Budapest, Hungary

[3] MTA ELTE Res Grp, Budapest, Hungary

关键词: Finite difference method;    Newton;    Picard;    Constant-slope;    Quasi-linearization;    Linear shooting;   

DOI  :  10.1016/j.cam.2019.03.004

来源: Elsevier

【 摘 要 】

This paper studies numerical solution of nonlinear two-point boundary value problems for second order ordinary differential equations. First, it establishes a connection between the finite difference method and the quasi-linearization method. We prove that using finite differences to discretize the sequence of linear differential equations arising from quasi-linearization (Newton method on operator level) leads to the usual iteration formula of the Newton finite difference method. From the provided derivation, it can easily be inferred that such a relation holds also for the Picard and the constant-slope methods. Based on this result, we propose a way of replacing the Newton, Picard, and constant-slope finite difference methods by respective successive application of the linear shooting method. This approach has a number of advantages. It removes the necessity of solving systems of algebraic equations, hence working with matrices, altogether. Compared to the usual finite difference method with general solver, it reduces the number of computational operations from O(N-3), where N is the number of mesh-points, to only O(N). (C) 2019 Elsevier B.V. All rights reserved.

【 授权许可】

Free