  We predict the shape of the transverse momentum p_T spectrum of Upsilon
production. The distribution at low p_T is dominated by the region of small
impact parameter b and may be computed reliably in perturbation theory. We
resum to all orders in the strong coupling alpha_s the process-independent
large logarithmic contributions that arise from initial-state gluon showers in
the small p_T (< M_Upsilon) region. The cross section at large p_T is
represented by the alpha_s^3 lowest-order non-vanishing perturbative
contribution.
