There is proposed a representation of the Drazin inverse of a given polynomial square matrix, based on the extension of the Leverrier--Faddeev algorithm. Also, an algorithm for symbolic computation of the Drazin inverse of polynomial matrix is established. This algorithm represents an extension of the papers [5], [7] and a continuation of the papers [8], [9], [10]. The implementation is developed in the symbolic package MATHEMATICA.