In this paper, we present a numerical method with guaranteed accuracy to solve initial value problems (IVPs) of normal form simultaneous first order ordinary differential equations (ODEs) which have wide domain. Our method is based on the algorithm proposed by Kashiwagi, by which we can obtain inclusions of exact values at several discrete points of the solution curve of ODEs. The method can be regarded as an extension of the Lohner's method. But the algorithm is not efficient for equations which have wide domain, because the error bounds become too wide from a practical point of view. Our purpose is to produce tight bounds even for such equations. We realize it by combining Kashiwagi's algorithm with the mean value form. We also consider the wrapping effects to obtain tighter bounds.
The copyright of the original papers published on this site belongs to IEICE. Unauthorized use of the original or translated papers is prohibited. See IEICE Provisions on Copyright for details.
Copy
Takatomi MIYATA, Yasutaka NAGATOMO, Masahide KASHIWAGI, "Long Time Integration for Initial Value Problems of Ordinary Differential Equations Using Power Series Arithmetic" in IEICE TRANSACTIONS on Fundamentals,
vol. E84-A, no. 9, pp. 2230-2237, September 2001, doi: .
Abstract: In this paper, we present a numerical method with guaranteed accuracy to solve initial value problems (IVPs) of normal form simultaneous first order ordinary differential equations (ODEs) which have wide domain. Our method is based on the algorithm proposed by Kashiwagi, by which we can obtain inclusions of exact values at several discrete points of the solution curve of ODEs. The method can be regarded as an extension of the Lohner's method. But the algorithm is not efficient for equations which have wide domain, because the error bounds become too wide from a practical point of view. Our purpose is to produce tight bounds even for such equations. We realize it by combining Kashiwagi's algorithm with the mean value form. We also consider the wrapping effects to obtain tighter bounds.
URL: https://global.ieice.org/en_transactions/fundamentals/10.1587/e84-a_9_2230/_p
Copy
@ARTICLE{e84-a_9_2230,
author={Takatomi MIYATA, Yasutaka NAGATOMO, Masahide KASHIWAGI, },
journal={IEICE TRANSACTIONS on Fundamentals},
title={Long Time Integration for Initial Value Problems of Ordinary Differential Equations Using Power Series Arithmetic},
year={2001},
volume={E84-A},
number={9},
pages={2230-2237},
abstract={In this paper, we present a numerical method with guaranteed accuracy to solve initial value problems (IVPs) of normal form simultaneous first order ordinary differential equations (ODEs) which have wide domain. Our method is based on the algorithm proposed by Kashiwagi, by which we can obtain inclusions of exact values at several discrete points of the solution curve of ODEs. The method can be regarded as an extension of the Lohner's method. But the algorithm is not efficient for equations which have wide domain, because the error bounds become too wide from a practical point of view. Our purpose is to produce tight bounds even for such equations. We realize it by combining Kashiwagi's algorithm with the mean value form. We also consider the wrapping effects to obtain tighter bounds.},
keywords={},
doi={},
ISSN={},
month={September},}
Copy
TY - JOUR
TI - Long Time Integration for Initial Value Problems of Ordinary Differential Equations Using Power Series Arithmetic
T2 - IEICE TRANSACTIONS on Fundamentals
SP - 2230
EP - 2237
AU - Takatomi MIYATA
AU - Yasutaka NAGATOMO
AU - Masahide KASHIWAGI
PY - 2001
DO -
JO - IEICE TRANSACTIONS on Fundamentals
SN -
VL - E84-A
IS - 9
JA - IEICE TRANSACTIONS on Fundamentals
Y1 - September 2001
AB - In this paper, we present a numerical method with guaranteed accuracy to solve initial value problems (IVPs) of normal form simultaneous first order ordinary differential equations (ODEs) which have wide domain. Our method is based on the algorithm proposed by Kashiwagi, by which we can obtain inclusions of exact values at several discrete points of the solution curve of ODEs. The method can be regarded as an extension of the Lohner's method. But the algorithm is not efficient for equations which have wide domain, because the error bounds become too wide from a practical point of view. Our purpose is to produce tight bounds even for such equations. We realize it by combining Kashiwagi's algorithm with the mean value form. We also consider the wrapping effects to obtain tighter bounds.
ER -