AbstractNonlinear systems often give rise to fractal boundaries in phase space, hindering predictability. When a single boundary separates three or more different basins of attraction, we say that the set of basins has theWada property and initial conditions near that boundary are even more unpredictable. Many physical systems of interest with this topological property appear in the literature. However, so far the only approach to study Wada basins has been restricted to two-dimensional phase spaces. Here we report a simple algorithm whose purpose is to look for the Wada property in a given dynamical system. Another benefit of this procedure is the possibility to classify and study intermediate situations known as partially Wada boundaries.