A general iterative method for the Drazin inverse of a squre matrix with complex spectrum and its error bound are derived. Euler-Knopp and Newton-Raphson methods for the real spectrum case are generalized. Semiiterative methods for the Drazin inverse is also developed. Numerical examples are given to illustrate the theoretic results.