Arbitrary high order A-stable and B-convergent numerical methods for ODEs via deferred correction

5 Mar 2019  ·  Saint-Cyr E. R. Koyaguerebo-Ime, Yves Bourgault ·

This paper presents a sequence of deferred correction (DC) schemes built recursively from the implicit midpoint scheme for the numerical solution of general first order ordinary differential equations (ODEs). It is proven that each scheme is A-stable, satisfies a B-convergence property, and that the correction on a scheme DC(2j) of order 2j of accuracy leads to a scheme DC2j+2 of order 2j+2. The order of accuracy is guaranteed by a deferred correction condition. Numerical experiments with standard stiff and non-stiff ODEs are performed with the DC2, ..., DC10 schemes. The results show a high accuracy of the method. The theoretical orders of accuracy are achieved together with a satisfactory stability.

PDF Abstract
No code implementations yet. Submit your code now

Categories


Numerical Analysis Numerical Analysis 65B05 \and 65L04 \and 65L05 \and 65L12 \and 65L20