The Drazin generalized inverse appears in a number of applications including the theory
of linear time invariant differential-algebraic equations (DAEs). In this paper we
consider its robust computation. Various methods are proposed all of them based on
the determination of bases of invariant subspaces connected with the Drazin inverse.
We also include comparisons of our methods to some of the other approaches
found in the literature.