This paper addresses the problem of bounding the trajectories of nonlinear systems (transient and ultimate bounds) from initial conditions in given sets, when subject to possibly nonvanishing disturbances constrained by some finite-interval integral bounds, with a suitable controller. The so-called robustly-inescapable sets are determined from such initial conditions and disturbance bounds. In order to get numerical results, the approach considers embedding the nonlinear dynamics in a convex combination of polynomials, and solving sum-of-squares (SOS) problems on them, optimising some inescapable-set size parameters. Determination of approximate (locally) optimal solutions usually requires an iterative evaluation of SOS problems, because of products of decision variables.