An Efficient Algorithm for Solving General Linear Two-Point BVP

A new method is described to compute the solutions of linear BVP in an efficient and stable way. The stability is achieved by decoupling the multiple shooting recursion; this means that the choice of output points can be made virtually without regard to restrictions. By fixing the number of integration steps per “shooting” interval and assembling as many of them as is needed to fit the user's requirements, high efficiency is gained. Apart from a mathematical description, we also give a stability analysis of the method. A large number of numerical examples confirm this analysis and illustrate the possibilities of the algorithm.

[1]  S. D. Conte.  The Numerical Solution of Linear Boundary Value Problems , 1966 .

[2]  S. M. Roberts,et al.  Two-point boundary value problems: shooting methods , 1972 .

[3]  P. Deuflhard.  A modified Newton method for the solution of ill-conditioned systems of nonlinear equations with application to multiple shooting , 1974 .

[4]  V. Pereyra,et al.  An adaptive finite difference solver for nonlinear two point boundary problems with mild boundary layers. , 1975 .

[5]  H. Keller.  Numerical Solution of Two Point Boundary Value Problems , 1976 .

[6]  J. Varah.  Alternate Row and Column Elimination for Solving Certain Linear Systems , 1976 .

[7]  H. A. Watts,et al.  Computational Solution of Linear Two-Point Boundary Value Problems via Orthonormalization , 1977 .

[8]  H. Schwetlick,et al.  Stoer, J. / Bulirsch, R., Einführung in die Numerische Mathematik II, IX, 286 S., 1973. DM 14,80, US $ 5.50. Berlin-Heidelberg-New York. Springer-Verlag , 1978 .

[9]  Michael R. Osborne.  THE STABILIZED MARCH IS STABLE , 1979 .

[10]  R.M.M. Mattheij,et al.  Characterizations of dominant and dominated solutions of linear recursions , 1980 .

[11]  Richard Weiss,et al.  SOLVEBLOK: A Package for Solving Almost Block Diagonal Linear Systems , 1980, TOMS.

[12]  R.M.M. Mattheij.  Stable computation of solutions of unstable linear initial value recursions , 1982 .

[13]  R. Mattheij.  The conditioning of linear boundary value problems , 1982 .

[14]  R.M.M. Mattheij.  The stability of LU-decompositions of block tridiagonal matrices , 1984, Bulletin of the Australian Mathematical Society.

[15]  R. Mattheij,et al.  On optimal shooting intervals , 1984 .