For Poisson's equation over a polygonal domain of general shape, the solution of which may have a singularity around re-entrant corners, we provide an explicit a priori error estimate for the approximate solution obtained by finite element methods of high degree. The method used herein is a direct extension of the one developed in preceding paper of the second and third listed authors, which provided a new approach to deal with the singularity by using linear finite elements. In the present paper, we also give a detailed discussion of the dependency of the convergence order on solution singularities, mesh sizes and degrees of the finite element method used.