Abstract
In this article, a variable step size strategy is adopted in formulating a new variable step hybrid block method (VSHBM) for the solution of the Kepler problem, which is known to be a rigid and stiff differential equation. To derive the VSHBM, the step size ratio r is left the same, halved, or doubled in order to optimize the total number of steps, minimize the number of formulae stored in the code, and ensure that the method is zero-stable. The method is formulated by integrating the Lagrange polynomial with limits of integration selected at special points. The article further analyzed the stability, order, consistency, and convergence properties of the VSHBM. The stability regions of the VSHBM at different values of the step size ratios were also plotted and plots showed that the method is fit for solving the Kepler problem. The results generated were then compared with some existing methods, including the MATLAB inbuilt stiff solver (ode 15 s), with respect to total number of failure steps, total number of steps, total function calls, maximum error, and computation time.
Subject
Statistics and Probability,Statistical and Nonlinear Physics,Analysis
Reference43 articles.
1. Variable step size selection methods for implicit integration schemes for ordinary differential equations;Holsapple;Int. J. Numer. Anal. Model.,2007
2. Numerical Solution of Ordinary Differential Equations;Atkinson,2009
3. Numerical integrators for stiff and highly oscillatory differential equations
4. Numerical Methods for Ordinary Differential Systems: The Initial Value Problem;Lambert,1991
5. A Pair of Three-Step Hybrid Block Methods for the Solutions of Linear and Nonlinear First-Order Systems
Cited by
17 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献