One-step methods and implicit extrapolation technique for index 1 differential-algebraic systems

Russian Journal of Numerical Analysis and Mathematical Modelling, 19(6) (2004), 527-553
http://dx.doi.org/10.1515/rnam.2004.19.6.527

In this paper we first formulate and prove a number of theorems concerning the convergence of combined numerical one-step methods for index 1 differential-algebraic systems. Then, we use these results to justify an implicit extrapolation technique and show their practical importance. Second, we give a theory of adjoint and symmetric one-step methods for differential-algebraic equations and we also determine symmetric methods among Runge–Kutta formulae. We prove that algebraically stable symmetric Runge–Kutta formulae are symplectic and they have a structure which is in some sense similar to the structure of Gauss methods. Finally, we come to the concept of quadratic extrapolation for index 1 differential-algebraic systems and develop an advanced version of the localglobal step size control based on the extrapolation technique. Numerical tests support the theoretical results of the paper.