A particular solution of a dynamical system is completely determined by its initial condition. When the omega limit set reduces to a point, the solution settles at steady state. The possible steady states of the system are completely determined by its parameters. However, with the same parameter set, it is possible that several steady states can originate from different initial conditions (multi-stability). In that case the outcome depends on the chosen initial condition. Therefore, it is important to assess the domain of attraction for each possible attractor. The algorithms presented here are general and robust enough so as to solve the problem of reconstructing the basin of attraction of each stable equilibrium point. In order to have a graphical representation of the separatrix manifolds, we focus on systems of two and three ordinary differential equations exhibiting bi- or tri-stability. For this purpose we have implemented several Matlab functions for the approximation of the points lying on the curves or on the surfaces determining the basins of attraction and for the reconstruction of such curves and surfaces. We approximate the latter with the implicit partition of unity method using radial basis functions as local approximants. Numerical results, obtained with a Matlab package made available to the scientific community, support our findings.
['Roberto Cavoretto', 'Alessandra De Rossi', 'Emma Perracchione', 'Ezio Venturino']
Robust Approximation Algorithms for the Detection of Attraction Basins in Dynamical Systems