Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems
Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems
On the construction of error estimators for implicit Runge-Kutta methods
Journal of Computational and Applied Mathematics
Order reduction of stiff solvers at elastic multibody systems
Applied Numerical Mathematics - Selected papers on eighth conference on the numerical treatment of differential equations 1-5 September 1997, Alexisbad, Germany
Stiff differential equations solved by Radau methods
Proceedings of the on Numerical methods for differential equations
Journal of Computational and Applied Mathematics
| Hi-index | 7.29 |
A variable stepsize implementation of a recently introduced one-parameter family of high order strongly A-stable Runge-Kutta collocation methods with the first internal stage of explicit type is presented. The so-called SAFERK(@a,s) methods with free parameter @a and s stages are well-suited for the integration of stiff and differential-algebraic systems, and they are computationally equivalent to the (s-1)-stage Radau IIA method, since they all have a similar amount of implicitness. For the same number of implicit stages, both SAFERK(@a,s) and Radau IIA(s-1) methods possess algebraic order 2s-3, whereas the stage order is one unit higher for SAFERK methods. Although there are no L-stable schemes in this method family, the free parameter @a can be selected in order to minimize the error coefficients or to maximize the numerical dissipation. Besides a general discussion of the method class, it is shown here how the 4-stage methods can be endowed with an embedded third order formula, and an implementation based on the perfected RADAU5 code with an adaptive stepsize controller proves to be competitive for a wide selection of test problems including electric circuit analysis, constrained mechanical systems, and time-dependent partial differential equations treated by the method of lines.