The Open UniversitySkip to content

Multiple shooting using a dichotomically stable integrator for solving differential-algebraic equations

England, Roland; Lamour, René and López-Estrada, Jesús (2002). Multiple shooting using a dichotomically stable integrator for solving differential-algebraic equations. Applied Numerical Mathematics, 42(1-3) pp. 117–131.

DOI (Digital Object Identifier) Link:
Google Scholar: Look up in Google Scholar


In previous work by the first author, it has been established that a dichotomically stable discretization is needed when solving a stiff boundary-value problem in ordinary differential equations (ODEs), when sharp boundary layers may occur at each end of the interval. A dichotomically stable implicit Runge–Kutta method, using the 3-stage, fourth-order, Lobatto IIIA formulae, has been implemented in a variable step-size initial-value integrator, which could be used in a multiple-shooting approach.

In the case of index-one differential–algebraic equations (DAEs) the use of the Lobatto IIIA formulae has an advantage, over a comparable Gaussian method, that the order is the same for both differential and algebraic variables, and there is no need to treat them separately.

The ODE integrator (SYMIRK [R. England, R.M.M. Mattheij, in: Lecture Notes in Math., Vol. 1230, Springer, 1986, pp. 221–234]) has been adapted for the solution of index-one DAEs, and the resulting integrator (SYMDAE) has been inserted into the multiple-shooting code (MSHDAE) previously developed by R. Lamour for differential–algebraic boundary-value problems. The standard version of MSHDAE uses a BDF integrator, which is not dichotomically stable, and for some stiff test problems this fails to integrate across the interval of interest, while the dichotomically stable integrator SYMDAE encounters no difficulty. Indeed, for such problems, the modified version of MSHDAE produces an accurate solution, and within limits imposed by computer word length, the efficiency of the solution process improves with increasing stiffness. For some nonstiff problems, the solution is also entirely satisfactory.

Item Type: Journal Item
Copyright Holders: 2001 IMACS
ISSN: 1873-5460
Extra Information: Numerical Solution of Differential and Differential-Algebraic Equations, 4-9 September 2000, Halle, Germany
Keywords: differential–algebraic equations; boundary-value problems; dichotomic stability; multiple shooting
Academic Unit/School: Faculty of Science, Technology, Engineering and Mathematics (STEM) > Mathematics and Statistics
Faculty of Science, Technology, Engineering and Mathematics (STEM)
Item ID: 2403
Depositing User: Roland England
Date Deposited: 12 Jun 2006
Last Modified: 07 Dec 2018 08:51
Share this page:


Altmetrics from Altmetric

Citations from Dimensions

Actions (login may be required)

Policies | Disclaimer

© The Open University   contact the OU