Abstract. We investigate the asymptotic behavior at +∞ of non-oscillatory solutions to differential equations y = G(t, y), t > a, where G : R 1+l → R l is definable in a polynomially bounded o-minimal structure. In particular, we show that the Pfaffian closure of a polynomially bounded o-minimal structure on the real field is levelled.A classical topic in asymptotic analysis is the study of the behavior of solutions at +∞ to equations P (t, y(t), . . . , y (n) (t)) = 0, where P : R 2+n → R is a nontrivial polynomial function. It was conjectured by E. Borel that the growth of such solutions should be bounded at +∞ by the (n + 1)-st compositional iterate of e x . This turned out to be false-see Boshernitzan [2] for a detailed account-unless infinitely oscillating solutions are somehow ruled out. But with such assumptions, sharper estimates are obtained that rule out transexponential growth as well as certain kinds of intermediate asymptotic behavior (e.g., that of solutions to the functional equation g • g = e x ; see Rosenlicht [23, §2]). In this paper, we extend these results to (systems of) first-order ODEs defined over polynomially bounded o-minimal structures on the real field. Examples of large classes of maps that generate such structures are given by van den Dries and Speissegger [9, 10] and Rolin et al. [20].Some of our results hold in o-minimal expansions of arbitrary ordered fields, while some appear to make sense only over the real numbers. All but one we prove by so-called standard methods (that is, without model theory). We organize the exposition as follows: First, results are stated over R in the form of an extended abstract, followed by proofs. Then we indicate those results that, after appropriate modification, hold in o-minimal expansions of arbitrary ordered fields.