A Novel Approach to the Initial Value Problem with a complete validated algorithm

Abstract

We consider the first order autonomous differential equation (ODE) x'= f( x) where f: Rn Rn is locally Lipschitz. For x0∈ Rn and h>0, the initial value problem (IVP) for ( f, x0,h) is to determine if there is a unique solution, i.e., a function x:[0,h] Rn that satisfies the ODE with x(0)= x0. Write x = IVP f( x0,h) for this unique solution. We pose a corresponding computational problem, called the End Enclosure Problem: given ( f,B0,h,0) where B0⊂eq Rn is a box and 0>0, to compute a pair of non-empty boxes (B0,B1) such that B0⊂eq B0, width of B1 is <0, and for all x0∈ B0, x= IVP f( x0,h) exists and x(h)∈ B1. We provide a complete validated algorithm for this problem. Under the assumption (promise) that for all x0∈ B0, IVP f( x0,h) exists, we prove the halting of our algorithm. This is the first halting algorithm for IVP problems in such a general setting. We also introduce novel techniques for subroutines such as StepA and StepB, and a scaffold datastructure to support our End Enclosure algorithm. Among the techniques are new ways refine full- and end-enclosures based on a radical transform combined with logarithm norms. Our preliminary implementation and experiments show considerable promise, and compare well with current validated algorithms.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…