  I report on an implementation of an algorithm for the automated numerical
calculation of spin- and colour-correlated Born matrix elements in QCD. These
spin- and colour-correlated matrix elements are needed for NLO calculations in
combination with the subtraction method. Both massless and massive quarks are
considered. There are no restrictions on the number of external particles. As a
trivial sub-case, the algorithm also applies to Born matrix elements without
any correlations. These are sufficient for leading order calculations.
