This paper treats the slow-motion approximation for radiating systems as a problem in singular perturbations. By using the method of matched asymptotic expansions, we can construct approximations valid both in the near zone and the wave zone. The outgoing-wave boundary condition applied to the wave-zone expansion leads, by matching, to a unique and easily calculable radiation resistance in the near zone. The method is developed and illustrated with model problems from mechanics and electromagnetism; these should form a useful and accessible introduction to the method of matched asymptotic expansions. The method is then applied to the general relativistic problem of gravitational radiation from gravitationally bound systems, where a significant part of the radiation can be attributed to nonlinear terms in the expansion of the metric. This analysis shows that the formulas derived from the standard linear approximation remain valid for gravitationally bound systems. In particular, it shows that, according to general relativity, bodies in free-fall motion do indeed radiate. These results do not depend upon any definition of gravitational field energy.