  Lower and upper bounds for the cumulative distribution function (cdf) of the
time to extinction in a subcritical two-sex branching process are derived. A
recursive procedure for approximating this cdf is also utilized. The results
are illustrated with some simulations.
