  Most recent calculations of quarkonium production are based on the NRQCD
factorization formalism. This formalism is reviewed. To make predictions about
specific cross section, universal NRQCD matrix elements need to be extracted
from experiments. Extractions from different experimental situations are
compared, with some emphasis on the extraction from LEP.
