Invariant Subspaces, Derivative Arrays, and the Computation of the Drazin Inverse