Abstract
When we solve IVPs of ODEs with guaranteed accuracy, we have some problems about explosion of interval width and calculating cost, especially in the case of ODEs which is sensitive with initial value such as Lorenz equation. In this paper, we propose an accurate and fast method for IVPs using power series arithmetic, mean value form and affne arithmetic. Moreover, we show the efficiency of our method through some numerical examples and comparison with Lohner's AWA package.
