Providing a Basin of Attraction to a Target Region by Computation of Lyapunov-like Functions

In this paper, we present a method for computing a basin of attraction to a target region for non-linear ordinary differential equations. This basin of attraction is ensured by a Lyapunov-like polynomial function that we compute using an interval based branch-and-relax algorithm. This algorithm relaxes the necessary conditions on the coefficients of the Lyapunov-like function to a system of linear interval inequalities that can then be solved exactly, and iteratively reduces the relaxation error by recursively decomposing the state space into hyper-rectangles. Tests on an implementation are promising.