Abstract
AbstractThe Alekseev–Gröbner lemma is combined with the theory of modified equations to obtain an a priori estimate for the global error of numerical integrators. This estimate is correct up to a remainder term of order h2p, where h denotes the step size and p the order of the method. It is applied to nonlinear oscillators whose behaviour is described by the Emden–Fowler equation y″+tνyn=0. The result shows explicitly that later terms sometimes blow up faster than the leading term of order hp, necessitating the whole computation. This is supported by numerical experiments.
Subject
Computational Theory and Mathematics,General Mathematics
Reference13 articles.
1. 5 Gragg W. , ‘Repeated extrapolation to the limit in the numerical solution of ordinary differential equations’, Ph.D. thesis, University of California, Los Angeles, 1964.
2. Error growth in the numerical integration of periodic orbits by multistep methods, with application to reversible systems
Cited by
1 articles.
订阅此论文施引文献
订阅此论文施引文献，注册后可以免费订阅5篇论文的施引文献，订阅后可以查看论文全部施引文献