Year: 2007
International Journal of Numerical Analysis and Modeling, Vol. 4 (2007), Iss. 2 : pp. 210–240
Abstract
Implicit integration schemes for ODEs, such as Runge-Kutta and Runge-Kutta-Nyström methods, are widely used in mathematics and engineering to numerically solve ordinary differential equations. Every integration method requires one to choose a step-size, $h$, for the integration. If $h$ is too large or too small the efficiency of an implicit scheme is relatively low. As every implicit integration scheme has a global error inherent to the scheme, we choose the total number of computations in order to achieve a prescribed global error as a measure of efficiency of the integration scheme. In this paper, we propose the idea of choosing $h$ by minimizing an efficiency function for general Runge-Kutta and Runge-Kutta-Nyström integration routines. This efficiency function is the critical component in making these methods variable step-size methods. We also investigate solving the intermediate stage values of these routines using both Newton's method and Picard iteration. We then show the efficacy of this approach on some standard problems found in the literature, including a well-known stiff system.
You do not have full access to this article.
Already a Subscriber? Sign in as an individual or via your institution
Journal Article Details
Publisher Name: Global Science Press
Language: English
DOI: https://doi.org/2007-IJNAM-860
International Journal of Numerical Analysis and Modeling, Vol. 4 (2007), Iss. 2 : pp. 210–240
Published online: 2007-01
AMS Subject Headings: Global Science Press
Copyright: COPYRIGHT: © Global Science Press
Pages: 31
Keywords: Runge-Kutta implicit integration methods variable step-size methods solving stiff systems.