Eshelby’s inclusion problem is implemented by boundary integral equation method (BIEM) and parametric dislocation dynamics (PDD), which are surface- and line-discretization approaches based on Green’s function method. In BIEM calculations, unknown constraint displacements at the inclusion boundary are deduced from surface integration associated with eigenstrain, and then unknown tractions are solved by a linear system of equations. Field quantities and strain energy are readily calculated by using the boundary tractions and displacements. In the meantime, inclusion stress by PDD is directly calculated by line integration over the dislocation elements, where the burgers vector satisfies geometrical argument of eigenstrain of inclusion. In particular, double line integral formula associated with dislocation interaction is utilized for elastic energy calculation. Computation efficiency of BIEM and PDD is mutually compared as a function of the number of Gauss integration points. Accuracy of field quantities and strain-energy of inclusion are verified by comparison of numerical results to analytical solutions.