Abstract. In this note we investigate the problem of computing the domain of attraction of a flow on R 2 for a given attractor. We consider an operator that takes two inputs, the description of the flow and a cover of the attractors, and outputs the domain of attraction for the given attractor. We show that: (i) if we consider only (structurally) stable systems, the operator is (strictly semi-)computable; (ii) if we allow all systems defined by C 1 -functions, the operator is not (semi-)computable. We also address the problem of computing limit cycles on these systems.Many problems about dynamical systems (DSs) are concerned with their long term behavior. For example, given some trajectory, where will it end up? Which are the invariant sets of a DS? Which are its attractors?Recently, with the advent of increasingly powerful digital computers, numerous new ideas and concepts related to these question have appeared (e.g. sensitive dependence on initial conditions, chaos, strange attractors, Mandelbrot set). However it is interesting to note that the computer, while being an invaluable tool to get some intuition about a DS, is rarely used to prove results. Usually the formal analysis of DSs is done analytically (but often relies on information provided by numerical simulations), using heavy mathematics with little reliance on the computer. A notable exception is the proof that the Lorenz strange attractor exists and is robust under small perturbations [1], [2].One of the reasons for this phenomenon is that computers introduce truncation errors which, in conjunction with other properties such as sensitive dependency on initial conditions, is likely to produce simulated trajectories that cannot be proved accurate. Of course, there are many results exhibiting that these simulations are valuable; the foremost of such results is perhaps the Shadowing Lemma [3], [4]. However the accuracy of a particular simulation, especially when we are interested in global properties, can usually be put into question.In this paper we deal with a particular type of the problems mentioned above: is it possible to conceive a computer program that, given an input that describes a dynamical system (DS) as well as an attractor of this DS, computes (rigourously) the basin of attraction of the given attractor?