  Using a subtraction method, we derive the formulae suitable for use in
Monte-Carlo event generators to give the angular distribution for the
gluon-quark induced NLO corrections in Drell-Yan lepton pair production. We
also give the corresponding helicity density matrix for W and Z boson
production.
