  I present a highly efficient method for evolving parton distributions in
perturbative QCD. The method allows evolving the parton distribution functions
according to any of the commonly-used truncations of the evolution equations
(which differ in their treatment of higher-order terms). I also give formul\ae\
for computing crossing functions within the method.
